Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 511,619 A→T 14.4% Q82L (CAA→CTA)  fetA → iron exporter, ATP‑binding subunit, ABC transporter FetAB subunit; peroxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273511,6190AT14.4% 702.2 / 90.8 248Q82L (CAA→CTA) fetAiron exporter, ATP‑binding subunit, ABC transporter FetAB subunit; peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base A (122/90);  new base T (16/20);  total (138/110)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

ATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTT  >  CP009273/511476‑511766
                                                                                                                                               |                                                                                                                                                   
  ttAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCTGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCtt                                                                                                                                             <  2:1186061/150‑1 (MQ=255)
  ttAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCtt                                                                                                                                             <  1:3489606/150‑1 (MQ=255)
     aGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTATCTTACt                                                                                                                                          <  1:256158/150‑1 (MQ=255)
          aTTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGcc                                                                                                                                     >  1:2693612/1‑150 (MQ=255)
             aCCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAg                                                                                                                                  >  2:506052/1‑150 (MQ=255)
               cGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGac                                                                                                                                <  1:3100310/150‑1 (MQ=255)
                ggTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGaca                                                                                                                               <  1:3561660/150‑1 (MQ=255)
                  tCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACAcc                                                                                                                             <  2:3054417/150‑1 (MQ=255)
                   ccTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCg                                                                                                                            <  2:3008004/150‑1 (MQ=255)
                    cTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGa                                                                                                                           <  1:1925528/150‑1 (MQ=255)
                    cTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGa                                                                                                                           <  2:2203563/150‑1 (MQ=255)
                        tGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCt                                                                                                                       >  2:1822812/1‑150 (MQ=255)
                         ggTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTg                                                                                                                      >  2:3478775/1‑150 (MQ=255)
                         ggTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTg                                                                                                                      >  1:998261/1‑150 (MQ=255)
                         ggTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTg                                                                                                                      >  2:1205906/1‑150 (MQ=255)
                          gTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGt                                                                                                                     <  1:3276045/150‑1 (MQ=255)
                          gTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGt                                                                                                                     >  2:2092730/1‑150 (MQ=255)
                             gtgGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTg                                                                                                                  >  2:2187583/1‑150 (MQ=255)
                             gtgGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTg                                                                                                                  <  1:3365728/150‑1 (MQ=255)
                               ggCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGc                                                                                                                >  1:2872858/1‑150 (MQ=255)
                                 cAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGa                                                                                                              <  1:2049483/150‑1 (MQ=255)
                                  aaaaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGAt                                                                                                             <  1:1519777/150‑1 (MQ=255)
                                   aaaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATa                                                                                                            <  2:4001057/150‑1 (MQ=255)
                                    aaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATAc                                                                                                           >  2:3577444/1‑150 (MQ=255)
                                    aaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATAc                                                                                                           >  1:1428661/1‑150 (MQ=255)
                                      gTACGCTCCTAAAAATAGTTGCTTCATTGATAAGCCCAACCAGCGGAACCAAACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAAACTACCCCCAACAAGACTCTTACAGCGCCCAGACACCGCCGCTGTTTGGCGATACgg                                                                                                         >  1:3308437/1‑150 (MQ=255)
                                        aCGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGta                                                                                                       >  2:3639479/1‑150 (MQ=255)
                                         cGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGtat                                                                                                      <  2:1328785/150‑1 (MQ=255)
                                          gctgctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGtata                                                                                                     >  2:4506930/1‑150 (MQ=255)
                                            tgctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACg                                                                                                   >  2:3161561/1‑150 (MQ=255)
                                              ctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGAt                                                                                                 <  2:3839577/150‑1 (MQ=255)
                                              ctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGAt                                                                                                 <  1:2756803/150‑1 (MQ=255)
                                              ctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGAt                                                                                                 <  2:1579525/150‑1 (MQ=255)
                                                aaaaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATaa                                                                                               >  1:2994935/1‑150 (MQ=255)
                                                aaaaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATaa                                                                                               <  2:4459816/150‑1 (MQ=255)
                                                 aaaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAAt                                                                                              >  2:1375129/1‑150 (MQ=255)
                                                 aaaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAAt                                                                                              <  1:1324027/150‑1 (MQ=255)
                                                   aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                            >  1:2767518/1‑150 (MQ=255)
                                                   aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                            >  2:4420554/1‑150 (MQ=255)
                                                   aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                            <  1:729100/150‑1 (MQ=255)
                                                   aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                            <  1:1986776/150‑1 (MQ=255)
                                                    aTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTg                                                                                           <  2:1808499/150‑1 (MQ=255)
                                                    aTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTg                                                                                           >  2:4177892/1‑150 (MQ=255)
                                                     tAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGa                                                                                          >  2:4329501/1‑150 (MQ=255)
                                                     tAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGa                                                                                          >  2:397600/1‑150 (MQ=255)
                                                      aGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGAt                                                                                         >  2:3185469/1‑150 (MQ=255)
                                                      aGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGAt                                                                                         <  1:2291534/150‑1 (MQ=255)
                                                      aGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGAt                                                                                         <  2:2431140/150‑1 (MQ=255)
                                                        ttGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCt                                                                                       <  2:2587832/150‑1 (MQ=255)
                                                        ttGCTTCATTGACCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCt                                                                                       <  2:4129073/150‑1 (MQ=255)
                                                         tGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCtt                                                                                      <  1:1798202/150‑1 (MQ=255)
                                                            ttCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTcc                                                                                   >  1:4043779/1‑150 (MQ=255)
                                                            ttCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTcc                                                                                   >  2:1376545/1‑150 (MQ=255)
                                                               aTTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTg                                                                                >  2:2211654/1‑150 (MQ=255)
                                                                 tGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGc                                                                              <  2:2726568/150‑1 (MQ=255)
                                                                  gATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCa                                                                             >  1:776654/1‑150 (MQ=255)
                                                                  gATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCa                                                                             <  1:4251097/150‑1 (MQ=255)
                                                                     cAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTAGACGATAATCTGATCTTTCCCTGGCAGAt                                                                          >  1:1041582/1‑150 (MQ=255)
                                                                     cAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGACTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGAt                                                                          >  1:1041358/1‑150 (MQ=255)
                                                                      aGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATc                                                                         <  2:1975218/150‑1 (MQ=255)
                                                                           aaCCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCACCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTaa                                                                    >  1:4449687/1‑150 (MQ=255)
                                                                                  ggAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAg                                                             >  1:505157/1‑150 (MQ=255)
                                                                                  ggAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAg                                                             >  1:3769448/1‑150 (MQ=255)
                                                                                    aaCGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGcc                                                           >  1:3233385/1‑150 (MQ=255)
                                                                                     aCGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCt                                                          <  1:1956116/150‑1 (MQ=255)
                                                                                     aCGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCt                                                          <  1:13955/150‑1 (MQ=255)
                                                                                     aCGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCt                                                          <  1:1058242/150‑1 (MQ=255)
                                                                                      cGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTg                                                         <  2:2626103/150‑1 (MQ=255)
                                                                                        ttACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGAc                                                       <  2:1101573/150‑1 (MQ=255)
                                                                                         tACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGAcc                                                      >  1:4388403/1‑150 (MQ=255)
                                                                                           cTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCa                                                    <  1:2187583/150‑1 (MQ=255)
                                                                                            tGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAg                                                   >  1:2931693/1‑150 (MQ=255)
                                                                                               ttGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCa                                                <  1:2264387/150‑1 (MQ=255)
                                                                                                   aGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAtttt                                            >  2:3218098/1‑150 (MQ=255)
                                                                                                    ggTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttttt                                           <  2:731310/150‑1 (MQ=255)
                                                                                                    ggTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttttt                                           >  1:3120039/1‑150 (MQ=255)
                                                                                                      tGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttct                                         <  1:40095/150‑1 (MQ=255)
                                                                                                      tGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttct                                         >  1:1668830/1‑150 (MQ=255)
                                                                                                      tGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttct                                         >  1:3963728/1‑150 (MQ=255)
                                                                                                       gAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctc                                        >  2:2921356/1‑150 (MQ=255)
                                                                                                             gTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatttt                                  >  2:2839381/1‑150 (MQ=255)
                                                                                                             gTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatttt                                  >  2:552442/1‑150 (MQ=255)
                                                                                                             gTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatttt                                  >  2:4340110/1‑150 (MQ=255)
                                                                                                               cAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgattttct                                >  1:3468614/1‑150 (MQ=255)
                                                                                                                aGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgattttctc                               <  1:47390/150‑1 (MQ=255)
                                                                                                                   acacTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGaa                            <  1:2669081/150‑1 (MQ=255)
                                                                                                                   acacTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGaa                            >  2:2870727/1‑150 (MQ=255)
                                                                                                                     acTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACg                          >  2:3760853/1‑150 (MQ=255)
                                                                                                                     acTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACg                          >  2:1500684/1‑150 (MQ=255)
                                                                                                                       tAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCt                        >  1:2513685/1‑150 (MQ=255)
                                                                                                                         aaGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTc                      >  2:483070/1‑150 (MQ=255)
                                                                                                                          aGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCg                     <  2:497278/150‑1 (MQ=255)
                                                                                                                              aGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCtt                 >  2:863828/1‑150 (MQ=255)
                                                                                                                                aaaTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGc               >  1:581064/1‑150 (MQ=255)
                                                                                                                                 aaTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGcc              >  1:1887650/1‑150 (MQ=255)
                                                                                                                                   tCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCgg            <  2:4049572/150‑1 (MQ=255)
                                                                                                                                           caacaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCGGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtt    >  1:2064550/1‑150 (MQ=255)
                                                                                                                                           caacaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtt    >  1:2142688/1‑150 (MQ=255)
                                                                                                                                            aacaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAttt   >  2:2946859/1‑150 (MQ=255)
                                                                                                                                             acaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtttt  >  2:280925/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
ATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTT  >  CP009273/511476‑511766

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: