Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,350,891 G→T 100% S159R (AGC→AGA)  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,350,8910GT97.4% 594.7 / 6.1 190S159R (AGC→AGA) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (96/89);  total (100/90)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 194 total aligned reads displayed.

GTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAAC  >  CP009273/1350746‑1351035
                                                                                                                                                 |                                                                                                                                                
gTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAATAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCCGTTTGCGGActct                                                                                                                                                >  1:695038/1‑148 (MQ=255)
gTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCt                                                                                                                                                >  2:104578/1‑148 (MQ=255)
 ttGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCtt                                                                                                                                               <  1:393103/148‑1 (MQ=255)
  tGCTCCTGGCGATGTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCttt                                                                                                                                              >  1:3245328/1‑148 (MQ=255)
   gCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCtttt                                                                                                                                             >  1:3258411/1‑148 (MQ=255)
    cTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTc                                                                                                                                            <  2:365658/148‑1 (MQ=255)
        tGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGt                                                                                                                                        >  1:2483707/1‑148 (MQ=255)
        tGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGt                                                                                                                                        >  1:2689012/1‑148 (MQ=255)
         ggCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGtt                                                                                                                                       <  1:2583838/148‑1 (MQ=255)
         ggCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGtt                                                                                                                                       <  1:2020130/148‑1 (MQ=255)
          gCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTa                                                                                                                                      <  2:1205691/148‑1 (MQ=255)
          gCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTa                                                                                                                                      >  2:1374528/1‑148 (MQ=255)
          gCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTa                                                                                                                                      <  2:922911/148‑1 (MQ=255)
          gCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGGTTGCGGACTCTTTTCACGTTa                                                                                                                                      >  2:887014/1‑148 (MQ=255)
             aTCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCg                                                                                                                                   >  2:2297614/1‑148 (MQ=255)
              tcttctTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCgg                                                                                                                                  >  2:1886021/1‑148 (MQ=255)
                 tctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaa                                                                                                                               <  2:1427140/148‑1 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              <  1:448706/148‑1 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              >  2:1918170/1‑148 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              >  2:1981603/1‑148 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              <  1:3049355/148‑1 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              <  1:3015718/148‑1 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              <  1:654821/148‑1 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              >  1:2191661/1‑148 (MQ=255)
                  ctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaaa                                                                                                                              >  1:1828723/1‑148 (MQ=255)
                           aCTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCt                                                                                                                     <  2:761082/148‑1 (MQ=255)
                           aCTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCt                                                                                                                     <  2:1653117/148‑1 (MQ=255)
                            cTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTg                                                                                                                    >  1:527319/1‑148 (MQ=255)
                            cTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTg                                                                                                                    >  1:601470/1‑148 (MQ=255)
                            cTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTg                                                                                                                    >  1:2901607/1‑148 (MQ=255)
                            cTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTg                                                                                                                    >  1:96234/1‑148 (MQ=255)
                               ccGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCg                                                                                                                 >  2:2902120/1‑148 (MQ=255)
                                  ggCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAg                                                                                                              <  1:63965/148‑1 (MQ=255)
                                   gCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGt                                                                                                             >  2:727121/1‑148 (MQ=255)
                                      tatTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAgg                                                                                                          <  1:2935938/148‑1 (MQ=255)
                                        ttCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAgggg                                                                                                        >  2:1170633/1‑148 (MQ=255)
                                         tCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGa                                                                                                       >  2:1627421/1‑148 (MQ=255)
                                           tGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAg                                                                                                     >  1:152389/1‑148 (MQ=255)
                                           tGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAg                                                                                                     >  1:2354162/1‑148 (MQ=255)
                                            gCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGt                                                                                                    >  2:854523/1‑148 (MQ=255)
                                              cGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtg                                                                                                  <  2:168767/148‑1 (MQ=255)
                                              cGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtg                                                                                                  <  2:253519/148‑1 (MQ=255)
                                              cGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtg                                                                                                  >  2:2818699/1‑148 (MQ=255)
                                                 cATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtgctg                                                                                               >  2:3107732/1‑148 (MQ=255)
                                                   tGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGcc                                                                                             >  2:2481435/1‑148 (MQ=255)
                                                    gACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCg                                                                                            <  1:104578/148‑1 (MQ=255)
                                                     aCCGAAGCAGAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGt                                                                                           <  2:2920657/148‑1 (MQ=255)
                                                       cGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTg                                                                                         <  2:838251/148‑1 (MQ=255)
                                                         aaGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGAc                                                                                       >  1:497342/1‑148 (MQ=255)
                                                         aaGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGAc                                                                                       >  1:3219563/1‑148 (MQ=255)
                                                           gCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGt                                                                                     <  2:714478/148‑1 (MQ=255)
                                                            cATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGtt                                                                                    <  1:2213000/148‑1 (MQ=255)
                                                            cATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGtt                                                                                    <  1:1384294/148‑1 (MQ=255)
                                                               aaTGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATg                                                                                 >  2:43906/1‑148 (MQ=255)
                                                                  gggTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCa                                                                              <  2:495005/148‑1 (MQ=255)
                                                                         gAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGgttg                                                                       <  1:312404/148‑1 (MQ=255)
                                                                           ggTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGgttgtt                                                                     >  2:2555998/1‑148 (MQ=255)
                                                                           ggTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGgttgtt                                                                     >  2:793772/1‑148 (MQ=255)
                                                                            gTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGttgttg                                                                    >  1:747389/1‑148 (MQ=255)
                                                                              gCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCg                                                                  <  1:2100436/148‑1 (MQ=255)
                                                                                 acaAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCCCGGGTTGTTGCGGTc                                                               >  1:352443/1‑148 (MQ=255)
                                                                                   aaaaaaGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCaa                                                             <  2:2211668/148‑1 (MQ=255)
                                                                                     aaaaGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCaaaa                                                           <  1:218716/148‑1 (MQ=255)
                                                                                      aaaGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCaaaaa                                                          <  1:431270/148‑1 (MQ=255)
                                                                                       aagaagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAAt                                                         >  1:2971816/1‑148 (MQ=255)
                                                                                        agaagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAAtt                                                        >  1:2882568/1‑148 (MQ=255)
                                                                                        agaagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAAtt                                                        >  2:718034/1‑148 (MQ=255)
                                                                                             cATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTg                                                   >  1:1785121/1‑148 (MQ=255)
                                                                                             cATCCGGCTGGGCGAGTCGGAACTAAACGGTATGATTATCCATTTTGCGGACTCTTTTCACGTTATAGGCAAATTGCAGGCTGTCTAAGTAGGGGAAGTTTCTGCCGTTGACGTTATTCCAACGGTTGTTGCGGTCAAAAATTAGCTg                                                   >  1:1787887/1‑148 (MQ=255)
                                                                                               tCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGa                                                 <  1:3380768/148‑1 (MQ=255)
                                                                                               tCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGa                                                 <  1:3499833/148‑1 (MQ=255)
                                                                                                  ggCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAgg                                              <  2:535379/148‑1 (MQ=255)
                                                                                                             tCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACcaca                                   <  1:2043217/148‑1 (MQ=255)
                                                                                                                gAACTCAACGGTATGATTATCCAGTTTGCGGACTATTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGATGCCGTTGACGTTATTCAAAGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACAACAtcg                                >  2:2067868/1‑148 (MQ=255)
                                                                                                                gAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcg                                >  1:1041502/1‑148 (MQ=255)
                                                                                                                     cAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGc                           <  2:3272228/148‑1 (MQ=255)
                                                                                                                      aaCGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGca                          <  2:2329884/148‑1 (MQ=255)
                                                                                                                       aCGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcat                         <  1:885601/148‑1 (MQ=255)
                                                                                                                        cGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcatt                        >  1:2243559/1‑148 (MQ=255)
                                                                                                                         ggTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattc                       <  2:2837460/148‑1 (MQ=255)
                                                                                                                          gTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattca                      >  2:419561/1‑148 (MQ=255)
                                                                                                                          gTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattca                      <  2:1815164/148‑1 (MQ=255)
                                                                                                                           tATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattcat                     >  2:2364739/1‑148 (MQ=255)
                                                                                                                            aTGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattcatt                    <  2:526729/148‑1 (MQ=255)
                                                                                                                            aTGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattcatt                    >  1:2770288/1‑148 (MQ=255)
                                                                                                                             tGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCAttt                   >  1:803458/1‑148 (MQ=255)
                                                                                                                              gATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCAtttt                  >  1:3106012/1‑148 (MQ=255)
                                                                                                                              gATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCAtttt                  <  2:2717744/148‑1 (MQ=255)
                                                                                                                               aTTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTa                 <  1:1537805/148‑1 (MQ=255)
                                                                                                                                 tATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTAct               >  1:2998890/1‑147 (MQ=255)
                                                                                                                                  aTCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTActa              <  1:2812691/148‑3 (MQ=255)
                                                                                                                                      aGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGGCGGCAGTCATTTTACTAGTgg          >  2:2809179/1‑148 (MQ=255)
                                                                                                                                       gTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTggg         <  1:2526726/148‑1 (MQ=255)
                                                                                                                                       gTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTggg         >  1:2656614/1‑148 (MQ=255)
                                                                                                                                        tttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGa        >  1:129893/1‑148 (MQ=255)
                                                                                                                                         ttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAg       >  1:163208/1‑148 (MQ=255)
                                                                                                                                         ttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAg       >  1:1844130/1‑148 (MQ=255)
                                                                                                                                         ttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAg       >  1:1443706/1‑148 (MQ=255)
                                                                                                                                         ttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAg       >  1:1203471/1‑148 (MQ=255)
                                                                                                                                          tGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGt      <  2:3245328/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAAC  >  CP009273/1350746‑1351035

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