Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 511,619 A→T 26.6% 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,6190AT26.6% 536.4 / 233.3 285Q82L (CAA→CTA) fetAiron exporter, ATP‑binding subunit, ABC transporter FetAB subunit; peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base A (109/100);  new base T (36/40);  total (145/140)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 297 total aligned reads displayed.

TGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTT  >  CP009273/511470‑511766
                                                                                                                                                     |                                                                                                                                                   
tGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACTCTGCTAAAAATAGTTGCTTCATTGATCATCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCcaaca                                                                                                                                                     >  2:2107109/1‑150 (MQ=255)
tGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCcaaca                                                                                                                                                     >  1:2818087/1‑150 (MQ=255)
 ggCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCcaacaa                                                                                                                                                    >  1:1654832/1‑150 (MQ=255)
   cGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGt                                                                                                                                                  <  1:2437319/150‑1 (MQ=255)
    gAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGtc                                                                                                                                                 <  2:810600/150‑1 (MQ=255)
    gAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGtc                                                                                                                                                 <  1:1045138/150‑1 (MQ=255)
     aaTTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGtct                                                                                                                                                >  1:263158/1‑150 (MQ=255)
      aTTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGtctc                                                                                                                                               <  1:4240223/150‑1 (MQ=255)
       tttAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTctct                                                                                                                                              <  1:4170951/150‑1 (MQ=255)
        ttAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCtt                                                                                                                                             <  2:2461846/150‑1 (MQ=255)
            gTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTg                                                                                                                                         >  2:1359809/1‑150 (MQ=255)
                  tACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCa                                                                                                                                   >  2:3433338/1‑150 (MQ=255)
                  tACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCa                                                                                                                                   <  1:3238635/150‑1 (MQ=255)
                    ccGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGa                                                                                                                                 >  2:3308546/1‑150 (MQ=255)
                      ggTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGaca                                                                                                                               >  2:483893/1‑150 (MQ=255)
                      ggTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGaca                                                                                                                               >  1:2013631/1‑150 (MQ=255)
                         ccTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCg                                                                                                                            <  1:2980711/150‑1 (MQ=255)
                           ttCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGAc                                                                                                                          <  1:1119414/150‑1 (MQ=255)
                           ttCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGAc                                                                                                                          <  2:2840924/150‑1 (MQ=255)
                            tCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACg                                                                                                                         >  2:1588591/1‑150 (MQ=255)
                              tGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCt                                                                                                                       <  2:616899/150‑1 (MQ=255)
                                gTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGt                                                                                                                     <  1:596140/150‑1 (MQ=255)
                                gTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGt                                                                                                                     >  2:122684/1‑150 (MQ=255)
                                 ttGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGtt                                                                                                                    <  2:1061053/150‑1 (MQ=255)
                                   gtgGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTg                                                                                                                  <  1:4398548/150‑1 (MQ=255)
                                     ggCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGc                                                                                                                >  1:792658/1‑150 (MQ=255)
                                     ggCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGc                                                                                                                <  1:510875/150‑1 (MQ=255)
                                      gCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACGGCGCCCAGACACCGACGCTGTTTGGCg                                                                                                               <  2:2600364/150‑1 (MQ=255)
                                        aaaaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGAt                                                                                                             <  1:2082190/150‑1 (MQ=255)
                                        aaaaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGAt                                                                                                             >  1:1155937/1‑150 (MQ=255)
                                           aGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACg                                                                                                          >  1:4847874/1‑150 (MQ=255)
                                           aGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACg                                                                                                          <  1:1082176/150‑1 (MQ=255)
                                             tACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGt                                                                                                        >  1:3813034/1‑150 (MQ=255)
                                              aCGCTGATAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCTGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGta                                                                                                       <  2:4435200/150‑1 (MQ=255)
                                               cGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGtat                                                                                                      <  2:3164699/150‑1 (MQ=255)
                                                  tgctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACg                                                                                                   <  1:1734930/150‑1 (MQ=255)
                                                      aaaaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATaa                                                                                               >  1:681186/1‑150 (MQ=255)
                                                      aaaaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATaa                                                                                               >  1:4744566/1‑150 (MQ=255)
                                                         aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                            >  1:3559272/1‑150 (MQ=255)
                                                         aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                            >  1:4443012/1‑150 (MQ=255)
                                                          aTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTg                                                                                           <  1:1079155/150‑1 (MQ=255)
                                                          aTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTg                                                                                           >  2:2751538/1‑150 (MQ=255)
                                                               tGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCtt                                                                                      >  1:1990520/1‑150 (MQ=255)
                                                                  ttCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTcc                                                                                   >  1:5029036/1‑150 (MQ=255)
                                                                  ttCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTcc                                                                                   <  1:564235/150‑1 (MQ=255)
                                                                     aTTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTg                                                                                <  2:197066/150‑1 (MQ=255)
                                                                     aTTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTg                                                                                <  1:1323420/150‑1 (MQ=255)
                                                                       tGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGc                                                                              >  2:3225093/1‑150 (MQ=255)
                                                                        gATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCa                                                                             >  2:4678940/1‑150 (MQ=255)
                                                                        gATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCa                                                                             >  1:2879088/1‑150 (MQ=255)
                                                                          tCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGa                                                                           >  1:1492893/1‑150 (MQ=255)
                                                                          tCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGa                                                                           >  2:1180221/1‑150 (MQ=255)
                                                                          tCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGa                                                                           >  1:2927684/1‑150 (MQ=255)
                                                                          tCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGa                                                                           >  2:1537919/1‑150 (MQ=255)
                                                                           cAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGAt                                                                          <  2:2362976/150‑1 (MQ=255)
                                                                           cAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGAt                                                                          <  2:2926811/150‑1 (MQ=255)
                                                                           cAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGAt                                                                          <  2:1453430/150‑1 (MQ=255)
                                                                            aGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATc                                                                         >  2:1790601/1‑150 (MQ=255)
                                                                            aGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATc                                                                         <  2:4262085/150‑1 (MQ=255)
                                                                                   ccAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAcc                                                                  >  1:4065257/1‑150 (MQ=255)
                                                                                   ccAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAcc                                                                  <  1:4355367/150‑1 (MQ=255)
                                                                                   ccAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAcc                                                                  <  2:4475973/150‑1 (MQ=255)
                                                                                     aGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCgg                                                                >  2:406238/1‑150 (MQ=255)
                                                                                      gCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGc                                                               <  2:549926/150‑1 (MQ=255)
                                                                                       cGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCa                                                              <  2:377722/150‑1 (MQ=255)
                                                                                        ggAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAg                                                             >  2:268425/1‑150 (MQ=255)
                                                                                            cGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTg                                                         >  1:2295429/1‑150 (MQ=255)
                                                                                            cGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTg                                                         <  1:71323/150‑1 (MQ=255)
                                                                                               tACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGAcc                                                      <  2:786488/150‑1 (MQ=255)
                                                                                                aCTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGAccc                                                     <  1:3114398/150‑1 (MQ=255)
                                                                                                  tGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAg                                                   <  1:4541180/150‑1 (MQ=255)
                                                                                                     ttGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCa                                                >  2:4308935/1‑150 (MQ=255)
                                                                                                         aGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAtttt                                            <  1:342276/150‑1 (MQ=255)
                                                                                                                 aTGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatt                                    >  1:537606/1‑150 (MQ=255)
                                                                                                                   gTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatttt                                  <  2:1876973/150‑1 (MQ=255)
                                                                                                                    tCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgattttc                                 <  2:3289108/150‑1 (MQ=255)
                                                                                                                       gCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgattttctcg                              <  1:4938344/150‑1 (MQ=255)
                                                                                                                         acacTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGaa                            >  1:667990/1‑150 (MQ=255)
                                                                                                                          cacTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAAc                           <  2:4517990/150‑1 (MQ=255)
                                                                                                                            cTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGc                         >  1:1845032/1‑150 (MQ=255)
                                                                                                                             tAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCt                        <  1:4463912/150‑1 (MQ=255)
                                                                                                                              aaaGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCtt                       <  2:3941495/150‑1 (MQ=255)
                                                                                                                               aaGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTc                      <  1:1815251/150‑1 (MQ=255)
                                                                                                                                aGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCg                     <  2:3693167/150‑1 (MQ=255)
                                                                                                                                  ccAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGcc                   >  1:2212343/1‑150 (MQ=255)
                                                                                                                                   cAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCt                  >  2:767204/1‑150 (MQ=255)
                                                                                                                                    aNAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGGTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCtt                 >  1:19649/1‑150 (MQ=255)
                                                                                                                                    aGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCtt                 >  1:3630490/1‑150 (MQ=255)
                                                                                                                                    aGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCtt                 >  1:4646342/1‑150 (MQ=255)
                                                                                                                                      aaaTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCGGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGc               <  2:1145861/150‑1 (MQ=255)
                                                                                                                                       aaTCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGcc              <  2:5058551/150‑1 (MQ=255)
                                                                                                                                       aaTCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGcc              <  1:3918535/150‑1 (MQ=255)
                                                                                                                                       aaTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGcc              <  2:1020948/150‑1 (MQ=255)
                                                                                                                                         tCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCgg            >  2:3027141/1‑150 (MQ=255)
                                                                                                                                          cTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGa           >  2:1752930/1‑150 (MQ=255)
                                                                                                                                            aCCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACa         <  2:4177479/150‑1 (MQ=255)
                                                                                                                                             ccGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAg        <  2:3912093/150‑1 (MQ=255)
                                                                                                                                             ccGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAg        >  1:2340775/1‑150 (MQ=255)
                                                                                                                                                 caacaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtt    >  1:4211407/1‑150 (MQ=255)
                                                                                                                                                 caacaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtt    >  2:2426630/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTT  >  CP009273/511470‑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: