Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,718,678 C→G 38.3% G143A (GGC→GCC)  kgtP ← alpha‑ketoglutarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,718,6780CG38.3% 51.8 / 75.7 68G143A (GGC→GCC) kgtPalpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base G (11/15);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GTGCAGCGTCTTCCATGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAT  >  CP009273/2718532‑2718824
                                                                                                                                                  |                                                                                                                                                  
gTGCAGCGTCTTCCATGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCa                                                                                                                                                  <  1:590411/149‑1 (MQ=255)
 tGCAGCGTCTTCCATGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCat                                                                                                                                                 >  2:389091/1‑149 (MQ=255)
           tCCATGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGcc                                                                                                                                       >  2:284213/1‑149 (MQ=255)
             cATGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGAAACTTCACACAGATAGGTGGCGCTGGTGCCATATTCTCCGCCaa                                                                                                                                     >  2:807711/1‑149 (MQ=255)
              aTGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAAc                                                                                                                                    >  1:375760/1‑149 (MQ=255)
                  tgtgtTGTAAAACTCCGACATCCATTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGAt                                                                                                                                <  1:782749/149‑1 (MQ=255)
                    tgttgtAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATGGGTGGCGCTGGTGGCGTATTCTGCGCCAACAGATaa                                                                                                                              >  2:151018/1‑149 (MQ=255)
                    tgttgtAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATaa                                                                                                                              >  2:30447/1‑149 (MQ=255)
                             aCCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGaaa                                                                                                                     >  2:835822/1‑149 (MQ=255)
                                    cAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAAcga                                                                                                              <  2:856313/149‑1 (MQ=255)
                                    cAACCAGTAGGGCTAGCAGTTGTCAGCAGCTCAAATCCTCATACTGACATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACGGACCATCCCTGAAATAAAcga                                                                                                              >  1:189085/1‑149 (MQ=255)
                                    cAACCAGTAGGGCGAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAAcga                                                                                                              <  1:25681/149‑1 (MQ=255)
                                     aaCCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAAcgag                                                                                                             <  2:515032/149‑1 (MQ=255)
                                           tAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGaa                                                                                                       >  2:768433/1‑149 (MQ=255)
                                               gCTAGCAGTTGGCCGCCGATCAACGTTACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCaa                                                                                                   <  2:383580/149‑1 (MQ=255)
                                                cTAGCAGTTGTCCGCCGCTCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGCCGCCGCTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCaat                                                                                                  <  2:678187/149‑1 (MQ=255)
                                                cTAGCAGTTGTCCGCCGATCAACGGCACATACGGAAATGATGCGTAAACACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGATCGCGCTGGTGCCATATTCTCCGCCCACAGGAAACCCCTGAAAGAAACGAGCGAGAAGCcat                                                                                                  >  1:294444/1‑149 (MQ=255)
                                                   ccagtgccccgccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCaataat                                                                                               <  1:710546/141‑1 (MQ=255)
                                                       ttGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCg                                                                                           >  1:42848/1‑149 (MQ=255)
                                                          tCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTCCCATATACGCCGCCAACAGAAAATCCGCGAAAGAAACGGGCGAGAAGCAATAAGGCCGGaa                                                                                        >  1:94799/1‑148 (MQ=255)
                                                            cgccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGcc                                                                                      <  1:284213/149‑1 (MQ=255)
                                                            cgccgATCAAACTCACACATTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCAGATAGGTGGCGCTAGTTCCATATTCTCCGCCAACAGATCCCCCCTGAAATAACCGCCCGAGAAGCAATAATGACGGAGcc                                                                                      <  1:239447/149‑1 (MQ=255)
                                                             gccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCCACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAACCCCTGCAATAAACGAGCGCGAAGCaagaacacccgcgccc                                                                                     >  1:759821/1‑138 (MQ=255)
                                                              ccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCGGAAATAAACGAGCAAGAAGGATTGATGCCGGGACCCa                                                                                    >  2:235986/1‑149 (MQ=255)
                                                                  tCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTa                                                                                <  1:643391/149‑1 (MQ=255)
                                                                  gcAACGTCACATACTGAACTGATGCGTAAAAACCGTTGCGCCCTTCAACGGCAACTACACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTa                                                                                <  1:502946/148‑1 (MQ=255)
                                                                   cAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTAc                                                                               <  1:706818/149‑1 (MQ=255)
                                                                   cAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGAAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTAc                                                                               <  2:715670/149‑1 (MQ=255)
                                                                       gTCACATACGCGAATGATGCGTAAGAACCTTCGCGCCCTTCCACGGCAACCTCGCTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGTTAATCCCTGAAATCAACGCGCGAGAAGCAATAATGCCGGAGCCCACGTACCtat                                                                           <  1:362557/149‑1 (MQ=255)
                                                                         cacaTACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAg                                                                         >  2:359689/1‑149 (MQ=255)
                                                                           caTACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAGTGTCGGAGCCCACGTACCTATAGtt                                                                       >  2:6044/1‑149 (MQ=255)
                                                                             tACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATAGAGGTGGCGCTGGTGCCATCTTCTCCGCCAACAGAAAATCCCTGAAATAAACGAGCGACAAGCACTAAACCCGTAGCCCACCTAcgactcgagtc                                                                     >  1:848053/1‑139 (MQ=255)
                                                                                 gAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATaa                                                                 >  2:141902/1‑149 (MQ=255)
                                                                                     tgatgCGTAAAACCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGAGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAACCGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTg                                                             >  1:685647/1‑149 (MQ=255)
                                                                                     tgatgCATAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACCGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTg                                                             <  2:207529/149‑1 (MQ=255)
                                                                                        tgCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGa                                                          <  1:277363/149‑1 (MQ=255)
                                                                                              aaacCCTTTTCTCCCTTCAACGGCAACGTCACTCATATAGGTGGCGGTGGTGCCATATTCTCCGCCAACAGATAACCCCTGAAATAAACGAGCCAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTCGGaggcagg                                                    <  2:404413/149‑1 (MQ=255)
                                                                                              aaaaCCTTTGCTCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGATTCATAACCTGGGaggcagg                                                    >  1:407069/1‑149 (MQ=255)
                                                                                              aaaaCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGaggcagg                                                    >  2:616020/1‑149 (MQ=255)
                                                                                              aaaaCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCAAGTACCTATAGTTTCATGACCTGGGaggaagg                                                    >  1:372095/1‑149 (MQ=255)
                                                                                              aaaaCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCCGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGaggcagg                                                    <  2:825224/149‑1 (MQ=255)
                                                                                                   cTTTGCGCCCTTCAACGGCAACTTCACTCATATCGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATa                                               >  2:635110/1‑149 (MQ=255)
                                                                                                   cTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATa                                               >  2:651662/1‑149 (MQ=255)
                                                                                                       gcgcCCTTCAACGGCAACTTCACTCATATAGGCGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATACTGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCa                                           <  1:477134/149‑1 (MQ=255)
                                                                                                        cgcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCTGACAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                          <  2:849828/149‑1 (MQ=255)
                                                                                                        cgcCCTTCAACCGCAACTTCACTCATTTAGGCGGCGCTGGCGCCATATTCCCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                          <  2:375760/149‑1 (MQ=255)
                                                                                                         gcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAACAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGc                                         >  2:253682/1‑149 (MQ=255)
                                                                                                         gcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCGGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGc                                         <  2:54558/149‑1 (MQ=255)
                                                                                                         ccccTTCAACGGCAACTTCTCTCATATAGGTGGCTCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGc                                         <  1:318840/148‑1 (MQ=255)
                                                                                                           ccTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGGGCGCGAGGAATCCGTTCCCGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGa                                       >  2:850707/1‑149 (MQ=255)
                                                                                                           ccTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGa                                       >  2:44767/1‑149 (MQ=255)
                                                                                                              tCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATcc                                    >  1:288505/1‑149 (MQ=255)
                                                                                                              tCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGCGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATcc                                    <  1:491572/149‑1 (MQ=255)
                                                                                                                  cGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaaa                                <  2:394020/149‑1 (MQ=255)
                                                                                                                  cGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTACGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGGTCCGaaa                                >  1:408435/1‑149 (MQ=255)
                                                                                                                    gCAACTTCACTCATATAGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAaca                              <  1:332746/149‑1 (MQ=255)
                                                                                                                      aaCTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAAcaca                            <  2:77486/149‑1 (MQ=255)
                                                                                                                       aCTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATACCCCGCCAGTCAGGCGATAACCAGCGATCCGAAACAcat                           <  1:835822/149‑1 (MQ=255)
                                                                                                                        cTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGAGCCGAAACAcatc                          >  2:653883/1‑149 (MQ=255)
                                                                                                                          tCACGCATATAGGTGGCGCGGGTGGCATAGTCTCCGCCAACAGATAACCTCTGAAAGAAACGTGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATTACCAGCGATCCGAAACAcatcat                        <  1:596145/149‑1 (MQ=255)
                                                                                                                            aCTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATac                      >  2:52442/1‑149 (MQ=255)
                                                                                                                              tCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATacac                    >  2:695954/1‑149 (MQ=255)
                                                                                                                                  ataGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCg                <  2:186410/149‑1 (MQ=255)
                                                                                                                                  ataGGTGGCGCTGGTGGCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCg                <  2:42848/149‑1 (MQ=255)
                                                                                                                                        ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGGGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGAACCGAAACACATAATAGAACCCGATaaaa          >  1:281899/1‑147 (MQ=255)
                                                                                                                                        ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACa          <  2:258486/149‑1 (MQ=255)
                                                                                                                                        ggCGCTGGTGCCAGATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACa          <  1:348676/149‑1 (MQ=255)
                                                                                                                                                tGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAt  >  1:311533/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGCAGCGTCTTCCATGGTGTGTTGTAAAACCACGACAACCAGTAGGGCTAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAT  >  CP009273/2718532‑2718824

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: