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

ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456956‑4457254
                                                                                                                                                      |                                                                                                                                                      
atatGCCGTGGCAGGTGGCAGAAGTGCCAGAAGCCTGCGTTCCGCCTGCACACTGCAGTGTCGGTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGGCGGACGACCGTTTTCTGGTGATCAGCAACGGCGGATTTGGTGGGaa                                                                                                                                                       <  2:157203/151‑2 (MQ=255)
  atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATc                                                                                                                                                      <  2:1139474/150‑1 (MQ=255)
     ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAt                                                                                                                                                   <  2:685387/150‑1 (MQ=255)
         ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGaa                                                                                                                                              <  2:761611/151‑1 (MQ=255)
           caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTGCCTCAGaaa                                                                                                                                             >  1:1330123/1‑150 (MQ=255)
            aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAc                                                                                                                                            <  1:926149/150‑1 (MQ=255)
             ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActg                                                                                                                                          <  1:1202626/151‑1 (MQ=255)
              gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATGCGCTGTCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgc                                                                                                                                         <  2:722979/151‑1 (MQ=255)
                   agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAAATGCTgg                                                                                                                                      >  1:1047734/1‑149 (MQ=255)
                    gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTGCATAGGGAAccgctg                                                                                                                                       >  1:1019994/1‑147 (MQ=255)
                           cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGTGGGGCCCTCAAAAAATGCTGGATGGGc                                                                                                                                >  1:990500/1‑147 (MQ=255)
                               agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaag                                                                                                                        <  2:287036/151‑1 (MQ=255)
                               agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaag                                                                                                                        <  1:991364/151‑1 (MQ=255)
                                gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCCACCGAAACTGCtgggtggggggg                                                                                                                             >  2:1110464/1‑141 (MQ=255)
                                 ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa                                                                                                                      <  2:819944/151‑1 (MQ=255)
                                    gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAgg                                                                                                                    <  1:583599/150‑1 (MQ=255)
                                      gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCgg                                                                                                                 <  1:869589/151‑1 (MQ=255)
                                      gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAACAGAAAACGGTGGTTGGGGTGGGGGAGAAAGCgg                                                                                                                 >  1:32262/1‑151 (MQ=255)
                                         cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                               <  1:738259/150‑1 (MQ=255)
                                           gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                            <  1:292095/151‑1 (MQ=255)
                                               gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg                                                                                                         <  2:342770/150‑1 (MQ=255)
                                                cacaCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGGGGAGAAGGCgggaacc                                                                                                            >  1:123970/1‑141 (MQ=255)
                                                      ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCAACCGAAACTGCTGGGTGGGCGGGCGGAAGAGGGGGGAGCCGCGCAgg                                                                                                     >  2:82695/1‑146 (MQ=255)
                                                        aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCg                                                                                              <  1:888461/151‑1 (MQ=255)
                                                             cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCtt                                                                                         <  1:1018789/151‑1 (MQ=255)
                                                             cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGGCCATCAGAAAATGCTGGGTGGGCTGGGGAAGAAAGGGGGAGCCGCGCAGGAAATCGGccct                                                                                         >  2:1106778/1‑149 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGGTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                        <  1:504060/151‑1 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCtt                                                                                         <  2:367792/150‑1 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                        <  2:446161/151‑1 (MQ=255)
                                                                   ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGcc                                                                                     >  1:1301160/1‑149 (MQ=255)
                                                                       aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAg                                                                                   >  1:278981/1‑147 (MQ=255)
                                                                        cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGata                                                                                >  1:1006962/1‑149 (MQ=255)
                                                                        cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGTCCTTAGCCAGAtata                                                                              <  2:1110130/151‑1 (MQ=255)
                                                                          cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGAtata                                                                              >  1:1148734/1‑149 (MQ=255)
                                                                            ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGGAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGt                                                                            <  2:1275048/149‑1 (MQ=255)
                                                                             gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGGTGGATGGGGTGGCGAAGAAAGCGGGAGCCGCGCCGGAAATCGGGCCTAGCCAGATATAGTa                                                                           >  2:490155/1‑149 (MQ=255)
                                                                              cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTa                                                                           >  1:1098443/1‑148 (MQ=255)
                                                                               aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGGCGGCGGAGAAGGGGGGAGCCGCGCAGTAATTTGGCCTTAGGCAGGTATAGTATGc                                                                        >  2:774140/1‑150 (MQ=255)
                                                                                 atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTa                                                                     <  1:1212304/151‑1 (MQ=255)
                                                                                  tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGc                                                                        >  1:1078437/1‑147 (MQ=255)
                                                                                   atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTGGTGGGGGTCATCAAAAACTGCTCGGTGGGGCGGCGCAGAAAGGGGAAGCCGCGCAGTAATTTCGGCTTAGCCAGGTATAatatatgca                                                                     >  1:163878/1‑144 (MQ=255)
                                                                                       tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAAAAGGCGGGAACCGCGCAGTAATTCGGCCCTTGCCAGAAATAGGATGCg                                                                       >  1:1317855/1‑143 (MQ=255)
                                                                                           gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATTCGTAGGCCgg                                                               >  1:418606/1‑147 (MQ=255)
                                                                                              aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAgg                                                         >  1:843469/1‑150 (MQ=255)
                                                                                                 cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCggt                                                     >  1:277622/1‑149 (MQ=255)
                                                                                                  gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                    <  2:344792/151‑1 (MQ=255)
                                                                                                      aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGATGGGGCGGCGAAGAAAGGGGGAGCCGCGCAGGAAATTGGCCTTTGCCAGGTATAGTATGCGTAGGGCGGGTAAGGCggt                                                     >  2:973350/1‑144 (MQ=255)
                                                                                                          tGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGGGGGATCCATAAAAAACTGCTGGGGGGGCTGGCGAAGAAGGCGGAAGCCGGGCAGGAAATCGGCCTTAGCCAGATATAGTAGGCGTAGGCCGGAAAAGGCGGTCAcgccgc                                             >  1:9031/1‑150 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                    >  1:265012/1‑151 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                      <  1:1160667/149‑1 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                    <  2:460971/151‑1 (MQ=255)
                                                                                                                    atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                                  <  1:221072/151‑1 (MQ=255)
                                                                                                                    atTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCCGAAAATGGTTGGTGGGCCGGGGGAGAAGGCGGGAGGCGCCCCGTAATTTGGCCTTAGCCAGATTAAGTATGGGTTGGGCGGGAAAAGCGTTCACCCCGCATCCGGCCttt                                  >  2:1222928/1‑151 (MQ=255)
                                                                                                                       cTGGTGATGAGCAACGGCGGTTTTGGGGGGGTTCATCAGAAAATGCTGGATGGGGTGGCGGAGAAAGGGGGAGCCGCGCCGTAATTTCGGCCTAGCCCGATATAGTATGCGGAGGCCGGGGAAGGGGTTCACGCCCCATCCGGCAttt                                  >  2:356628/1‑148 (MQ=255)
                                                                                                                        tGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCATTTGAg                               >  1:1041355/1‑150 (MQ=255)
                                                                                                                                         ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc             >  1:912157/1‑151 (MQ=255)
                                                                                                                                         ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc             <  1:85909/151‑1 (MQ=255)
                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  2:1213034/1‑151 (MQ=255)
                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  1:514151/1‑151 (MQ=255)
                                                                                                                                             tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCt          >  2:7874/1‑151 (MQ=255)
                                                                                                                                               tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg        >  1:16611/1‑151 (MQ=255)
                                                                                                                                                      aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGt  <  1:431602/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  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: