Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,407,969 C→A 100% N249K (AAC→AAA acrE → cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,407,9690CA98.4% 583.6 / 0.9 190N249K (AAC→AAAacrEcytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (106/81);  total (107/83)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-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 190 total aligned reads displayed.

ACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGA  >  CP009273/3407825‑3408112
                                                                                                                                                |                                                                                                                                               
aCTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGTGGCTTAAGAAATCCGTAGAGAAAGTAAATTTGCATAAGGAACACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGt                                                                                                                                              <  1:985290/148‑1 (MQ=38)
    gCGTCTTTCCATCAGCTCGTTCCTATCTACGTTGATTTGACCCAATCCAGCATCTTTTTTATGAGGCTGAAGCAATCCGTAGAGCATGGAAATTTGCATATGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCaaa                                                                                                                                          <  1:1017733/148‑1 (MQ=255)
    gCGACTGTCCAGCAGCTCGATCCTATCTACTTTGATGTGACCCAATCCAGCACCGTTTTTATGTGGCTGAAGCAATCCTTAGAGCAAGGAAATTTGAATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCaaa                                                                                                                                          <  1:1017731/148‑1 (MQ=255)
    gCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCaaa                                                                                                                                          <  2:3216419/148‑1 (MQ=255)
       aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCACGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGTAAAAAGGTCAAACCt                                                                                                                                       >  1:1708097/1‑148 (MQ=255)
       aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCt                                                                                                                                       >  2:1533138/1‑148 (MQ=255)
       aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCt                                                                                                                                       <  2:592092/148‑1 (MQ=255)
        cTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTa                                                                                                                                      >  1:3679456/1‑148 (MQ=255)
           tCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATcc                                                                                                                                   >  1:1101078/1‑148 (MQ=255)
           tCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATcc                                                                                                                                   >  1:18546/1‑148 (MQ=255)
           tCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATcc                                                                                                                                   >  1:1949611/1‑148 (MQ=255)
                   gtCGATGCTATGTACTTTGATTTGACACTTTACAGCACAGTTTTTATGAGGCTTAAGCAATCCGTATATCATGTAATTTTGGATAAGGAAATCTCCAACATCAACGTAGAGTTGGTCATGGAAAAAGGTCTAACCTATACCCTGAAAg                                                                                                                           <  2:2254236/147‑1 (MQ=255)
                   cTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAg                                                                                                                           >  1:898894/1‑148 (MQ=255)
                    tCGATCCTATCTACGTTGATGTGACCCAATCCCGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAgg                                                                                                                          <  2:3613145/148‑1 (MQ=255)
                    tCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAgg                                                                                                                          <  2:1251319/148‑1 (MQ=255)
                       aTCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTAc                                                                                                                       >  1:3237960/1‑148 (MQ=255)
                         ccTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGc                                                                                                                     >  2:3668266/1‑148 (MQ=255)
                           tATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTg                                                                                                                   >  1:3349508/1‑148 (MQ=255)
                             tCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCa                                                                                                                 >  1:3583958/1‑148 (MQ=255)
                              cTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCaa                                                                                                                >  1:3023572/1‑148 (MQ=255)
                                aCGTTGATGTGACCCAATCCAGCACCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAAtt                                                                                                              <  2:2629271/148‑1 (MQ=255)
                                   ttGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtctc                                                                                                           <  1:457452/148‑1 (MQ=255)
                                    tGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTcc                                                                                                          >  1:1832659/1‑148 (MQ=255)
                                     gATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCg                                                                                                         <  2:7969/148‑1 (MQ=255)
                                      aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGa                                                                                                        <  1:2841813/148‑1 (MQ=255)
                                      aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGa                                                                                                        <  2:131845/148‑1 (MQ=255)
                                       tgtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAt                                                                                                       >  1:2680227/1‑148 (MQ=255)
                                       tgtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAt                                                                                                       <  2:1866001/148‑1 (MQ=255)
                                       tgtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAt                                                                                                       >  1:3111412/1‑148 (MQ=255)
                                       tgtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAt                                                                                                       <  2:3755211/148‑1 (MQ=255)
                                          gACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgtg                                                                                                    >  2:941971/1‑148 (MQ=255)
                                              cAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCg                                                                                                >  2:1304102/1‑148 (MQ=255)
                                                aTCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGtt                                                                                              <  1:833564/148‑1 (MQ=255)
                                                  ccaccaACGATTTTATGAGGCTTAAGCAATCCGTAGAGCAATGAAATTTGAATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtga                                                                                            <  1:157860/148‑1 (MQ=255)
                                                      cAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaa                                                                                        >  2:1120761/1‑148 (MQ=255)
                                                        aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      <  2:1097764/148‑1 (MQ=255)
                                                        aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  1:1648806/1‑148 (MQ=255)
                                                         cGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGc                                                                                     <  2:2737141/148‑1 (MQ=255)
                                                         cGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGc                                                                                     <  2:3095692/148‑1 (MQ=255)
                                                         cGATTTTAGGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGc                                                                                     <  1:892441/148‑1 (MQ=255)
                                                          taTTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCa                                                                                    <  2:1064746/147‑1 (MQ=255)
                                                                 tGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTc                                                                             >  2:1798045/1‑148 (MQ=255)
                                                                        gAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAAcc                                                                      <  2:109458/148‑1 (MQ=255)
                                                                          aGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCt                                                                    <  2:2087830/148‑1 (MQ=255)
                                                                           gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  1:2264766/1‑148 (MQ=255)
                                                                           gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  1:804690/1‑148 (MQ=255)
                                                                           gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  1:807121/1‑148 (MQ=255)
                                                                           gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  1:2941790/1‑148 (MQ=255)
                                                                           gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  1:2945015/1‑148 (MQ=255)
                                                                            cAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTAc                                                                  >  2:1023424/1‑148 (MQ=255)
                                                                              aTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGt                                                                >  2:3150573/1‑148 (MQ=255)
                                                                               tCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTg                                                               >  2:598258/1‑148 (MQ=255)
                                                                                ccGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGc                                                              >  1:3326597/1‑148 (MQ=255)
                                                                                ccGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGc                                                              >  1:884118/1‑148 (MQ=255)
                                                                                  gTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTg                                                            >  1:661609/1‑148 (MQ=255)
                                                                                   tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                           <  1:181329/148‑1 (MQ=255)
                                                                                   tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                           <  1:527845/148‑1 (MQ=255)
                                                                                    agagCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTc                                                          <  2:395003/148‑1 (MQ=255)
                                                                                     gagCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCt                                                         >  2:2671439/1‑148 (MQ=255)
                                                                                      agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        >  1:11162/1‑148 (MQ=255)
                                                                                      agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        >  2:990550/1‑148 (MQ=255)
                                                                                        cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTcc                                                      <  2:3046407/148‑1 (MQ=255)
                                                                                        cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTcc                                                      <  2:2837921/148‑1 (MQ=255)
                                                                                        cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTcc                                                      <  2:1139993/148‑1 (MQ=255)
                                                                                           ggAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCAGAACCCTACGTGCTGTCTTCCCTa                                                   >  2:2632752/1‑148 (MQ=255)
                                                                                            gAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTaa                                                  <  1:1562497/148‑1 (MQ=255)
                                                                                             aaaTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAc                                                 >  1:3705692/1‑148 (MQ=255)
                                                                                               aTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAccc                                               >  1:2400045/1‑148 (MQ=255)
                                                                                               aTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAccc                                               >  1:3786132/1‑148 (MQ=255)
                                                                                                tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCg                                              >  1:3114333/1‑148 (MQ=255)
                                                                                                 ttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGc                                             >  1:1237262/1‑148 (MQ=255)
                                                                                                 ttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGc                                             >  1:1752706/1‑148 (MQ=255)
                                                                                                      tAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                        <  1:2458076/148‑1 (MQ=255)
                                                                                                      tAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                        <  2:3583958/148‑1 (MQ=255)
                                                                                                      tAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                        <  2:393414/148‑1 (MQ=255)
                                                                                                       aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       <  2:3224130/148‑1 (MQ=255)
                                                                                                       aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       >  1:47099/1‑148 (MQ=255)
                                                                                                        aGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATAc                                      >  1:1303165/1‑148 (MQ=255)
                                                                                                             aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCttt                                 <  2:2246927/148‑1 (MQ=255)
                                                                                                               cGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTg                               >  1:1445543/1‑148 (MQ=255)
                                                                                                                     caGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTATCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGt                         <  1:554600/148‑1 (MQ=255)
                                                                                                                       gCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTACCCCGCAACATACGCTTTTGCCGGGGAt                       >  1:1061069/1‑148 (MQ=255)
                                                                                                                        cAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATg                      <  1:1280545/148‑1 (MQ=255)
                                                                                                                        cAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACAGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGACGGGTATg                      >  2:449991/1‑148 (MQ=255)
                                                                                                                         aaCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGt                     >  2:3623139/1‑148 (MQ=255)
                                                                                                                           cGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGttt                   >  1:460781/1‑148 (MQ=255)
                                                                                                                           cGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGttt                   >  1:2546688/1‑148 (MQ=255)
                                                                                                                             tAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgt                 >  2:350115/1‑148 (MQ=255)
                                                                                                                             tAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgt                 >  2:3132718/1‑148 (MQ=255)
                                                                                                                             tAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgt                 >  2:2799830/1‑148 (MQ=255)
                                                                                                                                agTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcg              >  2:1307988/1‑148 (MQ=255)
                                                                                                                                agTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcg              >  1:1999089/1‑148 (MQ=255)
                                                                                                                                  ttGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcgtg            >  1:2397725/1‑148 (MQ=255)
                                                                                                                                  ttGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcgtg            >  1:514933/1‑148 (MQ=255)
                                                                                                                                     gTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCAc         >  2:2443869/1‑148 (MQ=255)
                                                                                                                                     gTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCAc         >  2:708939/1‑148 (MQ=255)
                                                                                                                                        aTGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGa      <  2:3462525/148‑1 (MQ=255)
                                                                                                                                            aaaaAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtga  <  2:3661803/148‑1 (MQ=255)
                                                                                                                                            aaaaAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtga  <  1:2808049/148‑1 (MQ=255)
                                                                                                                                            aaaaAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtga  <  2:807121/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGA  >  CP009273/3407825‑3408112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: