Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,901,514 C→A 38.2% P504P (CCG→CCT pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,901,5140CA38.2% 294.3 / 215.3 227P504P (CCG→CCTpyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (70/70);  new base A (45/42);  total (115/112)
Fisher's exact test for biased strand distribution p-value = 8.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

GCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGG  >  CP009273/2901369‑2901661
                                                                                                                                                 |                                                                                                                                                   
gCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGa                                                                                                                                                  >  2:59139/1‑149 (MQ=255)
  cGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAAc                                                                                                                                                <  1:3713890/149‑1 (MQ=255)
   gCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCAGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACt                                                                                                                                               >  1:683022/1‑149 (MQ=255)
    ctctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACtt                                                                                                                                              <  1:103484/149‑1 (MQ=255)
    ctctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACtt                                                                                                                                              >  2:3534885/1‑149 (MQ=255)
    ctctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACtt                                                                                                                                              >  1:4010101/1‑149 (MQ=255)
     tctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTc                                                                                                                                             <  1:1180224/149‑1 (MQ=255)
     tctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGCCACGCCACGAACCACGGGTGATTCGGAACTTc                                                                                                                                             >  2:3112218/1‑149 (MQ=255)
           ttttttACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCgatgat                                                                                                                                       >  1:108245/1‑149 (MQ=255)
               ttacttacTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCg                                                                                                                                   <  1:1560837/149‑1 (MQ=255)
                tacttactTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGa                                                                                                                                  <  1:238837/149‑1 (MQ=255)
                 acttactgCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGAc                                                                                                                                 <  1:2988562/149‑1 (MQ=255)
                   ttacttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCa                                                                                                                               >  2:3434636/1‑149 (MQ=255)
                    tacttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCaa                                                                                                                              >  2:3648515/1‑149 (MQ=255)
                     acttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAAc                                                                                                                             >  2:1237871/1‑149 (MQ=255)
                       ttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTg                                                                                                                           >  1:2845211/1‑149 (MQ=255)
                         cGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGat                                                                                                                         <  2:3818323/149‑1 (MQ=255)
                         cGCCTGACGTGTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGat                                                                                                                         <  2:782092/149‑1 (MQ=255)
                             tGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGatcatc                                                                                                                     >  2:1659886/1‑149 (MQ=255)
                              gACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATcc                                                                                                                    >  2:234759/1‑149 (MQ=255)
                                   ttCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGa                                                                                                               <  1:2151197/149‑1 (MQ=255)
                                   ttCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGa                                                                                                               >  2:29251/1‑149 (MQ=255)
                                       ggAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACgg                                                                                                           >  2:2205031/1‑149 (MQ=255)
                                            tCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGc                                                                                                      >  1:3792798/1‑149 (MQ=255)
                                            tCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGAGGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGc                                                                                                      <  1:743890/149‑1 (MQ=255)
                                             cGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCa                                                                                                     <  1:1948912/149‑1 (MQ=255)
                                                 ggcggcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAAcgc                                                                                                 <  1:1790531/149‑1 (MQ=255)
                                                   cggcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACgcgc                                                                                               <  1:2599502/149‑1 (MQ=255)
                                                   cggcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCGCGACCAACTGATCATCCCCGGAACGGCCCGCAACgcgc                                                                                               <  1:2173820/149‑1 (MQ=255)
                                                    ggcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCa                                                                                              >  1:3190637/1‑149 (MQ=255)
                                                     gcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAg                                                                                             >  1:216972/1‑149 (MQ=255)
                                                      cTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGa                                                                                            >  1:2046685/1‑149 (MQ=255)
                                                           acaAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGc                                                                                       <  1:2250505/149‑1 (MQ=255)
                                                             aaaGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCAt                                                                                     >  2:3743547/1‑149 (MQ=255)
                                                               aGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCt                                                                                   >  2:3254209/1‑149 (MQ=255)
                                                               aGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCt                                                                                   >  2:400640/1‑149 (MQ=255)
                                                                gCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCtt                                                                                  <  1:334348/149‑1 (MQ=255)
                                                                gCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCtt                                                                                  <  1:333891/149‑1 (MQ=255)
                                                                   tGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCaa                                                                               <  1:1645541/149‑1 (MQ=255)
                                                                   tGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCaa                                                                               >  2:1175020/1‑149 (MQ=255)
                                                                   tGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCaa                                                                               <  1:1783144/149‑1 (MQ=255)
                                                                    gCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAAt                                                                              <  2:3883533/149‑1 (MQ=255)
                                                                    gCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAAt                                                                              <  2:529156/149‑1 (MQ=255)
                                                                      aaaCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCt                                                                            <  1:3885520/149‑1 (MQ=255)
                                                                      aaaCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCt                                                                            <  1:3885600/149‑1 (MQ=255)
                                                                      aaaCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCt                                                                            >  2:902330/1‑149 (MQ=255)
                                                                       aaCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTg                                                                           >  2:2815233/1‑149 (MQ=255)
                                                                       aaCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTg                                                                           >  1:1429052/1‑149 (MQ=255)
                                                                        aCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGt                                                                          <  1:354094/149‑1 (MQ=255)
                                                                         cggCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGtt                                                                         <  1:2130581/147‑1 (MQ=255)
                                                                         ccgCGGGTGACCATCCCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGtt                                                                         <  1:51810/147‑1 (MQ=255)
                                                                           gCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTc                                                                       >  2:3833482/1‑149 (MQ=255)
                                                                                  aCCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCa                                                                >  1:980968/1‑149 (MQ=255)
                                                                                    cATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtg                                                              <  2:1199893/149‑1 (MQ=255)
                                                                                     aTCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCGAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgt                                                             <  1:3061503/149‑1 (MQ=255)
                                                                                     aTCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgt                                                             >  2:1259043/1‑149 (MQ=255)
                                                                                     aTCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgt                                                             <  2:3601469/149‑1 (MQ=255)
                                                                                      tCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgtt                                                            >  2:2380342/1‑149 (MQ=255)
                                                                                      tCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgtt                                                            >  1:1933603/1‑149 (MQ=255)
                                                                                        aCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttgt                                                          >  1:3280250/1‑149 (MQ=255)
                                                                                        aCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttgt                                                          >  1:1104256/1‑149 (MQ=255)
                                                                                        aCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttgt                                                          <  2:452336/149‑1 (MQ=255)
                                                                                        aCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCACAGACCTGCATCTTCAATCTGTTTCACCAGCCtgttgt                                                          >  2:3100601/1‑149 (MQ=255)
                                                                                          gTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTACGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGttgttg                                                        >  2:2458688/1‑149 (MQ=255)
                                                                                           tGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGa                                                       <  2:266195/149‑1 (MQ=255)
                                                                                              aGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACtt                                                    <  1:3434636/149‑1 (MQ=255)
                                                                                               gTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTc                                                   <  2:683022/149‑1 (MQ=255)
                                                                                                   aaGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                               <  2:2046685/149‑1 (MQ=255)
                                                                                                     gTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACg                                             >  2:1077121/1‑149 (MQ=255)
                                                                                                       aaaCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCACCTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGt                                           <  2:2647905/149‑1 (MQ=255)
                                                                                                        aaCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTg                                          <  2:3679878/149‑1 (MQ=255)
                                                                                                         aCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGa                                         >  1:3563974/1‑149 (MQ=255)
                                                                                                         aCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGa                                         <  1:635160/149‑1 (MQ=255)
                                                                                                          cTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGAc                                        <  1:275386/149‑1 (MQ=255)
                                                                                                          cTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGAc                                        >  1:2304064/1‑149 (MQ=255)
                                                                                                           tCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACg                                       >  2:321542/1‑149 (MQ=255)
                                                                                                            ccGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGa                                      >  2:1472928/1‑149 (MQ=255)
                                                                                                             cGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGAt                                     <  2:2236836/149‑1 (MQ=255)
                                                                                                              ggATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATg                                    >  1:3162404/1‑149 (MQ=255)
                                                                                                               gATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                   <  1:3743547/149‑1 (MQ=255)
                                                                                                                aTGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGAc                                  <  2:1866717/149‑1 (MQ=255)
                                                                                                                  ggAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGc                                >  1:994717/1‑149 (MQ=255)
                                                                                                                   gAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCt                               <  1:2479/149‑1 (MQ=255)
                                                                                                                   gAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCt                               <  1:2532/149‑1 (MQ=255)
                                                                                                                    aaCGGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTc                              <  2:2807684/149‑1 (MQ=255)
                                                                                                                      cTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGCCCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaa                            >  2:914637/1‑149 (MQ=255)
                                                                                                                           aaGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAtt                       >  1:2930969/1‑149 (MQ=255)
                                                                                                                             gCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACGTATCAGCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACTGTGACGATGACGCTCAACAATTGt                     <  2:1391877/149‑1 (MQ=255)
                                                                                                                                aCGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCgg                  >  2:1550461/1‑149 (MQ=255)
                                                                                                                                    aCCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCa              >  2:2857958/1‑149 (MQ=255)
                                                                                                                                     ccACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt             >  1:2928104/1‑149 (MQ=255)
                                                                                                                                      cACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAtt            >  1:9477/1‑149 (MQ=255)
                                                                                                                                       aCGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg           >  1:1514526/1‑149 (MQ=255)
                                                                                                                                          ggTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTAc        <  1:2791126/149‑1 (MQ=255)
                                                                                                                                           gTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACa       <  2:2319109/149‑1 (MQ=255)
                                                                                                                                             gATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGc     <  1:2254422/149‑1 (MQ=255)
                                                                                                                                               ttCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTg   <  2:1933900/149‑1 (MQ=255)
                                                                                                                                                tCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg  <  1:3659832/149‑1 (MQ=255)
                                                                                                                                                tCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg  <  1:3534885/149‑1 (MQ=255)
                                                                                                                                                tCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg  >  1:3117835/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGG  >  CP009273/2901369‑2901661

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