Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 511,619 A→T 33.7% 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,6190AT33.7% 289.2 / 244.7 225Q82L (CAA→CTA) fetAiron exporter, ATP‑binding subunit, ABC transporter FetAB subunit; peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base A (60/89);  new base T (46/30);  total (106/119)
Fisher's exact test for biased strand distribution p-value = 4.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 237 total aligned reads displayed.

TGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTTG  >  CP009273/511470‑511767
                                                                                                                                                     |                                                                                                                                                    
   cGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGt                                                                                                                                                   >  1:3210225/1‑150 (MQ=255)
   cGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGt                                                                                                                                                   >  1:2032828/1‑150 (MQ=255)
   cGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGt                                                                                                                                                   >  1:1637280/1‑150 (MQ=255)
    gAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGtc                                                                                                                                                  >  2:3862428/1‑150 (MQ=255)
     aaTTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGtct                                                                                                                                                 >  1:4104668/1‑150 (MQ=255)
     aaTTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACAAAAGCCAGAAATCTACCGCCAACTAGtct                                                                                                                                                 >  1:1386346/1‑150 (MQ=255)
         tAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTa                                                                                                                                             >  2:3909929/1‑150 (MQ=255)
          aaGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTAc                                                                                                                                            <  1:1018619/150‑1 (MQ=255)
           aGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACt                                                                                                                                           <  2:1520519/150‑1 (MQ=255)
               aaTTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTgcgc                                                                                                                                       >  1:2423028/1‑150 (MQ=255)
                 ttACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGccc                                                                                                                                     >  1:2239631/1‑150 (MQ=255)
                  tACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCa                                                                                                                                    >  1:3572341/1‑150 (MQ=255)
                     cGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGac                                                                                                                                 <  1:1771770/150‑1 (MQ=255)
                       gTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGacac                                                                                                                               <  2:320128/150‑1 (MQ=255)
                        tCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACAcc                                                                                                                              >  2:735775/1‑150 (MQ=255)
                        tCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACAcc                                                                                                                              >  1:505702/1‑150 (MQ=255)
                         ccTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCg                                                                                                                             <  2:1073101/150‑1 (MQ=255)
                         ccTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCg                                                                                                                             <  2:2597872/150‑1 (MQ=255)
                         ccTTCTGGTTGTGGAAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCg                                                                                                                             <  2:1834872/150‑1 (MQ=255)
                          cTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGa                                                                                                                            >  1:2329733/1‑150 (MQ=255)
                             cTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGc                                                                                                                         >  1:1956097/1‑150 (MQ=255)
                               ggTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTg                                                                                                                       >  2:4124573/1‑150 (MQ=255)
                                   gtgGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTg                                                                                                                   >  1:2404522/1‑150 (MQ=255)
                                         aaaGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATa                                                                                                             <  1:92716/150‑1 (MQ=255)
                                            gTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACgg                                                                                                          <  2:2039125/150‑1 (MQ=255)
                                             tACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGt                                                                                                         <  2:1267521/150‑1 (MQ=255)
                                             tACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGt                                                                                                         >  1:3411668/1‑150 (MQ=255)
                                               cGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGtat                                                                                                       <  2:1008064/150‑1 (MQ=255)
                                               cGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGtat                                                                                                       <  1:2450907/150‑1 (MQ=255)
                                                 ctgctAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATAc                                                                                                     >  2:3726162/1‑150 (MQ=255)
                                                     tAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATa                                                                                                 <  1:2896821/150‑1 (MQ=255)
                                                     tAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATa                                                                                                 <  2:2488929/150‑1 (MQ=255)
                                                        aaaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGCTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATc                                                                                              >  1:2533659/1‑150 (MQ=255)
                                                         aaTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCt                                                                                             >  1:167079/1‑150 (MQ=255)
                                                          aTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTg                                                                                            >  1:3613485/1‑150 (MQ=255)
                                                          aTAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTGTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTg                                                                                            >  1:267242/1‑150 (MQ=255)
                                                           tAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTCCTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGa                                                                                           >  2:3233001/1‑150 (MQ=255)
                                                            aGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGAt                                                                                          <  1:3046908/150‑1 (MQ=255)
                                                                    cATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCt                                                                                  <  2:2729892/150‑1 (MQ=255)
                                                                     aTTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTg                                                                                 <  2:1884354/150‑1 (MQ=255)
                                                                       tGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGc                                                                               >  1:765617/1‑150 (MQ=255)
                                                                       tGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGc                                                                               <  2:709307/150‑1 (MQ=255)
                                                                       tGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATCATCTGATCTTTCCCTGGc                                                                               >  1:2112384/1‑150 (MQ=255)
                                                                          tCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGa                                                                            >  2:3488399/1‑150 (MQ=255)
                                                                           cAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGAt                                                                           >  2:2297058/1‑150 (MQ=255)
                                                                           cAGCCCAAACAGCGGAACGGTTCTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACACGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATTCGATAATCTGATCTTGCCCTGGCAGAt                                                                           <  1:441649/150‑1 (MQ=255)
                                                                            aGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATc                                                                          >  2:3097104/1‑150 (MQ=255)
                                                                               ccAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGt                                                                       >  1:377911/1‑150 (MQ=255)
                                                                               ccAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGt                                                                       <  1:3826098/150‑1 (MQ=255)
                                                                               ccAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGt                                                                       <  1:3074598/150‑1 (MQ=255)
                                                                                 aaCCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTaa                                                                     >  2:3303514/1‑150 (MQ=255)
                                                                                 aaCCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTaa                                                                     >  2:3626218/1‑150 (MQ=255)
                                                                                  aCCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAc                                                                    >  2:610943/1‑150 (MQ=255)
                                                                                  aCCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAc                                                                    >  2:341474/1‑150 (MQ=255)
                                                                                  aCCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAc                                                                    >  1:650655/1‑150 (MQ=255)
                                                                                   ccAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAcc                                                                   <  2:3732244/150‑1 (MQ=255)
                                                                                        ggAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAg                                                              >  2:3615183/1‑150 (MQ=255)
                                                                                         gAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGc                                                             >  2:3868207/1‑150 (MQ=255)
                                                                                           aCGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCt                                                           <  2:3585218/150‑1 (MQ=255)
                                                                                               tACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGAcc                                                       <  2:869528/150‑1 (MQ=255)
                                                                                               tACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGAcc                                                       <  2:3839941/150‑1 (MQ=255)
                                                                                                 cTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCa                                                     <  1:1937289/150‑1 (MQ=255)
                                                                                                  tGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAg                                                    >  1:2092299/1‑150 (MQ=255)
                                                                                                   gTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAAACGGCAGCCTGACCCAGc                                                   >  1:3937769/1‑150 (MQ=255)
                                                                                                     ttGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCa                                                 <  2:3540012/150‑1 (MQ=255)
                                                                                                     ttGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCa                                                 >  2:1221007/1‑150 (MQ=255)
                                                                                                     ttGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCa                                                 >  1:4120341/1‑150 (MQ=255)
                                                                                                       gAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAtt                                               <  1:9398/150‑1 (MQ=255)
                                                                                                        aaGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttt                                              >  1:1569821/1‑150 (MQ=255)
                                                                                                        aaGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttt                                              >  2:4157008/1‑150 (MQ=255)
                                                                                                        aaGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttt                                              <  2:1119824/150‑1 (MQ=255)
                                                                                                        aaGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttt                                              <  2:2371725/150‑1 (MQ=255)
                                                                                                        aaGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttt                                              >  2:2165096/1‑150 (MQ=255)
                                                                                                         aGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAtttt                                             <  2:3857741/150‑1 (MQ=255)
                                                                                                          ggTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCTCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttttt                                            >  1:1557978/1‑150 (MQ=255)
                                                                                                          ggTGAGGATGTCAGCAAACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCAttttt                                            <  2:492436/150‑1 (MQ=255)
                                                                                                            tGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttct                                          <  2:2983347/150‑1 (MQ=255)
                                                                                                                 aTGTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatt                                     >  1:1848994/1‑150 (MQ=255)
                                                                                                                 aTGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatt                                     >  1:31475/1‑150 (MQ=255)
                                                                                                                 aTGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatt                                     >  1:593249/1‑150 (MQ=255)
                                                                                                                   gTCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgatttt                                   >  1:2773834/1‑150 (MQ=255)
                                                                                                                    tCAGCACACTAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgattttc                                  <  2:926809/150‑1 (MQ=255)
                                                                                                                      aGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATttttctcgattttctc                                <  2:1665094/150‑1 (MQ=255)
                                                                                                                          cacTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAAc                            <  2:1610807/150‑1 (MQ=255)
                                                                                                                          cacTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAAc                            <  1:1899766/150‑1 (MQ=255)
                                                                                                                           acTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACg                           >  2:1392760/1‑150 (MQ=255)
                                                                                                                           acTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACg                           <  2:2762918/150‑1 (MQ=255)
                                                                                                                           acTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACg                           <  2:660312/150‑1 (MQ=255)
                                                                                                                           acTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACg                           <  2:660591/150‑1 (MQ=255)
                                                                                                                             tAAAGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCt                         >  2:745090/1‑150 (MQ=255)
                                                                                                                               aaGCCAGAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTc                       <  1:2718348/150‑1 (MQ=255)
                                                                                                                                     gAAATCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTg                 >  1:1055607/1‑150 (MQ=255)
                                                                                                                                       aaTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGcc               <  1:423110/150‑1 (MQ=255)
                                                                                                                                       aaTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGcc               <  1:2923411/150‑1 (MQ=255)
                                                                                                                                        aTCTACCGCCAACTAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCg              <  1:160608/150‑1 (MQ=255)
                                                                                                                                        aTCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCg              <  1:2764158/150‑1 (MQ=255)
                                                                                                                                               gCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCa       >  2:2029872/1‑150 (MQ=255)
                                                                                                                                                   acaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtttt   <  2:1390034/150‑1 (MQ=255)
                                                                                                                                                   acaaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCAtttt   >  1:1138632/1‑150 (MQ=255)
                                                                                                                                                    caaGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTTg  <  2:1316919/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGATGTCAGCACACTAAAGCCAGAAATCTACCGCCAACAAGTCTCTTACTGCGCCCAGACACCGACGCTGTTTGGCGATACGGTATACGATAATCTGATCTTTCCCTGGCAGATCCGTAACCGGCAGCCTGACCCAGCCATTTTTCTCGATTTTCTCGAACGCTTCGCCTTGCCGGACAGCATTTTG  >  CP009273/511470‑511767

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: