Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,351,022 G→T 100% R116S (CGT→AGT)  sapA ← antimicrobial peptide transport ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,351,0220GT98.3% 543.8 / 1.1 178R116S (CGT→AGT) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (103/72);  total (104/74)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-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 178 total aligned reads displayed.

TCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTC  >  CP009273/1350877‑1351165
                                                                                                                                                 |                                                                                                                                               
  cAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTg                                                                                                                                             <  1:925095/148‑1 (MQ=255)
  cAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTg                                                                                                                                             >  1:2098103/1‑148 (MQ=255)
    gTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTggg                                                                                                                                           <  1:800526/148‑1 (MQ=255)
            actcTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAcc                                                                                                                                   >  2:1852528/4‑148 (MQ=255)
            actcTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAcc                                                                                                                                   >  2:309247/4‑148 (MQ=255)
            actcTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAcc                                                                                                                                   >  2:17070/4‑148 (MQ=255)
                      cGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGtttt                                                                                                                         >  2:2780322/1‑148 (MQ=255)
                       gTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGttttt                                                                                                                        <  1:2984547/148‑1 (MQ=255)
                        ttATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGtttttt                                                                                                                       >  2:3178902/1‑148 (MQ=255)
                            cGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGaaa                                                                                                                   <  1:691061/148‑1 (MQ=255)
                                aaaTTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGa                                                                                                               <  2:160238/148‑1 (MQ=255)
                                  aTTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAAc                                                                                                             >  2:986577/1‑148 (MQ=255)
                                  aTTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAAc                                                                                                             >  2:2829451/1‑148 (MQ=255)
                                    tGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACAt                                                                                                           >  2:1975513/1‑148 (MQ=255)
                                         gCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAAAAACGCGa                                                                                                      >  1:2032107/1‑148 (MQ=255)
                                                aaGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGt                                                                                               >  2:1142707/1‑148 (MQ=255)
                                                aaGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGt                                                                                               >  2:2060593/1‑148 (MQ=255)
                                                  gTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGTAAAGGAACATCGCGACGCAGGTgg                                                                                             >  1:2928095/1‑148 (MQ=255)
                                                   tAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGa                                                                                            <  2:1892087/148‑1 (MQ=255)
                                                    aGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGaa                                                                                           <  2:1842681/148‑1 (MQ=255)
                                                     gttgAAGTTGCTGCCGTTGACGTTATGCCTGGTTTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAg                                                                                          <  1:2566141/145‑1 (MQ=255)
                                                        gAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACATAGCGACGCAGGTGGAAGCGa                                                                                       >  1:579664/1‑148 (MQ=255)
                                                         aaGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGAt                                                                                      <  1:2560057/148‑1 (MQ=255)
                                                         aaGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGAt                                                                                      <  1:3090767/148‑1 (MQ=255)
                                                          aGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATa                                                                                     <  1:2235518/148‑1 (MQ=255)
                                                            ttGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAgg                                                                                   >  2:365302/1‑148 (MQ=255)
                                                            ttGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAgg                                                                                   >  2:275141/1‑148 (MQ=255)
                                                               ctgcCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTcg                                                                                <  1:3395767/148‑1 (MQ=255)
                                                                 gcCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTcgcg                                                                              >  2:2483312/1‑148 (MQ=255)
                                                                    gTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCg                                                                           <  2:2434494/148‑1 (MQ=255)
                                                                     ttGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGt                                                                          <  2:2098103/148‑1 (MQ=255)
                                                                     ttGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGt                                                                          >  1:3050087/1‑148 (MQ=255)
                                                                      tGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGtt                                                                         >  2:1951548/1‑148 (MQ=255)
                                                                      tGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGtt                                                                         <  1:1513760/148‑1 (MQ=255)
                                                                      tGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGtt                                                                         >  2:153746/1‑148 (MQ=255)
                                                                       gtcGTTATGCCAGGGGTTGTTGCGGTCTAAAATTCGCTGGATGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCTGTAGTTTTTTTGTAACGGAACTTCGCGTCGCAGGTGGAAGCGATAGGTCGCGCCGTTg                                                                        <  1:2262883/146‑1 (MQ=255)
                                                                          gTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCg                                                                     >  2:1879602/1‑148 (MQ=255)
                                                                           ttATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGGTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGGTGTCGa                                                                    >  1:1371253/1‑148 (MQ=255)
                                                                            tATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAg                                                                   >  2:3123528/1‑148 (MQ=255)
                                                                              tGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGGTGGCGAGTa                                                                 >  2:3510029/1‑148 (MQ=255)
                                                                                 cACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACATCGCGACGCAGGTGGAAGCGAGAGGTCGCGCCGTTGTCGAGTACtt                                                              >  1:3142553/1‑148 (MQ=255)
                                                                                      gttgttGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAg                                                         <  2:112144/148‑1 (MQ=255)
                                                                                      gttgttGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAg                                                         >  1:3068217/1‑148 (MQ=255)
                                                                                       ttgttgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGc                                                        >  2:3390291/1‑148 (MQ=255)
                                                                                       ttgttgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGc                                                        >  2:3390262/1‑148 (MQ=255)
                                                                                       ttgttgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGc                                                        <  1:1725755/148‑1 (MQ=255)
                                                                                        tgttgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCt                                                       >  2:438440/1‑148 (MQ=255)
                                                                                          ttgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCttt                                                     <  1:97165/148‑1 (MQ=255)
                                                                                           tgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTc                                                    >  1:3254486/1‑148 (MQ=255)
                                                                                           tgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTc                                                    >  1:1012269/1‑148 (MQ=255)
                                                                                           tgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCCATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTc                                                    >  1:2801020/1‑148 (MQ=255)
                                                                                               gTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGTAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGGTGTGGAGTACTTCCCAGCTTTAGGCa                                                >  1:1339674/1‑148 (MQ=255)
                                                                                                 cAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAg                                              <  2:3369930/148‑1 (MQ=255)
                                                                                                 cAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAg                                              >  2:1357059/1‑148 (MQ=255)
                                                                                                   aaaaTTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGtt                                            >  2:3046039/1‑148 (MQ=255)
                                                                                                   aaaaTTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGtt                                            >  1:2600007/1‑148 (MQ=255)
                                                                                                        tCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTTTAAAGGAACATAGCGACGCAGGTGGAAGCGATATGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTGGGCAAGTTGCGGc                                       >  1:1634160/1‑148 (MQ=255)
                                                                                                        tCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGTAACGGAACATCGCGACGCAGGTGGTAGCGATAGGTCGCGCCGTTGTCGAGTACTTACCACCTTTCGGCAAGTTCCGGc                                       >  1:2196911/1‑148 (MQ=255)
                                                                                                        tCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGc                                       >  1:550387/1‑148 (MQ=255)
                                                                                                           cTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATc                                    <  2:799622/148‑1 (MQ=255)
                                                                                                             ggAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACAACGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAg                                  >  2:2881273/1‑148 (MQ=255)
                                                                                                              gAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAgg                                 >  1:284122/1‑148 (MQ=255)
                                                                                                              gAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAgg                                 <  1:834824/148‑1 (MQ=255)
                                                                                                               aaGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTTAAACGGAAAATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGGTACGGCATCAGGc                                >  1:3032755/1‑148 (MQ=255)
                                                                                                               aaGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGTAAAGGAACAACGCCACGCAGGTGGAATCGAGAGGTCGCGCCGGTGTGCATTACTTCCCCGCTTTCGGCAAGTTCCGGCCTCTGGc                                >  1:3032712/1‑148 (MQ=255)
                                                                                                                 ggTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGa                              >  2:669878/1‑148 (MQ=255)
                                                                                                                 ggTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGa                              >  2:484391/1‑148 (MQ=255)
                                                                                                                 ggTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGa                              >  2:484388/1‑148 (MQ=255)
                                                                                                                   tAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTTTAAAGGAACATAGCGACCCAGGTGGGAGCGCTAGGTCGCGCCGGTGTGGAGTACTTACCCGCTTTGGGCAAGTTACGGCATCACGCGATa                            >  1:3335995/1‑148 (MQ=255)
                                                                                                                   tAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATa                            >  2:3154534/1‑148 (MQ=255)
                                                                                                                        accacATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGttt                       >  2:2128170/1‑146 (MQ=255)
                                                                                                                         ccacATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGtata                      <  1:1735094/148‑1 (MQ=255)
                                                                                                                         ccacATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCCTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGtata                      >  1:2639448/1‑148 (MQ=255)
                                                                                                                          cacaTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAg                     >  2:1161918/1‑148 (MQ=255)
                                                                                                                          cacaTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATAGCGACGCAGGTGGAAGCGCTAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCCTCAGGCGATAGGTATAg                     >  1:2443838/1‑148 (MQ=255)
                                                                                                                            caTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACCGAACAACCCGACGCACGTGGAAGCGATAGGTCGCGCCCTTGTCGAGTAATTACCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAggg                   >  1:800110/1‑148 (MQ=255)
                                                                                                                             aTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGa                  <  2:257408/148‑1 (MQ=255)
                                                                                                                             aTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGa                  <  2:256703/148‑1 (MQ=255)
                                                                                                                             aTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGa                  >  2:2375285/1‑148 (MQ=255)
                                                                                                                             aTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGa                  <  1:2241686/148‑1 (MQ=255)
                                                                                                                              tcgtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGAt                 <  1:393944/148‑1 (MQ=255)
                                                                                                                               cgtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTTAAACGGAACATCGCGACGCAGGTGGCACCGCTAGGTCGCGCCGTTGTGGAGTACTTCCCAGCTTTAGGCAAGTTCCGGCCTCAGGCGATAGGTATATGGATc                >  2:659167/1‑148 (MQ=255)
                                                                                                                               cgtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGTAAAGGAACAACGCGACGCAGGTGGAAGCGCTAGGTCGCGCCCTTGTCGAGTAAATCCCACCATTCGGCAAGGTACGGCATCAGGCGATAGGTATATGGATc                >  1:1117599/1‑148 (MQ=255)
                                                                                                                               cgtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATATGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATc                >  2:3461572/1‑148 (MQ=255)
                                                                                                                                 tcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTTAAAAGGAACATAGCGACGCAGGTGGAACCGATAGGTAGCGCCCGTGTCGAGTAATTACCACCTTTCGGCAAGTTCCGGCATAACGCGCTAGGTATAGGGATAGa              >  1:1244054/1‑148 (MQ=255)
                                                                                                                                  cgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGAc             <  2:1434108/148‑1 (MQ=255)
                                                                                                                                  cgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGAc             <  2:1434369/148‑1 (MQ=255)
                                                                                                                                  cgGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGAc             <  1:1026569/148‑1 (MQ=255)
                                                                                                                                    gCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACAACGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACAt           >  1:186236/1‑148 (MQ=255)
                                                                                                                                    gCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACAAAGCCACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACAt           >  1:185933/1‑148 (MQ=255)
                                                                                                                                      attcattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGNCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATcc         <  1:1802942/148‑1 (MQ=255)
                                                                                                                                      attcattTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATcc         >  2:631971/1‑148 (MQ=255)
                                                                                                                                        tcattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcag       <  1:309295/148‑1 (MQ=255)
                                                                                                                                           ttttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagcag    >  2:2835801/1‑148 (MQ=255)
                                                                                                                                           ttttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagcag    <  1:275141/148‑1 (MQ=255)
                                                                                                                                           ttttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCACCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagcag    >  2:403939/1‑148 (MQ=255)
                                                                                                                                            tttACTAGTGGGAGTAAACCAGTCGGTTGTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATATGTCGCGCCCTTGTCGAGTACTTCCCACCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGAACCACATCcagcag  >  2:3347086/1‑147 (MQ=255)
                                                                                                                                             ttACTAGTGGGAGTAAACCAGTCGGTTTTTTTAAAAGGAAAATAGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCGTTCGGCAAGTTCCGGCATCAGGCGATAGGTATATGGATCGACAAACAGCAGTc  >  1:1684665/1‑148 (MQ=255)
                                                                                                                                             ttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTc  >  2:2370700/1‑148 (MQ=255)
                                                                                                                                             ttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACAACGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGATTACTTCCCACCTTTAGGCAAGTTCCGGCATCAGGCGATAAGTATAGGGATCGACAACCAGCAGTc  >  1:1698419/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTC  >  CP009273/1350877‑1351165

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