Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG39.0% 92.2 / 14.1 60I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (9/27);  new base G (23/0);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 2.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456956‑4457254
                                                                                                                                                      |                                                                                                                                                      
atatGCCGTGGCAGGTGGCAGAAGTGCCAGAAGCCTGCGTTCCGCCTGCACACTGCAGTGTCGGTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGGCGGACGACCGTTTTCTGGTGATCAGCAACGGCGG‑ATTTGGTGGGaa                                                                                                                                                       <  2:157203/151‑2 (MQ=255)
  atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATc                                                                                                                                                      <  2:1139474/150‑1 (MQ=255)
     ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCAt                                                                                                                                                   <  2:685387/150‑1 (MQ=255)
         ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGaa                                                                                                                                              <  2:761611/151‑1 (MQ=255)
           caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTGCCTCAGaaa                                                                                                                                             >  1:1330123/1‑150 (MQ=255)
            aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAAc                                                                                                                                            <  1:926149/150‑1 (MQ=255)
             ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAActg                                                                                                                                          <  1:1202626/151‑1 (MQ=255)
              gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATGCGCTGTCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAActgc                                                                                                                                         <  2:722979/151‑1 (MQ=255)
                   agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCCTCAGAAAATGCTgg                                                                                                                                      >  1:1047734/1‑149 (MQ=255)
                    gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGGGGGGTGCATAGGGAAccgctg                                                                                                                                       >  1:1019994/1‑147 (MQ=255)
                           cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGTGGG‑GCCCTCAAAAAATGCTGGATGGGc                                                                                                                                >  1:990500/1‑147 (MQ=255)
                               agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaag                                                                                                                        <  2:287036/151‑1 (MQ=255)
                               agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaag                                                                                                                        <  1:991364/151‑1 (MQ=255)
                                gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGGCCCACCGAAACTGCtgggtggggggg                                                                                                                             >  2:1110464/1‑141 (MQ=255)
                                 ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa                                                                                                                      <  2:819944/151‑1 (MQ=255)
                                    gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAgg                                                                                                                    <  1:583599/150‑1 (MQ=255)
                                      gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCgg                                                                                                                 <  1:869589/151‑1 (MQ=255)
                                      gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCAACAGAAAACGGTGGTTGGGGTGGGGGAGAAAGCgg                                                                                                                 >  1:32262/1‑151 (MQ=255)
                                         cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                               <  1:738259/150‑1 (MQ=255)
                                           gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                            <  1:292095/151‑1 (MQ=255)
                                               gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg                                                                                                         <  2:342770/150‑1 (MQ=255)
                                                cacaCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGGGGAGAAGGCgggaacc                                                                                                            >  1:123970/1‑141 (MQ=255)
                                                      ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCAACCGAAACTGCTGGGTGGGCGGGCGGAAGAGGGGGGAGCCGCGCAgg                                                                                                     >  2:82695/1‑146 (MQ=255)
                                                        aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCg                                                                                              <  1:888461/151‑1 (MQ=255)
                                                             cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCtt                                                                                         <  1:1018789/151‑1 (MQ=255)
                                                             cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGGGGGGGCCATCAGAAAATGCTGGGTGGGCTGGGGAAGAAAGGGGGAGCCGCGCAGGAA‑ATCGGccct                                                                                         >  2:1106778/1‑149 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGGTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTa                                                                                        <  1:504060/151‑1 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCtt                                                                                         <  2:367792/150‑1 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTa                                                                                        <  2:446161/151‑1 (MQ=255)
                                                                   ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAA‑TTCGGCCTTAGcc                                                                                     >  1:1301160/1‑149 (MQ=255)
                                                                       aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGGAA‑TTCGGCCTTAGCCAg                                                                                   >  1:278981/1‑147 (MQ=255)
                                                                        cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGGAA‑TTCGGCCTTAGCCAGata                                                                                >  1:1006962/1‑149 (MQ=255)
                                                                        cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGTCCTTAGCCAGAtata                                                                              <  2:1110130/151‑1 (MQ=255)
                                                                          cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGAtata                                                                              >  1:1148734/1‑149 (MQ=255)
                                                                            ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGGAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGt                                                                            <  2:1275048/149‑1 (MQ=255)
                                                                             gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCAGAAACTGGTGGATGGGGTGGCGAAGAAAGCGGGAGCCGCGCCGGAA‑ATCGGGCCTAGCCAGATATAGTa                                                                           >  2:490155/1‑149 (MQ=255)
                                                                              cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAA‑TTCGGCCTTAGCCAGATATAGTa                                                                           >  1:1098443/1‑148 (MQ=255)
                                                                               aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGGCGGCGGAGAAGGGGGGAGCCGCGCAGTAA‑TTTGGCCTTAGGCAGGTATAGTATGc                                                                        >  2:774140/1‑150 (MQ=255)
                                                                                 atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTa                                                                     <  1:1212304/151‑1 (MQ=255)
                                                                                  tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGc                                                                        >  1:1078437/1‑147 (MQ=255)
                                                                                   atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑ATTTGGTGGGGGTCATCAAAAACTGCTCGGTGGGGCGGCGCAGAAAGGGGAAGCCGCGCAGTAA‑TTTCGGCTTAGCCAGGTATAatatatgca                                                                     >  1:163878/1‑144 (MQ=255)
                                                                                       tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAAAAGGCGGGAACCGCGCAGTAA‑TTCGGCCCTTGCCAGAAATAGGATGCg                                                                       >  1:1317855/1‑143 (MQ=255)
                                                                                           gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATTCGTAGGCCgg                                                               >  1:418606/1‑147 (MQ=255)
                                                                                              aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAgg                                                         >  1:843469/1‑150 (MQ=255)
                                                                                                 cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCggt                                                     >  1:277622/1‑149 (MQ=255)
                                                                                                  gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                    <  2:344792/151‑1 (MQ=255)
                                                                                                      aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGGGGGGTCCATCAGAAACTGCTGGATGGGGCGGCGAAGAAAGGGGGAGCCGCGCAGGAA‑ATTGGCCTTTGCCAGGTATAGTATGCGTAGGGCGGGTAAGGCggt                                                     >  2:973350/1‑144 (MQ=255)
                                                                                                          tGGCGACCATATTCTGGTGATGAGCAACGGCGG‑CTTTGGGGGGATCCATAAAAAACTGCTGGGGGGGCTGGCGAAGAAGGCGGAAGCCGGGCAGGAA‑ATCGGCCTTAGCCAGATATAGTAGGCGTAGGCCGGAAAAGGCGGTCAcgccgc                                             >  1:9031/1‑150 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                    >  1:265012/1‑151 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                      <  1:1160667/149‑1 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                    <  2:460971/151‑1 (MQ=255)
                                                                                                                    atTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                                  <  1:221072/151‑1 (MQ=255)
                                                                                                                    atTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTCCATCCGAAAATGGTTGGTGGGCCGGGGGAGAAGGCGGGAGGCGCCCCGTAA‑TTTGGCCTTAGCCAGATTAAGTATGGGTTGGGCGGGAAAAGCGTTCACCCCGCATCCGGCCttt                                  >  2:1222928/1‑151 (MQ=255)
                                                                                                                       cTGGTGATGAGCAACGGCGG‑TTTTGGGGGGGTTCATCAGAAAATGCTGGATGGGGTGGCGGAGAAAGGGGGAGCCGCGCCGTAATTTCGGCC‑TAGCCCGATATAGTATGCGGAGGCCGGGGAAGGGGTTCACGCCCCATCCGGCAttt                                  >  2:356628/1‑148 (MQ=255)
                                                                                                                        tGGTGATGAGCAACGGCGG‑TTTTTGTGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCATTTGAg                               >  1:1041355/1‑150 (MQ=255)
                                                                                                                                         gg‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc             >  1:912157/1‑151 (MQ=255)
                                                                                                                                         gg‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc             <  1:85909/151‑1 (MQ=255)
                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  2:1213034/1‑151 (MQ=255)
                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  1:514151/1‑151 (MQ=255)
                                                                                                                                             tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCt          >  2:7874/1‑151 (MQ=255)
                                                                                                                                               tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg        >  1:16611/1‑151 (MQ=255)
                                                                                                                                                      aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGt  <  1:431602/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAA‑TTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456956‑4457254

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