Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,751,5370AC47.8% 82.2 / 15.1 67V111G (GTT→GGT) yiaTputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/11);  new base C (0/32);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 4.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGG  >  NC_000913/3751390‑3751684
                                                                                                                                                       |                                                                                                                                                    
aCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑ccc                                                                                                                                                  >  2:541376/1‑151 (MQ=255)
  gCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑ccccc                                                                                                                                                >  1:946857/1‑151 (MQ=255)
   cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCg                                                                                                                                               <  1:608495/151‑1 (MQ=255)
   cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCg                                                                                                                                               >  2:709675/1‑151 (MQ=255)
      aGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGcc                                                                                                                                             >  2:26164/1‑150 (MQ=255)
        cgcAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCAATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCAt                                                                                                                                           >  1:396763/1‑150 (MQ=255)
               gtCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTAc                                                                                                                                   >  1:170323/1‑151 (MQ=255)
                      aCAGACGCCCTTTCTGCTTTTGGTGGAATCCCGAACGCCCCCCCCCCCGGCCGTTGCTGTTACCGAGAC‑CGCCCACGGCTGCAGCGCTTTTCACGGTCCCCCAC‑CG‑CGCG‑TGTTGACACCGGGCAC‑CCCCCGCCATAGCCGTACTATCAc                                                                                                                             <  1:769161/150‑1 (MQ=255)
                      aCAGACGCCCTACCTGCATTTTGGGAAATACCGATGGCCCCCCCCCCCAGCCTTTGCTTTTATCCAGAC‑CGTCCGCAGCGGCAGAGGCTTTCCCGGTTCCCCAC‑CG‑CTCG‑GGGTGATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCAc                                                                                                                             <  1:523222/150‑1 (MQ=255)
                       cagagacccTATCTTCATTTTGGGTAATCCCGATATCCCCCCCCCCCGCCCTTTCCTCTTTCCCAGAC‑CCTCCGCGCCGGCAGAGGCTTTC‑CGGTTCCCCCT‑CC‑CCCGTTGGTGACACCAGACAC‑CCCCCGCCATAGCCGTACTATCTCg                                                                                                                            <  1:1054673/143‑1 (MQ=255)
                          acccccTACCGGCATTTTGTGAAATCCGGTAACCCCCCCCACCCAGCGTTTTCTTTTATCCAAAT‑GGTCGGCAGCTCGAGAGCTTTCCCCGGTTCCCCAA‑CG‑CTCG‑TGTGTATGCCGGCCAA‑CCCCCGCCATAGCCGTACTATCACGCtt                                                                                                                         <  2:170901/146‑1 (MQ=255)
                           cgacCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTAt                                                                                                                       >  1:237512/1‑151 (MQ=255)
                             acCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATc                                                                                                                      >  1:5972/1‑150 (MQ=255)
                              ccTTTCTCCTTTTTGTGAAATCCCGAAGGCCCCCCCCCCCCCCCGTTGCTGTTACCAAAAC‑CGCCCGCGGCGGCGGAGCTTTTCCCGTTCCCCCAC‑CC‑CTGG‑TGGTAATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCa                                                                                                                     <  1:1076334/150‑1 (MQ=255)
                                  tCTGCATTTTGTGAAATCCCGATAGCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGc                                                                                                                  <  2:415086/149‑1 (MQ=255)
                                    tGCTTTTTGTGAAATCCCGATAGCCCCCCCCCCCAGCCTTTGCTTTTACCCAGAC‑CGCCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTAATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATc                                                                                                                      <  1:431252/143‑1 (MQ=255)
                                       aTTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                            >  2:1206913/1‑150 (MQ=255)
                                        ttttGTGGAATACGGAAACCCCCCCCCCCCAGCCGTTGCTGTTACCCACAC‑CGTCCGCAGCTGGAGAGCCTTTCACGTTTCCCCAA‑CA‑CTCG‑TGGGGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                            <  2:15793/149‑1 (MQ=255)
                                         tttGTGAAATCCCGACAGCCCCCCCCCCCAGCCGTTCCTGTTATCCAAAC‑CGCCCGCGCCGGCAGAGGTTTTCAGGGTCCCCCAA‑CC‑CTCG‑TGGTAATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGt                                                                                                               <  1:996650/145‑1 (MQ=255)
                                              ggAATCCCGATACCCCCCCCCCCCGGCCGTTGCTGTTACCCAGAG‑CGTCCGCAGCTGCAGGGGCTTTCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATCCCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCAt                                                                                                     <  2:828554/150‑1 (MQ=255)
                                                aaTACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAAt                                                                                                  >  1:1075744/1‑151 (MQ=255)
                                                 aTACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATc                                                                                                 >  2:677266/1‑151 (MQ=255)
                                                 aTACCAATACCCCCCCCCCCCAGCCGTTGCTGTTACCCAGAC‑CGTCCGCAGCTGCAGAGGCTTTCACGTTTCCCCAC‑CC‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAAt                                                                                                  <  2:1014451/150‑1 (MQ=255)
                                                    ccgatgcccccccccccCAGCCGTCGCTGTTATCCAGAA‑CGTCCGCAGCTGCGGTGCTTTTCCCGGTTCCCCAC‑CG‑CTCG‑TGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCg                                                                                                <  2:508338/140‑1 (MQ=255)
                                                    ccGATAGCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                              <  2:772263/151‑1 (MQ=255)
                                                    ccGATACCCCCCCCCCCCAGCCGTTGCTGTTACCCAGAC‑CGCCCGCAGCTGCAGGGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                              <  1:233976/151‑1 (MQ=255)
                                                    ccGAGAGCTCCCCCCCCCACCCGTTGCTGTTACCCAGAA‑CGCCCGCGGCTGCAGGGGCTTTCACGGTCCCCCAC‑CG‑CTCG‑TGGTGATACCGGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                              <  1:476764/151‑1 (MQ=255)
                                                    cagATACCCCCCCCCCCCACCCGTCGCTGTTATCCAGAT‑CGTCCGCACCTGCGGGGGCTTCCACGGTCCCCCAC‑CG‑CTCG‑TGGTGATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGt                                                                                               <  2:1083719/148‑1 (MQ=255)
                                                       atacccccccccccCACCCGTTGCTGTTACCCAGAA‑CGTCCGCGGCGGCAGAGGTTTCCACGTTCCCCCAC‑CG‑CTCG‑TGGTGATCCCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCAtt                                                                                           <  1:1137341/145‑1 (MQ=255)
                                                             ccccccccCAGCCGTTGCTGTTATCCAGAC‑CGTCCGCAGCTGCAGAGCCTTTCCCGGTCCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                     <  2:350345/151‑1 (MQ=255)
                                                              ccccACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                     >  2:380293/1‑150 (MQ=255)
                                                                 cccccaCCCGTTCCTGTTACCAAGAC‑CGGCCCCAGCGGCAGGGGCTTTCGCGGTCCCCAAC‑CG‑CTGG‑TGTTGACCCCGGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAggg                                                                                 <  1:376255/149‑1 (MQ=255)
                                                                   cccaGCCGTTGTTGTTATCCAGAC‑CTTCCGCAGCTGCAGAGGCTTTCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGt                                                                                <  2:552111/150‑1 (MQ=255)
                                                                        cctttGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGGGGCTTTCACGTTTCCCCAA‑CG‑CTCG‑TGGTGATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAAtt                                                                          <  2:1142754/148‑1 (MQ=255)
                                                                        ccGTGGCTGTTGTCAAAAT‑CGTCCGCGGGGGCGGAGGCTTTCCGGGTCCCCCAA‑CG‑CTGG‑TGGGGATACCGGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAAtt                                                                          <  1:241452/151‑1 (MQ=255)
                                                                           ttGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                        >  2:1006868/1‑150 (MQ=255)
                                                                           ttGCTGTTACCCAGAC‑CGCCCGCAGCTGCAGAGGCTTTCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATCCCGGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTc                                                                         <  1:162085/149‑1 (MQ=255)
                                                                            tGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                        >  2:929231/1‑149 (MQ=255)
                                                                             tcttttATTCAGAACCGTCCCGAGCGGAGGAGTCTTTCAC‑GTTCCCCCA‑CGTCCCG‑TGGTAATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTc                                                                         <  1:351228/144‑1 (MQ=255)
                                                                               tGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTGTCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCgg                                                                   >  1:37002/1‑151 (MQ=255)
                                                                                      ccaaacctcccccAGCGGCAGAGGTTTTCCCGTTCCCCCAA‑CC‑CTCG‑TGGTGACCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCa                                                             <  1:588120/139‑1 (MQ=255)
                                                                                       aGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAg                                                            <  2:780641/150‑1 (MQ=255)
                                                                                       aGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAgg                                                           >  1:976319/1‑151 (MQ=255)
                                                                                            cGTCCGCAGCTGCGGAGGCTTTCGCGTTCCCCCAC‑CG‑CTCG‑TGGTGATCCCAGCCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGAt                                                         <  1:8610/149‑1 (MQ=255)
                                                                                              tCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCgg                                                     >  1:755002/1‑151 (MQ=255)
                                                                                               ccccggCTGCGGAGCTTTTCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATACC‑GGCAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCg                                                      <  2:286841/144‑1 (MQ=255)
                                                                                                cGCAGCTGCGGAGGTTTCCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTc                                                   <  1:470293/151‑1 (MQ=255)
                                                                                                cGCAGCTCCGGAGGCTTTCCCGTTCCCCCAA‑CG‑CCCG‑TGTGGATACGAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATg                                                        <  2:993035/146‑1 (MQ=255)
                                                                                                     cTGCAGAGTCTTTCACGTTCCCCCAC‑CG‑CTCG‑TGGTGATACCAGCCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAg                                                 <  1:26163/148‑1 (MQ=255)
                                                                                                        cgggggCTTCCACGGTCCCCCAA‑CG‑CTCG‑TGGTGATCCCGGCCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCa                                                  <  2:698546/140‑1 (MQ=255)
                                                                                                          ggggCTTCCCCGCTCCCCCAG‑CG‑CTCG‑TGGTGATCCCAGGCAC‑CCCCCGCCATAGCCGTACTGTCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGtt                                         <  2:652423/149‑1 (MQ=255)
                                                                                                               tttCCCGTTCCCCCAC‑CG‑CTCG‑TGGGGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtt                                     <  1:576080/150‑1 (MQ=255)
                                                                                                                tctcccgTTCCCCAATCG‑CGCG‑TGGAGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTc                                        <  2:1068548/141‑1 (MQ=255)
                                                                                                                  cgcGGTTCCCCAA‑CG‑CTGG‑TGGTGATACCAGACAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtttt                                   <  1:541350/147‑1 (MQ=255)
                                                                                                                    cGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTc                               >  1:208664/1‑151 (MQ=255)
                                                                                                                     ggTCCCCCAA‑CG‑CGGG‑GGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttt                                    <  2:970526/145‑1 (MQ=255)
                                                                                                                       tcccccAC‑CG‑CTCG‑TGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCa                              <  2:170327/147‑1 (MQ=255)
                                                                                                                          cccAC‑CG‑CTCG‑GGGTGATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTg                         <  1:84358/151‑1 (MQ=255)
                                                                                                                           ccAA‑CG‑CTCG‑TGGTGATACCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGc                        <  2:367522/151‑1 (MQ=255)
                                                                                                                                 g‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                     >  2:271010/1‑149 (MQ=255)
                                                                                                                                 g‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                     >  2:1203475/1‑149 (MQ=255)
                                                                                                                                         ggtgATCCCAGGCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGaa                 <  2:557990/147‑1 (MQ=255)
                                                                                                                                             aTACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCt         >  2:4420/1‑151 (MQ=255)
                                                                                                                                               aCCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTa        >  1:678504/1‑150 (MQ=255)
                                                                                                                                                ccAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAAc      >  2:565804/1‑151 (MQ=255)
                                                                                                                                                   ggCAC‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGt    <  2:237518/150‑1 (MQ=255)
                                                                                                                                                    gCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTgg  >  2:415004/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAA‑CGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAA‑CG‑CTCG‑TGGTGATACCAGGCAA‑CCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGG  >  NC_000913/3751390‑3751684

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