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,5300CA98.0% 648.7 / 3.2 203E113* (GAG→TAG) fliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (104/95);  total (105/98)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
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 211 total aligned reads displayed.

TCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACC  >  CP009273/1993385‑1993674
                                                                                                                                                 |                                                                                                                                                
tCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaat                                                                                                                                                <  1:281343/148‑1 (MQ=255)
tCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaat                                                                                                                                                >  1:250464/1‑148 (MQ=255)
  aTAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCCTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaataa                                                                                                                                              <  1:433628/148‑3 (MQ=255)
     gTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGaa                                                                                                                                           <  2:2952939/148‑1 (MQ=255)
        ggTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCAACAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATg                                                                                                                                        <  1:915600/148‑1 (MQ=255)
             ccAGACCGACCCCCCCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGAc                                                                                                                                   <  2:172219/148‑1 (MQ=255)
                gaccgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTg                                                                                                                                >  2:1330565/1‑148 (MQ=255)
                gaccgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTg                                                                                                                                >  2:2404461/1‑148 (MQ=255)
                 accgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGa                                                                                                                               >  1:262538/1‑148 (MQ=255)
                  ccgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGAt                                                                                                                              <  2:2093164/148‑1 (MQ=255)
                   cgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGAtt                                                                                                                             >  2:1168488/1‑148 (MQ=255)
                   cgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGAtt                                                                                                                             >  1:1738114/1‑148 (MQ=255)
                    gaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTa                                                                                                                            <  1:941982/148‑1 (MQ=255)
                      cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAAt                                                                                                                          <  1:2835812/148‑1 (MQ=255)
                      cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAAt                                                                                                                          >  2:3234938/1‑148 (MQ=255)
                       ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATc                                                                                                                         <  2:1291231/148‑1 (MQ=255)
                        cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATca                                                                                                                        <  2:1076723/148‑1 (MQ=255)
                          cACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcacc                                                                                                                      <  2:250464/148‑1 (MQ=255)
                            cTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcaccac                                                                                                                    >  1:2486282/1‑148 (MQ=255)
                                 ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCaa                                                                                                               <  2:3037373/148‑1 (MQ=255)
                                   gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATa                                                                                                             <  1:794011/148‑1 (MQ=255)
                                      tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                                                                                                          <  1:2171930/148‑1 (MQ=255)
                                      tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                                                                                                          >  2:3390086/1‑148 (MQ=255)
                                      tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                                                                                                          <  2:1033815/148‑1 (MQ=255)
                                        tCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGttt                                                                                                        >  1:2949115/1‑148 (MQ=255)
                                           gATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTaga                                                                                                     >  1:2023188/1‑148 (MQ=255)
                                           gATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTaga                                                                                                     <  1:1161290/148‑1 (MQ=255)
                                              catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGACATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc                                                                                                  >  1:1902804/1‑148 (MQ=255)
                                              catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc                                                                                                  <  2:2253434/148‑1 (MQ=255)
                                              catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc                                                                                                  <  1:3506484/148‑1 (MQ=255)
                                               atcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTGGAGTcc                                                                                                 >  1:2952708/1‑148 (MQ=255)
                                               atcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc                                                                                                 <  2:297402/148‑1 (MQ=255)
                                               atcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc                                                                                                 >  2:2702594/1‑148 (MQ=255)
                                                tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa                                                                                                <  1:301420/148‑1 (MQ=255)
                                                tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa                                                                                                >  1:1011042/1‑148 (MQ=255)
                                                 cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAg                                                                                               >  1:3150937/1‑148 (MQ=255)
                                                 cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAg                                                                                               <  2:3540542/148‑1 (MQ=255)
                                                   gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCGAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCg                                                                                             <  2:1149903/148‑1 (MQ=255)
                                                    cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa                                                                                            >  2:3352603/1‑148 (MQ=255)
                                                       ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGc                                                                                         <  2:1365935/148‑1 (MQ=255)
                                                        tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGcc                                                                                        <  1:1592641/148‑1 (MQ=255)
                                                        tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGAGCTAATCAGAAATGGTGACCTGATTAAACACCACATCAATACGTTTAGATTCCAGCGACGcc                                                                                        >  1:2902043/1‑148 (MQ=255)
                                                          tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAg                                                                                      >  1:177955/1‑148 (MQ=255)
                                                             tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt                                                                                   >  1:1087916/1‑148 (MQ=255)
                                                                tGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATAcc                                                                                >  1:775867/1‑148 (MQ=255)
                                                                  ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGt                                                                              >  2:610981/1‑148 (MQ=255)
                                                                   cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTc                                                                             >  1:1370885/1‑148 (MQ=255)
                                                                    ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTcc                                                                            >  1:3538832/1‑148 (MQ=255)
                                                                    ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTcc                                                                            <  1:2298847/148‑1 (MQ=255)
                                                                      cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCa                                                                          >  1:1770915/1‑148 (MQ=255)
                                                                       gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAt                                                                         <  2:1531092/148‑1 (MQ=255)
                                                                       gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAt                                                                         <  2:492550/148‑1 (MQ=255)
                                                                       gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAt                                                                         <  2:1528407/148‑1 (MQ=255)
                                                                         tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAttt                                                                       >  2:3215012/1‑148 (MQ=255)
                                                                          aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTg                                                                      >  2:435247/1‑148 (MQ=255)
                                                                           ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTgg                                                                     >  1:712657/1‑148 (MQ=255)
                                                                           ccTTTTGTCACCAGCGCCTGAAGACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTgg                                                                     <  2:1850079/148‑1 (MQ=255)
                                                                             ttttttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTc                                                                   >  1:1051272/1‑148 (MQ=255)
                                                                             ttttttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTc                                                                   >  2:124013/1‑148 (MQ=255)
                                                                              tttttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCg                                                                  >  1:628626/1‑148 (MQ=255)
                                                                                    aCCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTa                                                            >  1:703076/1‑148 (MQ=255)
                                                                                       aGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTg                                                         <  1:1330565/148‑1 (MQ=255)
                                                                                          gcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACg                                                      <  2:1632046/148‑1 (MQ=255)
                                                                                           ccTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGc                                                     <  2:177955/148‑1 (MQ=255)
                                                                                            cTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGcc                                                    <  2:734543/148‑1 (MQ=255)
                                                                                            cTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGcc                                                    <  1:1168488/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCTTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   <  1:330020/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   <  1:220260/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   <  2:1624610/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   <  2:712657/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   <  1:666999/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   >  2:2387433/1‑148 (MQ=255)
                                                                                              gAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTc                                                  >  1:1151675/1‑148 (MQ=255)
                                                                                               aaTACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCa                                                 >  2:110395/1‑148 (MQ=255)
                                                                                                aTACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCaa                                                <  1:876507/148‑1 (MQ=255)
                                                                                                 tACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAAc                                               >  2:1574722/1‑148 (MQ=255)
                                                                                                  aCCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCACTACGTTTAGAGTCCAGCTACGCCAGCATACCGTCCCATGTGGTCGGTGTTAGTGACGCCTCATCg                                              >  1:3336028/1‑148 (MQ=255)
                                                                                                  aCCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACg                                              >  2:1066201/1‑148 (MQ=255)
                                                                                                     aGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCa                                           >  1:1249077/1‑148 (MQ=255)
                                                                                                      gAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCaa                                          >  1:3278770/1‑148 (MQ=255)
                                                                                                      gAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCaa                                          >  1:2910510/1‑148 (MQ=255)
                                                                                                        aaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGa                                        <  2:1738114/148‑1 (MQ=255)
                                                                                                         aTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGAt                                       <  2:1779764/148‑1 (MQ=255)
                                                                                                          tGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATg                                      >  1:3142682/1‑148 (MQ=255)
                                                                                                           ggTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGt                                     <  2:65634/148‑1 (MQ=255)
                                                                                                           ggTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGt                                     >  2:1971936/1‑148 (MQ=255)
                                                                                                               tACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTg                                 <  2:597966/148‑1 (MQ=255)
                                                                                                                 cGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGcc                               <  2:2486282/148‑1 (MQ=255)
                                                                                                                    ggTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGc                            >  2:2946662/1‑148 (MQ=255)
                                                                                                                       tGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGt                         >  1:2393666/1‑148 (MQ=255)
                                                                                                                       tGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGt                         >  1:3344915/1‑148 (MQ=255)
                                                                                                                        gagaAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGtt                        <  2:3150937/148‑1 (MQ=255)
                                                                                                                             aTCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCa                   >  2:1222146/1‑148 (MQ=255)
                                                                                                                               cGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCaaa                 >  1:504773/1‑148 (MQ=255)
                                                                                                                                  aTTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTc              >  2:3519486/1‑148 (MQ=255)
                                                                                                                                      tttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt          <  2:139347/148‑1 (MQ=255)
                                                                                                                                           tGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCaaa     >  1:926851/1‑148 (MQ=255)
                                                                                                                                           tGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCaaa     <  2:2553906/148‑1 (MQ=255)
                                                                                                                                            gcgcTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGATGTTGGGCAAATTCCACTTCaaaa    <  1:1104357/148‑1 (MQ=255)
                                                                                                                                              gcTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAcc  <  1:2554567/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACC  >  CP009273/1993385‑1993674

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: