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,5300CA97.1% 535.2 / 6.6 171E113* (GAG→TAG) fliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (76/90);  total (79/92)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 177 total aligned reads displayed.

TCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAACC  >  CP009273/1993385‑1993674
                                                                                                                                                 |                                                                                                                                                
tCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaat                                                                                                                                                <  2:1387098/148‑1 (MQ=255)
 cATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaatc                                                                                                                                               >  2:1547408/1‑148 (MQ=255)
 cATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaatc                                                                                                                                               <  1:3322388/148‑1 (MQ=255)
 cATAGTTGGTGCCCAGACCGACCCCCAATTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtaatc                                                                                                                                               <  1:2636201/148‑1 (MQ=255)
    aGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGa                                                                                                                                            >  1:531549/1‑148 (MQ=255)
    aGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGa                                                                                                                                            >  1:955082/1‑148 (MQ=255)
     gTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGaa                                                                                                                                           <  1:3487067/148‑1 (MQ=255)
     gTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGaa                                                                                                                                           >  1:1486258/1‑148 (MQ=255)
      ttGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGaaa                                                                                                                                          <  2:319881/148‑1 (MQ=255)
       tGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAAt                                                                                                                                         <  1:1028854/148‑1 (MQ=255)
        ggTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATg                                                                                                                                        <  1:1008134/148‑1 (MQ=255)
         gTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATgg                                                                                                                                       >  1:1099577/1‑148 (MQ=255)
         gTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATgg                                                                                                                                       >  1:1503995/1‑148 (MQ=255)
         gTGCCCAGACCGACCCCCACTTTATTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATgg                                                                                                                                       <  2:2564092/148‑1 (MQ=255)
           gCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTg                                                                                                                                     <  1:2111661/148‑1 (MQ=255)
            cccAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGa                                                                                                                                    <  2:1606126/148‑1 (MQ=255)
             ccAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGAc                                                                                                                                   <  1:202889/148‑1 (MQ=255)
              cAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCATGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATTGTGAcc                                                                                                                                  >  2:2524702/1‑148 (MQ=255)
                gaccgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                                                                                                                <  2:397427/148‑1 (MQ=255)
                gaccgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTg                                                                                                                                <  2:1587822/148‑1 (MQ=255)
                 accgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGa                                                                                                                               <  1:1021435/148‑1 (MQ=255)
                   cgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGAtt                                                                                                                             <  1:252401/148‑1 (MQ=255)
                    gaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTa                                                                                                                            <  1:292472/148‑1 (MQ=255)
                       ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATc                                                                                                                         <  1:3387467/148‑1 (MQ=255)
                        cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATca                                                                                                                        <  1:3547922/148‑1 (MQ=255)
                          cACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcacc                                                                                                                      <  1:227895/148‑1 (MQ=255)
                          cACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATaacc                                                                                                                      >  1:1836654/1‑148 (MQ=255)
                           aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcacca                                                                                                                     <  1:2087921/148‑1 (MQ=255)
                            cTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcaccac                                                                                                                    >  2:2079442/1‑148 (MQ=255)
                            cTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATcaccac                                                                                                                    <  2:1663685/148‑1 (MQ=255)
                             ttttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACcaca                                                                                                                   >  2:1911698/1‑148 (MQ=255)
                               ttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGATAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATc                                                                                                                 <  2:2629197/148‑1 (MQ=255)
                               ttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATc                                                                                                                 <  2:3513954/148‑1 (MQ=255)
                                 ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCaa                                                                                                               >  1:826945/1‑148 (MQ=255)
                                      tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                                                                                                          <  1:3206531/148‑1 (MQ=255)
                                       ttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGtt                                                                                                         <  1:3511092/148‑1 (MQ=255)
                                          aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTag                                                                                                      >  2:187450/1‑148 (MQ=255)
                                          aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTag                                                                                                      >  2:927401/1‑148 (MQ=255)
                                           gATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTaga                                                                                                     <  1:187633/148‑1 (MQ=255)
                                            atcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag                                                                                                    >  1:2214456/1‑148 (MQ=255)
                                            atcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag                                                                                                    <  1:385995/148‑1 (MQ=255)
                                             tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGt                                                                                                   >  2:1393912/1‑148 (MQ=255)
                                                    cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa                                                                                            <  2:1317805/148‑1 (MQ=255)
                                                    cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa                                                                                            >  1:2771694/1‑148 (MQ=255)
                                                    cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa                                                                                            >  1:2673494/1‑148 (MQ=255)
                                                      gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACg                                                                                          <  2:2344095/148‑1 (MQ=255)
                                                      gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACg                                                                                          >  1:57250/1‑148 (MQ=255)
                                                      gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACg                                                                                          >  1:2798227/1‑148 (MQ=255)
                                                      gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACg                                                                                          >  1:3383992/1‑148 (MQ=255)
                                                        tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGcc                                                                                        >  1:792859/1‑148 (MQ=255)
                                                           aaTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCaac                                                                                     <  2:3270420/148‑3 (MQ=255)
                                                            aTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCa                                                                                    >  1:3462233/1‑148 (MQ=255)
                                                                  ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGt                                                                              <  2:2742797/148‑1 (MQ=255)
                                                                   cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTc                                                                             >  1:2367766/1‑148 (MQ=255)
                                                                    ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTcc                                                                            >  1:2191207/1‑148 (MQ=255)
                                                                      cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATCCGTTTAGCGTCCAGCGACGCCAGCATACCGTCCCa                                                                          <  1:2994058/148‑1 (MQ=255)
                                                                        ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAtt                                                                        >  2:476605/1‑148 (MQ=255)
                                                                        ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAtt                                                                        >  2:2444231/1‑148 (MQ=255)
                                                                          aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTg                                                                      >  2:2711871/1‑148 (MQ=255)
                                                                          aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTg                                                                      >  1:2623231/1‑148 (MQ=255)
                                                                             ttttttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTc                                                                   >  2:278824/1‑148 (MQ=255)
                                                                              tttttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCg                                                                  >  1:1244148/1‑148 (MQ=255)
                                                                              tttttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCg                                                                  <  2:757763/148‑1 (MQ=255)
                                                                                tttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGt                                                                <  1:391318/148‑1 (MQ=255)
                                                                                  tCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGttt                                                              <  1:2011137/148‑1 (MQ=255)
                                                                                   cACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGtttt                                                             <  2:1989249/148‑1 (MQ=255)
                                                                                    aCCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTa                                                            >  2:242179/1‑148 (MQ=255)
                                                                                     ccAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAg                                                           <  2:3543234/148‑1 (MQ=255)
                                                                                        gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGa                                                        <  2:945767/148‑1 (MQ=255)
                                                                                           ccTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGc                                                     <  2:1626770/148‑1 (MQ=255)
                                                                                             tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCt                                                   <  1:926930/148‑1 (MQ=255)
                                                                                              gAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTc                                                  >  2:2956472/1‑148 (MQ=255)
                                                                                               aaTACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCa                                                 <  1:1086530/148‑1 (MQ=255)
                                                                                               aaTACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCa                                                 <  2:1377050/148‑1 (MQ=255)
                                                                                               aaTACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCa                                                 >  2:912714/1‑148 (MQ=255)
                                                                                                     aGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCa                                           <  1:1393030/148‑1 (MQ=255)
                                                                                                       aaaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAg                                         <  1:3093025/148‑1 (MQ=255)
                                                                                                        aaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGa                                        <  2:1503995/148‑1 (MQ=255)
                                                                                                        aaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGa                                        <  2:1099577/148‑1 (MQ=255)
                                                                                                          tGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATCTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATg                                      <  1:3314994/148‑1 (MQ=255)
                                                                                                           ggTGTACGTGGTTGTGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTTATTAATCACCACATCAATCCGTTTATAGTCCAGAGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGt                                     <  1:634233/148‑1 (MQ=255)
                                                                                                                    ggTTGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGc                            >  2:3314116/1‑148 (MQ=255)
                                                                                                                      ttGAGAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTg                          >  2:434271/1‑148 (MQ=255)
                                                                                                                       tGATAAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGt                         <  1:1549075/148‑1 (MQ=255)
                                                                                                                        gagaAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGtt                        <  2:2388843/148‑1 (MQ=255)
                                                                                                                         agaAATCGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTg                       <  1:2862296/148‑1 (MQ=255)
                                                                                                                               cGTATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCaaa                 >  1:1735308/1‑148 (MQ=255)
                                                                                                                                 tATTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAAtt               >  2:620015/1‑148 (MQ=255)
                                                                                                                                  aTTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTc              <  1:3517251/148‑1 (MQ=255)
                                                                                                                                  aTTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTc              <  2:3396795/148‑1 (MQ=255)
                                                                                                                                  aTTTTTTCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTc              >  2:3112838/1‑148 (MQ=255)
                                                                                                                                    tttttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTCTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCa            >  2:1002069/1‑148 (MQ=255)
                                                                                                                                      tttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt          <  1:1102772/148‑1 (MQ=255)
                                                                                                                                      tttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt          >  1:3435793/1‑148 (MQ=255)
                                                                                                                                      tttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt          <  1:27133/148‑1 (MQ=255)
                                                                                                                                       ttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACtt         >  2:2812992/1‑148 (MQ=255)
                                                                                                                                       ttCTTGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACtt         >  2:1037622/1‑148 (MQ=255)
                                                                                                                                           tGCGCTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCACATTCCACTTCaaa     >  1:1733659/1‑148 (MQ=255)
                                                                                                                                              gcTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAcc  >  1:1653373/1‑148 (MQ=255)
                                                                                                                                              gcTAATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAcc  >  1:1653432/1‑148 (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: