Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,526,607 G→A 100% D179N (GAT→AAT)  ydcC → H repeat‑associated putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,526,6070GA98.6% 337.7 / ‑2.3 148D179N (GAT→AAT) ydcCH repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (88/57);  total (90/58)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 150 total aligned reads displayed.

CTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGTG  >  CP009273/1526518‑1526751
                                                                                         |                                                                                                                                                
cTAATGAGATTACAGCTATCCCAGACCTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGACAAATATTGCAGAGAAGATACCAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGaaa                                                                                        >  2:2251957/1‑148 (MQ=2)
cTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGaaa                                                                                        >  1:1146213/1‑148 (MQ=11)
  aaTGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAAcc                                                                                      >  1:1036623/1‑148 (MQ=17)
    tGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAg                                                                                    <  1:1597755/148‑1 (MQ=17)
     gagaTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACCGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAgg                                                                                   <  1:1229552/148‑1 (MQ=17)
     gagaTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAgg                                                                                   <  1:2250242/148‑1 (MQ=17)
     gagaTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAgg                                                                                   <  1:2248307/148‑1 (MQ=17)
      agaTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAggg                                                                                  >  2:697290/1‑148 (MQ=17)
       gaTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAgggg                                                                                 <  1:1393009/148‑1 (MQ=17)
        aTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGggc                                                                                >  2:3335819/1‑148 (MQ=32)
        aTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGggc                                                                                >  1:1724860/1‑148 (MQ=17)
         ttACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGggcg                                                                               >  1:2829335/1‑148 (MQ=17)
          tACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGggcgg                                                                              >  1:623575/1‑148 (MQ=17)
              gCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaa                                                                          >  1:2565716/1‑148 (MQ=17)
                tATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaat                                                                        <  2:61380/148‑1 (MQ=17)
                 aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaata                                                                       >  1:556217/1‑148 (MQ=17)
                 aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaata                                                                       >  1:556179/1‑148 (MQ=17)
                  tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaataa                                                                      >  1:1388877/1‑148 (MQ=17)
                  tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaataa                                                                      <  1:3313926/148‑1 (MQ=17)
                   cccAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaataaa                                                                     >  2:1354176/1‑148 (MQ=17)
                   cccAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCtaaataaa                                                                     >  2:78265/1‑148 (MQ=17)
                     cAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGc                                                                   >  1:1584281/1‑148 (MQ=17)
                     cAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGc                                                                   >  2:470842/1‑148 (MQ=17)
                      aGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGcc                                                                  >  2:2862533/1‑148 (MQ=17)
                      aGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGcc                                                                  >  1:2954567/1‑148 (MQ=17)
                       gAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCt                                                                 <  2:464540/148‑1 (MQ=17)
                       gAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCt                                                                 <  2:1371598/148‑1 (MQ=17)
                        aaCTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCtt                                                                <  2:2548593/148‑1 (MQ=17)
                        aaCTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCtt                                                                <  2:2548514/148‑1 (MQ=17)
                        aaCTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCtt                                                                >  1:765173/1‑148 (MQ=17)
                        aaCTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCtt                                                                >  1:2587042/1‑148 (MQ=17)
                          cttcttAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTg                                                              <  1:1825393/148‑1 (MQ=17)
                          cttcttAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTg                                                              <  1:3020163/148‑1 (MQ=17)
                          cttcttAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTg                                                              >  1:2101658/1‑148 (MQ=17)
                           ttcttAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGa                                                             >  1:105561/1‑148 (MQ=17)
                           ttcttAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGa                                                             >  1:1760035/1‑148 (MQ=17)
                               tAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaa                                                         >  1:3400296/1‑148 (MQ=17)
                               tAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaa                                                         >  2:1224634/1‑148 (MQ=17)
                               tAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaa                                                         >  1:493103/1‑148 (MQ=17)
                                 aCATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaa                                                       <  1:3145912/148‑1 (MQ=17)
                                  cATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaaa                                                      <  1:1844633/148‑1 (MQ=17)
                                  cATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaaa                                                      >  1:1087297/1‑148 (MQ=17)
                                  cATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaaa                                                      >  1:462824/1‑148 (MQ=17)
                                  cATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaaa                                                      >  1:419049/1‑148 (MQ=17)
                                  cATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaaa                                                      >  1:2511030/1‑148 (MQ=17)
                                     gCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAAttt                                                   <  1:1022062/148‑1 (MQ=17)
                                     gCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAAttt                                                   >  2:296610/1‑148 (MQ=17)
                                     gCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAAttt                                                   >  1:3394942/1‑148 (MQ=17)
                                     gCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGAGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAAttt                                                   <  2:2484218/148‑1 (MQ=17)
                                      cTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTc                                                  <  1:1316355/148‑1 (MQ=17)
                                        ggATATTAAAGGAAAAATCATCACCACTGATGCGAGGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCg                                                <  2:2380561/148‑1 (MQ=12)
                                         gATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGc                                               >  2:527897/1‑148 (MQ=17)
                                         gATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGc                                               <  2:2742890/148‑1 (MQ=17)
                                         gATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGc                                               <  1:57372/148‑1 (MQ=17)
                                            atTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGa                                            >  1:670655/1‑148 (MQ=17)
                                             ttAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaa                                           >  2:411102/1‑148 (MQ=17)
                                                  ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAt                                      <  2:3302041/148‑1 (MQ=17)
                                                  ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAt                                      >  2:2276271/1‑148 (MQ=17)
                                                  ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAt                                      <  2:3302031/148‑1 (MQ=17)
                                                   gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAtt                                     >  1:1905377/1‑148 (MQ=17)
                                                    aaaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTa                                    >  1:1678528/1‑148 (MQ=17)
                                                    aaaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTa                                    >  2:2491247/1‑148 (MQ=17)
                                                    aaaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTa                                    >  1:1678482/1‑148 (MQ=17)
                                                      aaaTCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTaaa                                  >  1:3931/1‑148 (MQ=17)
                                                      aaaTCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTaaa                                  >  1:572652/1‑148 (MQ=17)
                                                       aaTCATCACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaat                                 <  2:2954567/148‑1 (MQ=17)
                                                        atcatcACAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaata                                <  1:2310036/148‑1 (MQ=17)
                                                          catcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataat                              >  1:2650696/1‑148 (MQ=17)
                                                          catcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataat                              >  2:658224/1‑148 (MQ=17)
                                                           atcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATc                             <  2:1033106/148‑1 (MQ=17)
                                                            tcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                            <  2:397177/148‑1 (MQ=17)
                                                            tcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                            <  2:1970905/148‑1 (MQ=17)
                                                            tcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                            >  2:326705/1‑148 (MQ=17)
                                                            tcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                            >  1:781437/1‑148 (MQ=17)
                                                            tcaCAACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATAAAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                            >  2:326688/1‑148 (MQ=17)
                                                             cacaACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGTTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTAGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCa                           >  1:879953/1‑148 (MQ=12)
                                                             cacaACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCa                           <  2:151614/148‑1 (MQ=17)
                                                             cacaACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCa                           <  1:1910954/148‑1 (MQ=17)
                                                              acaACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAg                          >  1:1582282/1‑148 (MQ=17)
                                                              acaACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAg                          >  1:1582285/1‑148 (MQ=17)
                                                               caACTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAgc                         >  2:1772940/1‑148 (MQ=17)
                                                                aaCTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAgcg                        >  2:1061132/1‑148 (MQ=18)
                                                                aaCTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAgcg                        >  2:1060985/1‑148 (MQ=18)
                                                                  cTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCa                      <  2:128691/148‑1 (MQ=21)
                                                                  cTGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCa                      >  1:2625798/1‑148 (MQ=21)
                                                                   tGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCAt                     >  2:1429926/1‑148 (MQ=21)
                                                                   tGATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCAt                     >  2:1429890/1‑148 (MQ=21)
                                                                    gATGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATg                    <  2:169842/148‑1 (MQ=21)
                                                                      tGCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGAc                  >  2:1554613/1‑148 (MQ=21)
                                                                       gCGATGGGTTGCCAGAAAAATGTTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACa                 >  2:1267253/1‑148 (MQ=21)
                                                                       gCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACa                 <  2:3051203/148‑1 (MQ=21)
                                                                       gCGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACa                 <  2:2515659/148‑1 (MQ=21)
                                                                        cGATGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAg                <  2:2071662/148‑1 (MQ=21)
                                                                          aTGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGtt              <  2:343300/148‑1 (MQ=21)
                                                                          aTGGGTTGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGtt              <  1:1329132/148‑1 (MQ=21)
                                                                               ttGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCa         >  2:825000/1‑148 (MQ=21)
                                                                                tGCCAGAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCaa        >  1:892171/1‑148 (MQ=21)
                                                                                     gAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGt   >  1:1961727/1‑148 (MQ=37)
                                                                                     gAAAAATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGt   >  1:2332172/1‑148 (MQ=37)
                                                                                      aaaaATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGAGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGTg  >  2:30083/1‑148 (MQ=255)
                                                                                         |                                                                                                                                                
CTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGTG  >  CP009273/1526518‑1526751

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