Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,019,203 G→T 20.7% E27* (GAG→TAG)  guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,019,2030GT20.7% 636.7 / 111.4 235E27* (GAG→TAG) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (103/83);  new base T (23/26);  total (126/109)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

GGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAC  >  CP009273/3019054‑3019350
                                                                                                                                                     |                                                                                                                                                   
ggCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAag                                                                                                                                                     >  2:323105/1‑150 (MQ=255)
   tACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAgaga                                                                                                                                                  >  1:305840/1‑150 (MQ=255)
     ccGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGAtt                                                                                                                                                <  1:805828/150‑1 (MQ=255)
     ccGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGAtt                                                                                                                                                <  1:279643/150‑1 (MQ=255)
      cGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGATAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTg                                                                                                                                               >  1:1553122/1‑150 (MQ=255)
      cGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTg                                                                                                                                               <  2:3404988/150‑1 (MQ=255)
      cGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTg                                                                                                                                               >  1:1333606/1‑150 (MQ=255)
       gACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGc                                                                                                                                              >  2:1297218/1‑150 (MQ=255)
       gACAGGAAAACGTGCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGCGGCGGTTGTTTTGATGTCACCCGTAGGATCGTTCACCCGGAAGAGATTGc                                                                                                                                              >  2:4077821/1‑150 (MQ=255)
         cAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCct                                                                                                                                            >  2:4010537/1‑150 (MQ=255)
         cAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCct                                                                                                                                            >  1:2868239/1‑150 (MQ=255)
         cAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCct                                                                                                                                            >  1:464983/1‑150 (MQ=255)
          aGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCctc                                                                                                                                           >  1:2130564/1‑150 (MQ=255)
          aGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctc                                                                                                                                           >  1:4530855/1‑150 (MQ=255)
              aaaCGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcg                                                                                                                                       <  2:4036771/150‑1 (MQ=255)
                 cGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctg                                                                                                                                    <  1:333862/150‑1 (MQ=255)
                  gTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctgc                                                                                                                                   >  1:3706915/1‑150 (MQ=255)
                   ttCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctgcg                                                                                                                                  <  1:2226293/150‑1 (MQ=255)
                   ttCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctgcg                                                                                                                                  <  2:463092/150‑1 (MQ=255)
                     cTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGt                                                                                                                                <  2:1916583/150‑1 (MQ=255)
                      tGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGtt                                                                                                                               >  1:2178497/1‑150 (MQ=255)
                      tGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGtt                                                                                                                               >  1:3390929/1‑150 (MQ=255)
                       gCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGttt                                                                                                                              <  2:2935217/150‑1 (MQ=255)
                              aTTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAgg                                                                                                                       >  2:2223738/1‑150 (MQ=255)
                                     cAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGttt                                                                                                                <  1:1474393/150‑1 (MQ=255)
                                       gCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttat                                                                                                              <  2:2354869/150‑1 (MQ=255)
                                       gCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttat                                                                                                              <  2:210235/150‑1 (MQ=255)
                                       gCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttat                                                                                                              <  2:3685372/150‑1 (MQ=255)
                                         ggAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                            >  2:3581430/1‑150 (MQ=255)
                                          gAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTAc                                                                                                           <  1:1716170/150‑1 (MQ=255)
                                           aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACt                                                                                                          >  2:3001131/1‑150 (MQ=255)
                                              agaTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCAt                                                                                                       <  2:3542844/150‑1 (MQ=255)
                                                  ggATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTaaa                                                                                                   >  2:1191595/1‑150 (MQ=255)
                                                   gATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAAc                                                                                                  <  1:4345420/150‑1 (MQ=255)
                                                    aTTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACa                                                                                                 >  2:3881408/1‑150 (MQ=255)
                                                       aaCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAggg                                                                                              >  1:3089145/1‑150 (MQ=255)
                                                        aCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTCTTGAGGATGGTTTATTACTCATTAAACAGGGa                                                                                             >  1:889237/1‑150 (MQ=255)
                                                        aCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGa                                                                                             >  1:889138/1‑150 (MQ=255)
                                                             aacggAGTCAATGATGTCAGGAGAACACACGTTATAAGCGGTACGAGGCTGTTTTATTGATGTCACCCGTACGATCGATAACGCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGt                                                                                        <  2:4428298/147‑1 (MQ=255)
                                                             aaaGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGt                                                                                        >  1:1240864/1‑150 (MQ=255)
                                                              aaGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTg                                                                                       <  1:1003386/150‑1 (MQ=255)
                                                              aaGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTg                                                                                       >  1:3760207/1‑150 (MQ=255)
                                                               aGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTgg                                                                                      >  2:668090/1‑150 (MQ=255)
                                                                   gTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATg                                                                                  >  1:2121403/1‑150 (MQ=255)
                                                                        tgatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                             <  2:582914/150‑1 (MQ=255)
                                                                         gatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTgg                                                                            <  2:934437/150‑1 (MQ=255)
                                                                         gatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTgg                                                                            <  2:1748714/150‑1 (MQ=255)
                                                                           tgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCg                                                                          <  1:2584592/150‑1 (MQ=255)
                                                                             tCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGaa                                                                        >  2:3380651/1‑150 (MQ=255)
                                                                             tCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGaa                                                                        <  2:2546091/150‑1 (MQ=255)
                                                                              cAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAAt                                                                       <  2:412777/150‑1 (MQ=255)
                                                                              cAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAAt                                                                       <  1:999998/150‑1 (MQ=255)
                                                                                 gagaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATggg                                                                    >  1:1958354/1‑150 (MQ=255)
                                                                                   gaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAAAAGGGAAAAGTGGAATGGTTTTGCGAATGGTaa                                                                  >  2:3309383/1‑150 (MQ=255)
                                                                                      cacacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAAc                                                               <  2:2522070/150‑1 (MQ=255)
                                                                                       acacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACg                                                              >  2:1660166/1‑150 (MQ=255)
                                                                                       acacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACg                                                              >  2:1921790/1‑150 (MQ=255)
                                                                                       acacGTTAAAAGCGGTACGAGGCAGTTTTAGTGATGTCACCCGTACGAACGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTAATGAGGATGGTTTATTACTCATAAAACAGGGAAAAGTGGAATGGTTTTGCGTATGGGAAAACg                                                              >  2:3921489/1‑150 (MQ=255)
                                                                                        cacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACgg                                                             <  1:51692/150‑1 (MQ=255)
                                                                                           gTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGaaa                                                          >  2:2601561/1‑150 (MQ=255)
                                                                                            ttAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAg                                                         >  2:1467898/1‑150 (MQ=255)
                                                                                              aaaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCa                                                       >  2:3462234/1‑150 (MQ=255)
                                                                                              aaaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCa                                                       >  2:2346031/1‑150 (MQ=255)
                                                                                               aaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCAt                                                      >  1:2678233/1‑150 (MQ=255)
                                                                                               aaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCAt                                                      >  1:729065/1‑150 (MQ=255)
                                                                                                 aGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCa                                                    >  1:2795059/1‑150 (MQ=255)
                                                                                                   cGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTGTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCaaa                                                  <  2:1862126/150‑1 (MQ=255)
                                                                                                   cGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCaaa                                                  >  1:3538856/1‑150 (MQ=255)
                                                                                                       aCGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTcc                                              <  1:59066/150‑1 (MQ=255)
                                                                                                       aCGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTcc                                              <  1:2359091/150‑1 (MQ=255)
                                                                                                         gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            >  2:2146681/1‑150 (MQ=255)
                                                                                                         gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATAGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            <  1:2744290/150‑1 (MQ=255)
                                                                                                         gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGAGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            >  2:2238097/1‑150 (MQ=255)
                                                                                                         gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGGTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            >  2:161077/1‑150 (MQ=255)
                                                                                                          aGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGa                                           >  2:744649/1‑150 (MQ=255)
                                                                                                          aGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGa                                           >  1:1094431/1‑150 (MQ=255)
                                                                                                          aGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGa                                           <  2:4464234/150‑1 (MQ=255)
                                                                                                           ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGac                                          >  1:2944553/1‑150 (MQ=255)
                                                                                                            gCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGaca                                         <  2:2888327/150‑1 (MQ=255)
                                                                                                              aGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACAcc                                       >  2:659811/1‑150 (MQ=255)
                                                                                                                  ttATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTc                                   <  2:1553122/150‑1 (MQ=255)
                                                                                                                     ttGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTcgcg                                >  2:2747174/1‑150 (MQ=255)
                                                                                                                      tGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGt                               >  2:1863063/1‑150 (MQ=255)
                                                                                                                      tGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGt                               >  2:564495/1‑150 (MQ=255)
                                                                                                                               ccGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTa                      <  2:1106780/150‑1 (MQ=255)
                                                                                                                                cGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTAt                     <  1:402178/150‑1 (MQ=255)
                                                                                                                                 gTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATc                    <  2:1806212/150‑1 (MQ=255)
                                                                                                                                 gTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATc                    <  2:2130286/150‑1 (MQ=255)
                                                                                                                                   aCGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgc                  <  1:3292200/150‑1 (MQ=255)
                                                                                                                                   aCGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgc                  <  1:1059523/150‑1 (MQ=255)
                                                                                                                                    cgatcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgcg                 <  2:2487253/150‑1 (MQ=255)
                                                                                                                                      atcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGc               >  1:2907969/1‑150 (MQ=255)
                                                                                                                                        cgatAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaa             <  2:3736034/150‑1 (MQ=255)
                                                                                                                                           tAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTTCTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGTGCGACTATCGCGGCAAACt          <  2:2178497/150‑1 (MQ=255)
                                                                                                                                               ccGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATa      <  1:3581430/150‑1 (MQ=255)
                                                                                                                                               ccGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATa      <  2:4530855/150‑1 (MQ=255)
                                                                                                                                                  gAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGACAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTa   >  1:1950585/1‑150 (MQ=255)
                                                                                                                                                  gAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTa   >  1:4382540/1‑150 (MQ=255)
                                                                                                                                                   aaGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTTACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAc  >  2:602456/1‑150 (MQ=255)
                                                                                                                                                   aaGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAc  >  2:4253940/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAC  >  CP009273/3019054‑3019350

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: ‑