Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,4310AC43.8% 97.8 / 23.4 65G34G (GGT→GGGlysPlysine transporter
Reads supporting (aligned to +/- strand):  ref base A (26/10);  new base C (0/28);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 4.44e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGC  >  NC_000913/2248281‑2248574
                                                                                                                                                        |                                                                                                                                                 
aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGa                                                                                                                                                   <  1:430875/151‑1 (MQ=255)
   ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                 <  2:1088514/150‑1 (MQ=255)
    gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                                              <  1:1039560/151‑1 (MQ=255)
    gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                                              >  1:792925/1‑151 (MQ=255)
     cAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                             <  1:56149/151‑1 (MQ=255)
          gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCg                                                                                                                                        >  2:1182779/1‑150 (MQ=255)
            aCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATg                                                                                                                                     <  2:1065673/151‑1 (MQ=255)
             ccGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATgg                                                                                                                                    <  1:1012264/151‑1 (MQ=255)
              cGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATgg                                                                                                                                    >  2:110463/1‑150 (MQ=255)
                gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                                                 <  1:1036376/151‑1 (MQ=255)
                 aaaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAcc                                                                                                                                <  1:763430/151‑1 (MQ=255)
                         atatatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATgg                                                                                                                        >  1:984882/1‑151 (MQ=255)
                            tatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCa                                                                                                                      >  1:69127/1‑150 (MQ=255)
                                      gtccacggggctgggTCACAAGAAAGTAACCCACCAGGCAAACCAGCAAATACGAGAGCAATCCCCCGCCCGGCCCTGCCGGAAAAATCGTGGCGCAAGGGGAAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAatca                                                                                                                  <  1:1118139/131‑1 (MQ=255)
                                      gTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTc                                                                                                             >  2:905154/1‑149 (MQ=255)
                                            cGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcg                                                                                                     >  1:86653/1‑151 (MQ=255)
                                              agaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                   >  1:143575/1‑151 (MQ=255)
                                                   ggTCATCAGAAAGAAAACCATCAGGCAAATCAGCATATACGAGAGCAACCCCCCGCCCGGCCCTGCCTGAAAATTCGTGGCGCAAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                              <  2:436299/151‑1 (MQ=255)
                                                    gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                              >  2:189512/1‑150 (MQ=255)
                                                       atcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                          >  1:1031161/1‑151 (MQ=255)
                                                       acaagaaaGAAACCCACCAGCCCAACCAGCATATACGAGACAAATCCCCCCCCCGGCCCTGCCTGAAAAATCGTGGCGCAAGAGGCAAAAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTa                                                                                             <  1:65591/142‑1 (MQ=255)
                                                           ggAAGAAAACCATCAGCCAAATCAGCATATACGAGAGAAACCCCCCGCCCGGCCCTGCCTGAAAAATCGTTGCGCAAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCa                                                                                        <  2:629611/149‑1 (MQ=255)
                                                             aaacaacccaccaGGCAAATCACAATAAACGAGACAAACCCCCCGCCCGCCCCTGCGGAAAAAATCGTGGCGCGAGAGGCAAAAAAAAGCCCGGTACCAATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                     <  1:824694/143‑1 (MQ=255)
                                                                 aaacctgacgccaaaTCAGAATAACCGAGAGCAACCCCCCGCCCGGCCCTGCCTGAAAAATCTTTCCGCGAGAGCCAACAAAAAGCCCTATCCCGATGGATCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                 <  2:570819/137‑1 (MQ=255)
                                                                 aaaCCACCAGCCAAATCAGCATATACGAGAGAAACCCCCCCCCGGGCCTTGCCTGAAAAACCGTGGCGCCAAGGAAAACAAAAGCCCTTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGt                                                                                  <  2:238888/149‑1 (MQ=255)
                                                                          gccaaaTCAGCATATAGGAGAGCATCCCCCCGCCGGGCCCTGCCTGAGAATTGTTGGGGCCGGAGGCAAAAAAAAGCCCGGCCCTGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGcc                                                                             <  2:1129546/141‑1 (MQ=255)
                                                                          gaccaaaaagaatatACGAGAGCAATCCCCCGCCGGGACCGCCTGGAGAATTGTTGGCCCAAGGGGCAAAAAAAAGCCCGTCCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                         <  2:1089610/138‑1 (MQ=255)
                                                                           gCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                                                      >  2:124073/1‑151 (MQ=255)
                                                                            ccAATCAGAATATACGAGAGCAATCCCCCCCCCGGCCCGGCCTGAAAAATCGTTGCGCAAGAGCCAAAAAAAAGACCTGTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                                                      <  1:1142987/150‑1 (MQ=255)
                                                                              aaccagataataCGAGAGCAATCCCCCCCCCGGCCCTGCCTGAAAATTGGTTGCCCCAGAGGAAAAAAAAGGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                   <  2:1164288/142‑1 (MQ=255)
                                                                              aaccAGCATATGGGAGAGCAACCCCCCCCCCGGCCCTGCCGGAGAAATCGTGGCGCAAGGGGAAAAAAAAAGCCCTTTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCtt                                                                     <  2:451509/146‑1 (MQ=255)
                                                                              aaccAGAATATACGAGAGAAATCCCCCCCCGGGCCCTCCTTGAAAAATCGTTGCGCGAGAGCCAAAAAAAAGCCCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                   <  2:361873/148‑1 (MQ=255)
                                                                              aaTAAGCATATAGGAAAGCAGCCCCCCGCCCGCACCAGCGTGAAAATTGGTGGGGCAAGAGAAAAAAAAAAGACCTGCCCAGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                    <  1:521892/150‑1 (MQ=255)
                                                                                 cAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgtg                                                                >  2:331125/1‑151 (MQ=255)
                                                                                      atagcaaagagcacgCCCCCGCCGGGCCTGGCCGAAAAATTGTTGGCGCCAGAGAAAAAAAAAACCCTTCTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                 <  2:1158819/132‑1 (MQ=255)
                                                                                         aCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                        >  2:818207/1‑151 (MQ=255)
                                                                                          cGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGtt                                                       >  2:785100/1‑151 (MQ=255)
                                                                                                cacccccccGCCCGGACCTGCCGGAGAAATGTCGGCGCAAGGGGCAAAAAAAGGCCCGTTCCGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                             <  1:31713/134‑1 (MQ=255)
                                                                                                caaccccccGCCGGGCCCTGCCAAAAAAATCGTTGCCCAAGAGCAAAAAAAAAGCCCTGTCCAGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                                                    <  1:233274/143‑1 (MQ=255)
                                                                                                 aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                 >  1:963138/1‑150 (MQ=255)
                                                                                                  accccccGCCCGCCCCTGCCTGAAAAATCGTCGCGCCAGAGGAAAAAAAAAGCCCTTTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                      <  1:547175/141‑1 (MQ=255)
                                                                                                    gCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                                              >  1:1107713/1‑150 (MQ=255)
                                                                                                      cccccccgCCCCTGCCGGAAAAATCGTTGCGCCAGAGCCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATa                                           <  2:416463/151‑1 (MQ=255)
                                                                                                          cccgcccctgcctaAAAAATGGTGCCGCCGAAGCCAAAAAAAAGCCCTGCCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                         <  1:1057308/143‑1 (MQ=255)
                                                                                                          cccgGCCCTGCCTGAAAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                        <  2:219148/150‑1 (MQ=255)
                                                                                                                    cctgAAAAACCTTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaaa                                     <  1:992253/143‑1 (MQ=255)
                                                                                                                     ctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                             >  2:145797/1‑150 (MQ=255)
                                                                                                                        gagaAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATc                          >  1:408047/1‑151 (MQ=255)
                                                                                                                        gaaaaaTCTTTGCACCAGAGCCAAAAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                              <  2:567710/144‑1 (MQ=255)
                                                                                                                         aaaaaTCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                           <  2:736901/147‑1 (MQ=255)
                                                                                                                         aaaaaTCGTTGCGCCAGAGAAAAAAAAAAGCCCTGTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                             <  1:1016826/145‑1 (MQ=255)
                                                                                                                         aaaaaCCGTGGCGCCAGAGGAAAAAAAAAGCCCGGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCa                            <  1:238714/146‑1 (MQ=255)
                                                                                                                          gaAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                        >  1:787269/1‑151 (MQ=255)
                                                                                                                           aaaGCGTTGCGCCAGAGCAAAAAAAAAGACCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAtat                         <  1:1141827/149‑3 (MQ=255)
                                                                                                                           aaaGCGTTGCGCCAGAGAAAACAAAAAGCCCTGTCCAGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                        <  2:193385/150‑1 (MQ=255)
                                                                                                                            aaTCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                       >  1:953537/1‑150 (MQ=255)
                                                                                                                                 ttGCGCAAGAGGCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaa                    <  2:63661/148‑1 (MQ=255)
                                                                                                                                  tggcgcCAGAGGAAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAAt                  <  1:230280/148‑1 (MQ=255)
                                                                                                                                     cgcCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAg               >  2:207296/1‑150 (MQ=255)
                                                                                                                                       ccAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg            >  1:970646/1‑151 (MQ=255)
                                                                                                                                       ccAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg            >  2:455553/1‑151 (MQ=255)
                                                                                                                                              caacaaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt     >  1:551697/1‑151 (MQ=255)
                                                                                                                                               aacaaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt     >  1:18317/1‑150 (MQ=255)
                                                                                                                                               aaaaaaaaGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt       <  2:692306/145‑1 (MQ=255)
                                                                                                                                                aaaaaaaGCCCTGCCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgc             <  1:1037363/139‑1 (MQ=255)
                                                                                                                                                   aaaaGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc  >  1:1177791/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
AGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGC  >  NC_000913/2248281‑2248574

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