Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,993,530 C→A 100% E113* (GAG→TAG)  fliY ← cystine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,993,5300CA99.2% 779.4 / ‑2.0 242E113* (GAG→TAG) fliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (130/110);  total (132/110)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 254 total aligned reads displayed.

CTTCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACCGGTT  >  CP009273/1993383‑1993678
                                                                                                                                                   |                                                                                                                                                    
 ttCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaatc                                                                                                                                                   >  2:4316180/1‑150 (MQ=255)
  tCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaatca                                                                                                                                                  <  1:4437691/150‑1 (MQ=255)
  tCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaatca                                                                                                                                                  <  2:731596/150‑1 (MQ=255)
   cATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAg                                                                                                                                                 <  2:2928697/150‑1 (MQ=255)
    aTAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGa                                                                                                                                                >  2:470577/1‑150 (MQ=255)
      aGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGaaa                                                                                                                                              >  1:4129533/1‑150 (MQ=255)
       gTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACTAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAAt                                                                                                                                             <  2:517928/150‑1 (MQ=255)
        ttGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATg                                                                                                                                            <  1:4567575/150‑1 (MQ=255)
          ggTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGt                                                                                                                                          <  1:904812/150‑1 (MQ=255)
           gTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTg                                                                                                                                         <  1:3930157/150‑1 (MQ=255)
           gTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTg                                                                                                                                         <  1:2476597/150‑1 (MQ=255)
             gCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGAc                                                                                                                                       <  2:1459478/150‑1 (MQ=255)
                   accgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGAtt                                                                                                                                 >  2:1258480/1‑150 (MQ=255)
                      gaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGACATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAAt                                                                                                                              >  2:1307363/1‑150 (MQ=255)
                      gaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAAt                                                                                                                              <  2:715657/150‑1 (MQ=255)
                       accCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATc                                                                                                                             <  2:1939325/150‑1 (MQ=255)
                          cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcacc                                                                                                                          >  1:2044419/1‑150 (MQ=255)
                          cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcacc                                                                                                                          <  2:4424777/150‑1 (MQ=255)
                           ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcacca                                                                                                                         <  1:4513384/150‑1 (MQ=255)
                             aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGGTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACcaca                                                                                                                       >  1:1396682/1‑150 (MQ=255)
                               ttttttGCCTGTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATc                                                                                                                     <  2:2585047/150‑1 (MQ=255)
                                  tttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAAt                                                                                                                  <  1:2704096/150‑1 (MQ=255)
                                  tttGCCTTTCAGATCATCGGATGTTTTAATGGTGCATTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAAt                                                                                                                  <  1:1295580/150‑1 (MQ=255)
                                   ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATa                                                                                                                 >  2:5236953/1‑150 (MQ=255)
                                    tGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATAc                                                                                                                <  2:170920/150‑1 (MQ=255)
                                     gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACg                                                                                                               <  2:4476638/150‑1 (MQ=255)
                                     gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACg                                                                                                               <  1:2175947/150‑1 (MQ=255)
                                       cTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGtt                                                                                                             <  1:2902504/150‑1 (MQ=255)
                                        tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGttt                                                                                                            <  1:3690128/150‑1 (MQ=255)
                                           cAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTaga                                                                                                         >  1:3940902/1‑150 (MQ=255)
                                            aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag                                                                                                        >  1:2630917/1‑150 (MQ=255)
                                            aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag                                                                                                        >  2:3069381/1‑150 (MQ=255)
                                                  tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGc                                                                                                  <  1:4966575/150‑1 (MQ=255)
                                                   cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCg                                                                                                 <  2:1540531/150‑1 (MQ=255)
                                                   cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCg                                                                                                 >  2:4251144/1‑150 (MQ=255)
                                                     gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGATTACCTGAAAGGGTGTACGGGGTTGAGAAATCGTTTTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGAc                                                                                               <  1:4317514/150‑1 (MQ=255)
                                                        gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGcc                                                                                            >  1:4894516/1‑150 (MQ=255)
                                                        gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGcc                                                                                            >  2:4172155/1‑150 (MQ=255)
                                                         ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCa                                                                                           <  2:4783949/150‑1 (MQ=255)
                                                         ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCa                                                                                           <  2:1963354/150‑1 (MQ=255)
                                                         ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCa                                                                                           <  1:974587/150‑1 (MQ=255)
                                                         ttttAAGGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCa                                                                                           <  1:609192/150‑1 (MQ=255)
                                                             aaTGGTGCCTTCTTTACATTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGATAATCAGAACTGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt                                                                                       <  1:188128/150‑1 (MQ=255)
                                                             aaTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt                                                                                       >  2:3219971/1‑150 (MQ=255)
                                                             aaTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt                                                                                       >  2:1352288/1‑150 (MQ=255)
                                                             aaTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt                                                                                       >  1:869221/1‑150 (MQ=255)
                                                             aaTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt                                                                                       <  1:5127502/150‑1 (MQ=255)
                                                                 gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCg                                                                                   >  1:1391388/1‑150 (MQ=255)
                                                                  tGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCTCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGt                                                                                  >  1:3426988/1‑150 (MQ=255)
                                                                    ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTcc                                                                                >  2:2836288/1‑150 (MQ=255)
                                                                     cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTccc                                                                               <  2:4828737/150‑1 (MQ=255)
                                                                      ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCa                                                                              >  2:646036/1‑150 (MQ=255)
                                                                        cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAtt                                                                            >  2:2787189/1‑150 (MQ=255)
                                                                        cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAtt                                                                            >  2:2962360/1‑150 (MQ=255)
                                                                         gtgaCCTTTTTTCACCCGCGCCTGAAGACCAGACATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAttt                                                                           <  2:3815541/147‑1 (MQ=255)
                                                                         gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAttt                                                                           >  1:1219760/1‑150 (MQ=255)
                                                                            aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGt                                                                        >  2:4425944/1‑150 (MQ=255)
                                                                            aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGt                                                                        >  2:213781/1‑150 (MQ=255)
                                                                             ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTc                                                                       <  1:4367820/150‑1 (MQ=255)
                                                                                   ttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGtttt                                                                 <  1:3450705/150‑1 (MQ=255)
                                                                                   ttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGtttt                                                                 >  2:333047/1‑150 (MQ=255)
                                                                                     ccccAGCGCCTGTATACCATAAATGTTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAATTGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAg                                                               <  2:2368273/148‑1 (MQ=255)
                                                                                      aCCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGt                                                              <  1:172334/150‑1 (MQ=255)
                                                                                       ccAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTg                                                             <  2:3722515/150‑1 (MQ=255)
                                                                                       ccAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTg                                                             <  2:4624034/150‑1 (MQ=255)
                                                                                         aGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGAc                                                           <  1:3394877/150‑1 (MQ=255)
                                                                                         aGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGAc                                                           <  2:11822/150‑1 (MQ=255)
                                                                                          gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACg                                                          <  1:3056134/150‑1 (MQ=255)
                                                                                          gcgcCTGAAGACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACg                                                          <  2:5048398/150‑1 (MQ=255)
                                                                                               tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCa                                                     >  1:4346918/1‑150 (MQ=255)
                                                                                                gAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCaa                                                    >  2:3343167/1‑150 (MQ=255)
                                                                                                  aTACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACg                                                  >  1:568866/1‑150 (MQ=255)
                                                                                                     ccAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCa                                               >  1:1007917/1‑150 (MQ=255)
                                                                                                     ccAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCa                                               >  2:1556431/1‑150 (MQ=255)
                                                                                                         aaaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGAt                                           <  2:264428/150‑1 (MQ=255)
                                                                                                         aaaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGAt                                           >  1:2907957/1‑150 (MQ=255)
                                                                                                            tGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGtt                                        >  1:1448161/1‑150 (MQ=255)
                                                                                                              gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGtttt                                      <  2:1781341/150‑1 (MQ=255)
                                                                                                               tgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTg                                     >  1:2842339/1‑150 (MQ=255)
                                                                                                                gtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACTTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGc                                    <  1:2264550/150‑1 (MQ=255)
                                                                                                                  aCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCa                                  >  1:1138865/1‑150 (MQ=255)
                                                                                                                       gTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGt                             >  2:204080/1‑150 (MQ=255)
                                                                                                                        ttGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGtt                            <  1:667804/150‑1 (MQ=255)
                                                                                                                          gagaAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTgg                          >  2:4082091/1‑150 (MQ=255)
                                                                                                                          gagaAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTgg                          >  2:1420097/1‑150 (MQ=255)
                                                                                                                           agaAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTggg                         >  2:572953/1‑150 (MQ=255)
                                                                                                                             aaaTCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCa                       >  1:1471727/1‑150 (MQ=255)
                                                                                                                              aaTCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCaa                      >  1:2723592/1‑150 (MQ=255)
                                                                                                                                tCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAAt                    <  2:869221/150‑1 (MQ=255)
                                                                                                                                 cGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAAtt                   <  2:4129520/150‑1 (MQ=255)
                                                                                                                                  gTATTTTTTCTTGCGCTAATCAGAAATGTTGACCTGATTAATCACCAAATCAATACGTTTAGAGTCCATCGACGCCAGCATACCTTACCATTTGGTCGGTTTTAGTGACGCCGCAACGCCAAGATGTTTTGCCATCTGTTGGGCAAATTc                  >  2:3946941/1‑150 (MQ=255)
                                                                                                                                      tttttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt              >  1:4259200/1‑150 (MQ=255)
                                                                                                                                       ttttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCACATTCCACtt             >  1:734004/1‑150 (MQ=255)
                                                                                                                                        tttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTc            <  1:2024314/150‑1 (MQ=255)
                                                                                                                                          tCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCaa          >  1:4662684/1‑150 (MQ=255)
                                                                                                                                           cTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCaaa         <  1:470577/150‑1 (MQ=255)
                                                                                                                                             tGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAc       <  1:375670/150‑1 (MQ=255)
                                                                                                                                             tGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAc       <  2:3633778/150‑1 (MQ=255)
                                                                                                                                                gcTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACCgg    <  1:333047/150‑1 (MQ=255)
                                                                                                                                                  tcatcaGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACCGGtt  <  2:4893738/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACCGGTT  >  CP009273/1993383‑1993678

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: ‑