Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,761,9690TC51.5% 20.7 / 123.6 97intergenic (‑14/‑214)yibF/rhsAglutathione S‑transferase homolog/Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base C (22/28);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGA  >  NC_000913/3761834‑3762101
                                                                                                                                       |                                                                                                                                    
gcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGcgcg                                                                                                                                   <  1:1296849/139‑1 (MQ=255)
   aCGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAgg                                                                                                                                <  1:1005188/139‑1 (MQ=255)
    cGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGa                                                                                                                               >  1:229728/1‑139 (MQ=255)
     gCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGaa                                                                                                                              <  2:229728/139‑1 (MQ=255)
             ccGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCa                                                                                                                      <  1:50669/139‑1 (MQ=255)
                tgTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTa                                                                                                                   <  1:1206957/138‑1 (MQ=255)
                tgTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGt                                                                                                                    <  2:1474780/137‑1 (MQ=255)
                tgTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGt                                                                                                                    >  1:1474780/2‑138 (MQ=255)
                tgTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   <  1:1399859/138‑1 (MQ=255)
                tgTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   <  1:30899/138‑1 (MQ=255)
                tgTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   <  1:539806/138‑1 (MQ=255)
                  ttATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATg                                                                                                                 <  2:446843/139‑1 (MQ=255)
                          cAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTg                                                                                                         >  2:224618/1‑139 (MQ=255)
                          cAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTg                                                                                                         >  1:1382658/1‑139 (MQ=255)
                            gTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAatat                                                                                                                        <  1:1014639/122‑1 (MQ=255)
                            gTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGAt                                                                                                       >  1:1139606/1‑139 (MQ=255)
                            gTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGAGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAatat                                                                                                                        >  2:1014639/1‑122 (MQ=255)
                             ttCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATg                                                                                                      >  2:1613617/1‑139 (MQ=255)
                             ttCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATg                                                                                                      >  2:1027555/1‑139 (MQ=255)
                                  taataaATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGaaa                                                                                                 <  1:224618/139‑1 (MQ=255)
                                        aTTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGc                                                                                           <  2:1219/139‑1 (MQ=255)
                                         ttCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGaaaaa                                                                                                                           <  1:1471595/106‑1 (MQ=255)
                                         ttCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGaaaaa                                                                                                                           >  2:1471595/1‑106 (MQ=255)
                                         ttCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCa                                                                                          >  2:969153/1‑139 (MQ=255)
                                            gAAAGTTATGCCCTTTTCTAACAACAGAATAGAACGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTa                                                                                                                   <  2:1243819/111‑1 (MQ=255)
                                            gAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTa                                                                                                                   >  1:1243819/1‑111 (MQ=255)
                                             aaaGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGaaa                                                                                      <  1:1613617/139‑1 (MQ=255)
                                                           ttCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               >  2:1729378/1‑132 (MQ=255)
                                                           ttCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               <  1:1729378/132‑1 (MQ=255)
                                                                     aGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               <  1:702616/122‑1 (MQ=255)
                                                                     aGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               >  2:702616/1‑122 (MQ=255)
                                                                     aGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaa                                                              >  2:337563/1‑139 (MQ=255)
                                                                       aatggAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaat                                                            >  2:1358829/5‑139 (MQ=255)
                                                                       aaTGGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaat                                                            <  1:1358829/139‑1 (MQ=255)
                                                                       aaTAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGt                                                                                                                    >  1:287939/1‑83 (MQ=255)
                                                                       aaTAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGt                                                                                                                    <  2:287939/83‑1 (MQ=255)
                                                                       aaTAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAatat                                                                                                                        >  1:1243787/1‑79 (MQ=255)
                                                                       aaTAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAatat                                                                                                                        <  2:1243787/79‑1 (MQ=255)
                                                                           gAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaataaat                                                        <  1:969153/139‑1 (MQ=255)
                                                                            aaaGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   >  2:467471/1‑79 (MQ=255)
                                                                            aaaGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   <  1:467471/79‑1 (MQ=255)
                                                                            aaaGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGtt                                                                                 >  2:1038673/1‑113 (MQ=255)
                                                                            aaaGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGtt                                                                                 <  1:1038673/113‑1 (MQ=255)
                                                                             aaGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATca                                                      <  2:214969/139‑1 (MQ=255)
                                                                              aGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATcac                                                     <  1:337563/139‑1 (MQ=255)
                                                                                tttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACat                                                   <  1:1662605/139‑1 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTa                                                                                                                   <  2:729770/74‑1 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTa                                                                                                                   >  1:729770/1‑74 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGtt                                                                                 <  2:1686273/108‑1 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGtt                                                                                 >  1:1686273/1‑108 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGt                                                                                                                    >  1:556245/1‑73 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGt                                                                                                                    <  2:556245/73‑1 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaataa                                                          >  2:1701458/1‑131 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaataa                                                          <  1:1701458/131‑1 (MQ=255)
                                                                                 ttGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACata                                                  >  2:26314/1‑139 (MQ=255)
                                                                                  tGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACat                                                   <  1:1007245/137‑1 (MQ=255)
                                                                                  tGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACat                                                   >  2:1007245/1‑137 (MQ=255)
                                                                                   gCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGt                                                                                                                    <  1:1121633/71‑1 (MQ=255)
                                                                                   gCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGt                                                                                                                    >  2:1121633/1‑71 (MQ=255)
                                                                                   gCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   <  1:696233/72‑1 (MQ=255)
                                                                                   gCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTa                                                                                                                   >  2:696233/1‑72 (MQ=255)
                                                                                   gCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaataaat                                                        <  1:918025/131‑1 (MQ=255)
                                                                                   gCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCaaataaat                                                        >  2:918025/1‑131 (MQ=255)
                                                                                          aaCGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATg                                         <  1:1027555/139‑1 (MQ=255)
                                                                                          aaCGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATg                                         <  2:465493/139‑1 (MQ=255)
                                                                                                   gtgtAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATc                                                       >  2:10928/1‑116 (MQ=255)
                                                                                                   gtgtAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATc                                                       <  1:10928/116‑1 (MQ=255)
                                                                                                      tAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCa                                           <  1:1241934/125‑1 (MQ=255)
                                                                                                      tAGCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGAAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCa                                           >  2:1241934/1‑125 (MQ=255)
                                                                                                        gCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGaa                                                                                       <  1:1003635/79‑1 (MQ=255)
                                                                                                        gCTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATCGCCTGGTGATGGGAAAAAGGGCAGaa                                                                                       >  2:1003635/1‑79 (MQ=255)
                                                                                                         cTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGcc                                                                                                              <  1:311354/55‑1 (MQ=255)
                                                                                                         cTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGcc                                                                                                              >  2:311354/1‑55 (MQ=255)
                                                                                                         cTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGcc                                                                                                              <  1:1322788/55‑1 (MQ=255)
                                                                                                         cTACCGACGAGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGcc                                                                                                              >  2:1322788/1‑55 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGaaa                                                                                                                             <  1:1164608/36‑1 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGaaa                                                                                                                             >  2:1164608/1‑36 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAGAAGGGCa                                                                                          >  2:1127524/1‑71 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCa                                                                                          <  1:1127524/71‑1 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCt                      >  2:336498/1‑139 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGtt                                                                                 >  1:1186292/1‑80 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGtt                                                                                 <  2:1186292/80‑1 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATgg                                        >  2:404728/1‑121 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATgg                                        <  1:404728/121‑1 (MQ=255)
                                                                                                             cgacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCt                      >  2:1516312/1‑139 (MQ=255)
                                                                                                              gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATg                                         >  2:247581/1‑119 (MQ=255)
                                                                                                              gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATg                                         >  2:254372/1‑119 (MQ=255)
                                                                                                              gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATg                                         <  1:247581/119‑1 (MQ=255)
                                                                                                              gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATg                                         <  1:254372/119‑1 (MQ=255)
                                                                                                              gacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               <  1:582460/81‑1 (MQ=255)
                                                                                                              gacgaGTTTCATGCCGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               >  2:582460/1‑81 (MQ=255)
                                                                                                                            cGAGTCCTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTg       <  1:1516312/139‑1 (MQ=255)
                                                                                                                              aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTGAt     <  2:1139606/139‑1 (MQ=255)
                                                                                                                                 ccTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGa  <  1:336498/139‑1 (MQ=255)
                                                                                                                                 ccTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGa  <  2:1382658/139‑1 (MQ=255)
                                                                                                                                 ccTTTGCGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGa  <  1:26314/139‑1 (MQ=255)
                                                                                                                                    ttGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               <  2:517403/59‑1 (MQ=255)
                                                                                                                                    ttGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGa                                                                               >  1:517403/1‑59 (MQ=255)
                                                                                                                                     tgtgCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACat                                                   <  2:67289/86‑1 (MQ=255)
                                                                                                                                     tgtgCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACat                                                   >  1:67289/1‑86 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAAATAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGA  >  NC_000913/3761834‑3762101

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