Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,464,061 G→T 7.2% E41* (GAA→TAA)  bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,464,0610GT 7.2% 1003.7 / 34.3 263E41* (GAA→TAA) bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base G (122/122);  new base T (14/5);  total (136/127)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 267 total aligned reads displayed.

GTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGC  >  CP009273/4463912‑4464210
                                                                                                                                                     |                                                                                                                                                     
 ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGa                                                                                                                                                      >  1:3624843/1‑150 (MQ=255)
   ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAg                                                                                                                                                    <  2:1714602/150‑1 (MQ=255)
    tAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                                                                                                                                                   >  2:472730/1‑150 (MQ=255)
     aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTCCAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGtt                                                                                                                                                  <  2:3255149/150‑1 (MQ=255)
     aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTTAAGtt                                                                                                                                                  >  2:162359/1‑150 (MQ=255)
     aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTTAAGtt                                                                                                                                                  >  2:3707631/1‑150 (MQ=255)
     aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTTAAGtt                                                                                                                                                  >  2:3072116/1‑150 (MQ=255)
     aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGtt                                                                                                                                                  >  1:34603/1‑150 (MQ=255)
     aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGtt                                                                                                                                                  <  2:134196/150‑1 (MQ=255)
      aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                                                                                                                                                 <  2:3393410/150‑1 (MQ=255)
           cTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt                                                                                                                                            <  2:1196794/150‑1 (MQ=255)
           cTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt                                                                                                                                            >  1:1663291/1‑150 (MQ=255)
           cTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt                                                                                                                                            >  1:1284864/1‑150 (MQ=255)
              tCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTAt                                                                                                                                         >  2:3647639/1‑150 (MQ=255)
              tCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTAt                                                                                                                                         >  2:1583974/1‑150 (MQ=255)
               cAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATc                                                                                                                                        >  1:4407343/1‑150 (MQ=255)
                  cAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTg                                                                                                                                     >  1:1141617/1‑150 (MQ=255)
                  cAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTg                                                                                                                                     >  2:1601040/1‑150 (MQ=255)
                    ggTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt                                                                                                                                   <  2:356748/150‑1 (MQ=255)
                    ggTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt                                                                                                                                   <  2:4013650/150‑1 (MQ=255)
                    ggTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt                                                                                                                                   <  1:2904725/150‑1 (MQ=255)
                     gTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTTAAGTTACTGATTATCTTGGTa                                                                                                                                  >  2:4050316/1‑150 (MQ=255)
                     gTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTTAAGTTACTGATTATCTTGGTa                                                                                                                                  >  2:1282982/1‑150 (MQ=255)
                       tttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTAtt                                                                                                                                >  2:4462106/1‑150 (MQ=255)
                        ttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTa                                                                                                                               <  1:2053890/150‑1 (MQ=255)
                         tGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTaa                                                                                                                              >  2:2150992/1‑150 (MQ=255)
                          gTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAc                                                                                                                             >  2:245141/1‑150 (MQ=255)
                          gTAATGGTCACTAAAAAACAATATCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATATTGGTATTAAc                                                                                                                             >  1:2967510/1‑150 (MQ=255)
                            aaTGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAccc                                                                                                                           <  1:728111/150‑1 (MQ=255)
                              tGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAccccc                                                                                                                         <  1:972050/150‑1 (MQ=255)
                              tGGACACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCACTGCGGCAAAAGTGCTATTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAGTCTCTTGGTATTAAccccc                                                                                                                         <  2:207577/150‑1 (MQ=255)
                                gTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCg                                                                                                                       <  2:2498446/150‑1 (MQ=255)
                                  cACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAg                                                                                                                     >  2:559083/1‑150 (MQ=255)
                                  cACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAg                                                                                                                     <  2:2396983/150‑1 (MQ=255)
                                       aaaaaCAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCctct                                                                                                                >  1:144540/1‑150 (MQ=255)
                                         aaaCAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTAc                                                                                                              >  2:1873568/1‑150 (MQ=255)
                                         aaaCAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTAc                                                                                                              >  1:878420/1‑150 (MQ=255)
                                          aaCAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTAcg                                                                                                             >  1:3246463/1‑150 (MQ=255)
                                             aaTCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCgg                                                                                                          >  2:3965984/1‑150 (MQ=255)
                                                 tcGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCttt                                                                                                      <  1:3909548/150‑1 (MQ=255)
                                                  cgcgTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCtttt                                                                                                     <  2:3522823/150‑1 (MQ=255)
                                                  cgcgTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCtttt                                                                                                     <  1:2002144/150‑1 (MQ=255)
                                                   gcgTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTg                                                                                                    <  1:3914803/150‑1 (MQ=255)
                                                    cgTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTgg                                                                                                   >  2:2662934/1‑150 (MQ=255)
                                                      ttCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCa                                                                                                 <  1:1638576/150‑1 (MQ=255)
                                                          aGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTaa                                                                                             >  2:3876615/1‑150 (MQ=255)
                                                          aGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTaa                                                                                             >  2:772237/1‑150 (MQ=255)
                                                          aGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTaa                                                                                             >  2:1312460/1‑150 (MQ=255)
                                                           ggTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTaaa                                                                                            <  2:3298735/150‑1 (MQ=255)
                                                              cgtcCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCt                                                                                         <  1:1837531/150‑1 (MQ=255)
                                                               gtcCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTg                                                                                        <  1:1735244/150‑1 (MQ=255)
                                                               gtcCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTg                                                                                        <  1:2098177/150‑1 (MQ=255)
                                                                 cccAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTggg                                                                                      <  1:3602120/150‑1 (MQ=255)
                                                                   cAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGtt                                                                                    >  1:4384650/1‑150 (MQ=255)
                                                                          ttCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGc                                                                             <  1:2384882/150‑1 (MQ=255)
                                                                            cGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCg                                                                           >  2:2967780/1‑150 (MQ=255)
                                                                             gCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCgt                                                                          >  1:4311410/1‑150 (MQ=255)
                                                                                ccTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTa                                                                       >  2:750805/1‑150 (MQ=255)
                                                                                   gAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTAATc                                                                    >  2:812057/1‑150 (MQ=255)
                                                                                     gCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCaa                                                                  >  2:4505148/1‑150 (MQ=255)
                                                                                           aaTCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATa                                                            >  2:1783721/1‑150 (MQ=255)
                                                                                           aaTCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATa                                                            >  1:2492331/1‑150 (MQ=255)
                                                                                            aTCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATAc                                                           >  1:4305258/1‑150 (MQ=255)
                                                                                            aTCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATAc                                                           >  1:32742/1‑150 (MQ=255)
                                                                                            aTCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATAc                                                           >  1:1384428/1‑150 (MQ=255)
                                                                                                 tGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCgg                                                      <  1:1567593/150‑1 (MQ=255)
                                                                                                   cGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTa                                                    <  2:3857493/150‑1 (MQ=255)
                                                                                                    ggCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTAc                                                   <  2:226081/150‑1 (MQ=255)
                                                                                                     gCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACg                                                  >  1:2080871/1‑150 (MQ=255)
                                                                                                      cAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTTAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACgg                                                 >  2:2403227/1‑150 (MQ=255)
                                                                                                        aaaGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGaa                                               <  1:1601040/150‑1 (MQ=255)
                                                                                                              cTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTAtt                                         <  2:456905/150‑1 (MQ=255)
                                                                                                               tttttCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTc                                        >  2:630139/1‑150 (MQ=255)
                                                                                                                 tttCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCg                                      <  1:2827660/150‑1 (MQ=255)
                                                                                                                 tttCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCg                                      >  2:2049265/1‑150 (MQ=255)
                                                                                                                    cACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCtt                                   <  1:2042891/150‑1 (MQ=255)
                                                                                                                       cAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGcc                                >  1:4041763/1‑150 (MQ=255)
                                                                                                                       cAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGcc                                <  1:785862/150‑1 (MQ=255)
                                                                                                                       cAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGcc                                >  1:4552550/1‑150 (MQ=255)
                                                                                                                        aaaaaGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCg                               >  1:3480953/1‑150 (MQ=255)
                                                                                                                          aaaGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGTAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGat                             <  1:338983/150‑1 (MQ=255)
                                                                                                                             ggTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATAtt                          >  1:2899417/1‑150 (MQ=255)
                                                                                                                              gTTTCGATGCTGTCAGTGTTGCTTAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATAttc                         >  2:4498370/1‑150 (MQ=255)
                                                                                                                              gTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATAttc                         >  1:2271936/1‑150 (MQ=255)
                                                                                                                               tttCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATAttct                        >  2:4329571/1‑150 (MQ=255)
                                                                                                                                ttCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATAttctt                       <  2:1652685/150‑1 (MQ=255)
                                                                                                                                 tCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATAttcttc                      >  2:3639020/1‑150 (MQ=255)
                                                                                                                                    aTGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGtg                   <  1:1050843/150‑1 (MQ=255)
                                                                                                                                       cTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGtgatg                <  1:832360/150‑1 (MQ=255)
                                                                                                                                        tGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGtgatga               >  2:1485030/1‑150 (MQ=255)
                                                                                                                                         gTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTgatgat              <  1:559083/150‑1 (MQ=255)
                                                                                                                                            aGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGt           >  1:4504505/1‑150 (MQ=255)
                                                                                                                                               gtTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCAGATATTCTTCGTGATGATCGTCCg        >  2:2777418/1‑150 (MQ=255)
                                                                                                                                                 tGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGt      <  1:3813702/150‑1 (MQ=255)
                                                                                                                                                 tGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGt      <  2:2644108/150‑1 (MQ=255)
                                                                                                                                                  gCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTa     >  1:2139936/1‑150 (MQ=255)
                                                                                                                                                   cTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAg    <  2:3204209/150‑1 (MQ=255)
                                                                                                                                                   cTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAg    <  2:1139610/150‑1 (MQ=255)
                                                                                                                                                   cTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAg    <  2:2513407/150‑1 (MQ=255)
                                                                                                                                                     gAAGTTAATGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGc  <  1:1508285/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGC  >  CP009273/4463912‑4464210

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: