Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,704,060 A→T 6.6% Q221L (CAA→CTA)  rsxD → SoxR iron‑sulfur cluster reduction factor component; putative membrane protein of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,704,0600AT 6.6% 1164.0 / 54.8 317Q221L (CAA→CTA) rsxDSoxR iron‑sulfur cluster reduction factor component; putative membrane protein of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (150/145);  new base T (15/6);  total (166/151)
Fisher's exact test for biased strand distribution p-value = 7.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 324 total aligned reads displayed.

TGATATGAACACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGC  >  CP009273/1703911‑1704209
                                                                                                                                                     |                                                                                                                                                     
tgATAAGAACACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAGCCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCa                                                                                                                                                       <  1:2950333/150‑1 (MQ=255)
   tatGAACACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATg                                                                                                                                                    <  1:1118276/150‑1 (MQ=255)
     tGAACACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCAGATCTACAGCGGTATTCTGGAGGGCGCTGGTTGGCAATggg                                                                                                                                                  >  2:1644416/1‑150 (MQ=255)
      gAACACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGt                                                                                                                                                 >  2:3150701/1‑150 (MQ=255)
       aaCACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTa                                                                                                                                                >  1:4665475/1‑150 (MQ=255)
           cacTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAAtc                                                                                                                                            <  2:1496367/150‑1 (MQ=255)
            acTACGCTGAGGAATTGCTGGCATTAGTCAGGCGACACCGCTGGATACAATTACAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAAtct                                                                                                                                           <  2:3355349/150‑1 (MQ=255)
                 gCTTAGGTATTTATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCAGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGTGTAAATCTCTCct                                                                                                                                      <  1:5063742/150‑1 (MQ=255)
                 gCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCct                                                                                                                                      <  1:4315608/150‑1 (MQ=255)
                 gCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCct                                                                                                                                      >  2:5317093/1‑150 (MQ=255)
                  cTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctg                                                                                                                                     >  1:1055085/1‑150 (MQ=255)
                  cTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctg                                                                                                                                     >  1:3792601/1‑150 (MQ=255)
                   ttAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctgg                                                                                                                                    <  2:2492928/150‑1 (MQ=255)
                    tAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCctggc                                                                                                                                   <  2:1156635/150‑1 (MQ=255)
                    tAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctggc                                                                                                                                   <  1:1284653/150‑1 (MQ=255)
                     aGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctggct                                                                                                                                  <  1:5334118/150‑1 (MQ=255)
                      ggTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctggctg                                                                                                                                 <  1:2723478/150‑1 (MQ=255)
                       gTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctggctgg                                                                                                                                <  2:3051604/150‑1 (MQ=255)
                       gTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctggctgg                                                                                                                                >  1:1763959/1‑150 (MQ=255)
                         aTTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCctggctggct                                                                                                                              <  1:3176845/150‑1 (MQ=255)
                             aTGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTggcg                                                                                                                          >  1:4711010/1‑150 (MQ=255)
                             aTGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTggcg                                                                                                                          >  2:759700/1‑150 (MQ=255)
                               ggCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTggcggc                                                                                                                        >  2:5264028/1‑150 (MQ=255)
                               ggCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTggcggc                                                                                                                        >  1:1961098/1‑150 (MQ=255)
                               ggCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCAGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCATGTCTGGCTggcggc                                                                                                                        >  2:2554362/1‑150 (MQ=255)
                                gCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGgcggcg                                                                                                                       <  1:1688177/150‑1 (MQ=255)
                                   ttAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCGGCGGCGTATg                                                                                                                   <  2:1886529/150‑1 (MQ=255)
                                      gTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGt                                                                                                                 >  2:515707/1‑150 (MQ=255)
                                       tCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGtt                                                                                                                <  2:3930284/150‑1 (MQ=255)
                                              acacCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATgg                                                                                                         <  2:3484994/150‑1 (MQ=255)
                                                acCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTCTGGTTGCTATGGCa                                                                                                       >  2:2833404/1‑150 (MQ=255)
                                                acCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCa                                                                                                       >  1:1491693/1‑150 (MQ=255)
                                                acCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCa                                                                                                       >  1:1198138/1‑150 (MQ=255)
                                                 ccGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAg                                                                                                      >  1:3643943/1‑150 (MQ=255)
                                                     tGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAg                                                                                                  <  1:981604/150‑1 (MQ=255)
                                                       gATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCg                                                                                                >  1:3943629/1‑150 (MQ=255)
                                                        aTACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGa                                                                                               >  1:4520609/1‑150 (MQ=255)
                                                        aTACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGa                                                                                               >  2:1979207/1‑150 (MQ=255)
                                                          aCATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGAtt                                                                                             >  1:5425202/1‑150 (MQ=255)
                                                            aTTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTct                                                                                           >  1:2480641/1‑149 (MQ=255)
                                                                  aaCCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCa                                                                                     >  2:3100801/1‑150 (MQ=255)
                                                                  aaCCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCa                                                                                     >  2:1210717/1‑150 (MQ=255)
                                                                    ccTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCata                                                                                   <  2:4584191/150‑1 (MQ=255)
                                                                     ctctGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCatat                                                                                  <  2:892646/150‑1 (MQ=255)
                                                                     ctctGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCatat                                                                                  >  2:2303534/1‑150 (MQ=255)
                                                                      tctGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATAtt                                                                                 <  2:448873/150‑1 (MQ=255)
                                                                      tctGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATAtt                                                                                 <  1:3013644/150‑1 (MQ=255)
                                                                        tGTCCGTGCCGGTCATTCGGTTGAACAGATGATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTcc                                                                               <  2:3728412/150‑1 (MQ=255)
                                                                          tCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTcccc                                                                             >  1:2363602/1‑150 (MQ=255)
                                                                             gTGCCGGTCATTCGGTTGAACCGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCa                                                                          <  2:414599/150‑1 (MQ=255)
                                                                             gTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCa                                                                          >  2:4088008/1‑150 (MQ=255)
                                                                             gTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCa                                                                          >  2:1063707/1‑150 (MQ=255)
                                                                              tGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAg                                                                         <  2:5025332/150‑1 (MQ=255)
                                                                                ccGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGAGTCGCTGGCATATTCCCCTCAGct                                                                       >  2:3788836/1‑150 (MQ=255)
                                                                                     cATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTa                                                                  <  1:3735111/150‑1 (MQ=255)
                                                                                        tCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTa                                                               <  2:4791716/150‑1 (MQ=255)
                                                                                           gTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGGGCTGGTTGGCAATGGGTTAATCTCGCCTGGCTGGCTGGCGGAGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACg                                                            <  2:2223368/150‑1 (MQ=255)
                                                                                             tGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCt                                                          >  1:401168/1‑150 (MQ=255)
                                                                                              gAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTg                                                         <  1:2887783/150‑1 (MQ=255)
                                                                                                   gATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGtt                                                    >  2:1550049/1‑150 (MQ=255)
                                                                                                   gATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGtt                                                    >  1:4252132/1‑150 (MQ=255)
                                                                                                   gATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGtt                                                    <  1:485795/150‑1 (MQ=255)
                                                                                                      tATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATg                                                 <  2:2375001/150‑1 (MQ=255)
                                                                                                       aTGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATgc                                                >  2:1801639/1‑150 (MQ=255)
                                                                                                       aTGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATgc                                                >  1:3909551/1‑150 (MQ=255)
                                                                                                         gCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATgcgc                                              >  1:2600787/1‑150 (MQ=255)
                                                                                                          cAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCa                                             >  2:1419229/1‑150 (MQ=255)
                                                                                                           aaTATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCaa                                            >  2:3726641/1‑150 (MQ=255)
                                                                                                            atatCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAAt                                           <  1:3552082/150‑1 (MQ=255)
                                                                                                                 cGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTgg                                      >  2:4959703/1‑150 (MQ=255)
                                                                                                                 cGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTgg                                      <  1:2038613/150‑1 (MQ=255)
                                                                                                                  gATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTggg                                     >  1:5165442/1‑150 (MQ=255)
                                                                                                                  gATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCTATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTggg                                     >  1:2156451/1‑150 (MQ=255)
                                                                                                                   aTCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGc                                    <  2:1559625/150‑1 (MQ=255)
                                                                                                                    tCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCt                                   <  2:2259958/150‑1 (MQ=255)
                                                                                                                    tCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCt                                   <  1:4971747/150‑1 (MQ=255)
                                                                                                                     cTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTg                                  >  1:2330336/1‑150 (MQ=255)
                                                                                                                      tACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTgg                                 >  2:2069826/1‑150 (MQ=255)
                                                                                                                      tACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTgg                                 >  2:383817/1‑150 (MQ=255)
                                                                                                                       aCAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgt                                >  1:2777974/1‑150 (MQ=255)
                                                                                                                        cAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgtt                               >  2:4317766/1‑150 (MQ=255)
                                                                                                                         aGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgttg                              >  1:1911536/1‑150 (MQ=255)
                                                                                                                         aGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgttg                              >  2:5031289/1‑150 (MQ=255)
                                                                                                                         aGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCGGGCATATTCCCCTCAGCTTCTTAGTACCGCTGGCGTTATGCGCAATGTTGGGCTGgttg                              <  1:556092/150‑1 (MQ=255)
                                                                                                                          gCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgttgt                             >  2:4715026/1‑150 (MQ=255)
                                                                                                                            ggTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgttgtac                           >  1:4114598/1‑148 (MQ=255)
                                                                                                                            ggTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTtc                           >  1:3033918/1‑150 (MQ=255)
                                                                                                                             gTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTtct                          <  1:726645/150‑1 (MQ=255)
                                                                                                                              tATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTtctc                         <  2:2652194/150‑1 (MQ=255)
                                                                                                                               aTTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCa                        >  1:3450993/1‑150 (MQ=255)
                                                                                                                               aTTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCa                        >  2:4440317/1‑150 (MQ=255)
                                                                                                                               aTTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCa                        >  1:4645980/1‑150 (MQ=255)
                                                                                                                                       gggcgCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAac                >  2:4778259/1‑150 (MQ=255)
                                                                                                                                       gggcgCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAac                <  1:1578593/150‑1 (MQ=255)
                                                                                                                                          cgcTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACt             <  2:217504/150‑1 (MQ=255)
                                                                                                                                              ggTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGgca         >  1:5033598/1‑150 (MQ=255)
                                                                                                                                              ggTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGgca         >  1:1019966/1‑150 (MQ=255)
                                                                                                                                                 tGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGgcagca      >  1:2077782/1‑150 (MQ=255)
                                                                                                                                                    cAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCg   >  1:2266801/1‑150 (MQ=255)
                                                                                                                                                     aaTGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGc  >  2:5453886/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATATGAACACACTACGCTTAGGTATTGATGGCATTAGTCAGGCGACACCGCTGGATACATTTAAAACCTCTGTCCGTGCCGGTCATTCGGTTGAACAGATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGC  >  CP009273/1703911‑1704209

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: