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

CGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456962‑4457255
                                                                                                                                               |                                                                                                                                                       
cGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGTTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCa                                                                                                                                                  <  2:42318/151‑1 (MQ=255)
cGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTcc                                                                                                                                                      >  2:379579/1‑147 (MQ=255)
   ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGaa                                                                                                                                               <  2:126203/151‑1 (MQ=255)
     caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGAGATGAGCAACGGCGGCTTTCGAGGGACGCATCACaaa                                                                                                                                              >  2:633051/1‑150 (MQ=255)
       ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGaaa                                                                                                                                              >  2:720536/1‑148 (MQ=255)
       ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct                                                                                                                                            <  1:211783/150‑1 (MQ=255)
             agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGGACATCCGAAActgtgg                                                                                                                                        >  2:627355/1‑145 (MQ=255)
              gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                                    <  1:822506/151‑1 (MQ=255)
               aagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                                    <  2:647397/150‑1 (MQ=255)
                 gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATggg                                                                                                                                  >  1:550970/1‑150 (MQ=255)
                   ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCCTCCGAAAATGCtgggtggg                                                                                                                                  >  1:727539/1‑148 (MQ=255)
                    gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTgg                                                                                                                              <  1:557040/151‑1 (MQ=255)
                     cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAACGCTGGGTGGGCTg                                                                                                                               >  2:297931/1‑149 (MQ=255)
                      agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGACCTCAGAAAATGCTgggagggggggcg                                                                                                                            >  1:98223/1‑140 (MQ=255)
                             tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCATCAGAAAATGCTGGATGGGGTGGCggagaa                                                                                                                      >  1:525001/1‑149 (MQ=255)
                                gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGAGCCATCAGAAACTGGTGGATGGGTGGGCGAAGAAGGCgg                                                                                                                 >  2:743477/1‑151 (MQ=255)
                                  tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAGGCgg                                                                                                                 >  1:100913/1‑149 (MQ=255)
                                   cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCgg                                                                                                                 >  1:70001/1‑148 (MQ=255)
                                    aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTTTATCTGAAACTGCTGGGTGGGGTGGGGCAGAAGGGGGGAg                                                                                                              >  1:95392/1‑150 (MQ=255)
                                       cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGcc                                                                                                            >  1:554861/1‑149 (MQ=255)
                                            cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTGCATAAGGAACTGCTGGGTGGGGTGGCGAAGAAGGGGGGAGCCGCGCAg                                                                                                      >  2:801501/1‑150 (MQ=255)
                                                ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTaa                                                                                                   <  1:720528/149‑1 (MQ=255)
                                                 gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                                <  2:635949/151‑1 (MQ=255)
                                                     ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCAACAGAAACTGCTGGAAGGGCTGGCGGAAGAAGCGGGAACCGCGCAGGAATTCgggc                                                                                            >  2:341433/1‑149 (MQ=255)
                                                        gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                         <  2:509375/151‑1 (MQ=255)
                                                        gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAAGGGGGAGCCGCGCAGGAATTTGGCCtt                                                                                          >  1:806409/1‑150 (MQ=255)
                                                         aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGGTGGGGGAGAAAGCGGGAGCCGCGCCGGAAATTGGCCtt                                                                                          >  1:938378/1‑149 (MQ=255)
                                                         aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTGGTGGGGTCCATCAGAAACAGCTGGGTGGGGGGGGGGAGAAGGCGGAAGCCGCGCAGTAATTCGGCCCTa                                                                                         >  1:380207/1‑150 (MQ=255)
                                                                tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTCGGCCCTAGcc                                                                                      >  1:112914/1‑146 (MQ=255)
                                                                 aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGatat                                                                                <  2:494027/151‑1 (MQ=255)
                                                                    cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATAAAAAACGGCTGGGTTGGCTGGCCCAGAAGGCGGATGCCGCGCCGTAATTAGGCCTTTGCCAGAGATAGt                                                                             >  2:772910/1‑151 (MQ=255)
                                                                     tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCATCAGAAAATGATGGATGGGCTGGGGAAGAAGGGGGAAGCCGCGCAGGAATTCGGCCTTAGTCCGAtata                                                                               >  1:311825/1‑148 (MQ=255)
                                                                      ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGGGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGAtata                                                                               >  2:875327/1‑147 (MQ=255)
                                                                         aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGATGGGGTGGGGGAGAAAGCGGAAGCCGCGCCGTAATTCGGCCTTTGCCAGATATAATATGCg                                                                        >  1:122969/1‑151 (MQ=255)
                                                                           atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGAGCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGt                                                                       >  1:670114/1‑150 (MQ=255)
                                                                             atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGACCATCAGAAACTGCTGGAAGGGGTGGGGGAGAAGGGCGGAGGCGCGCAGGAAATTCGCCTTAGCCAGATATAGTATGCGTggg                                                                    >  1:373272/1‑151 (MQ=255)
                                                                               ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGc                                                                   <  2:220414/150‑1 (MQ=255)
                                                                                  ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                                <  1:210334/150‑1 (MQ=255)
                                                                                  ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGa                                                               <  2:431835/151‑1 (MQ=255)
                                                                                    tgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAACAGCTGGGTGGGGTGGGGAAAAAAGCGGGAGCCGCGCAGTAAATCCGCCTTAGCCCGATATAGTATGCGCAgggccg                                                                >  1:459077/1‑144 (MQ=255)
                                                                                     gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATaa                                                            <  2:293393/151‑1 (MQ=255)
                                                                                           cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAAAAAATGCTGGATGGGCTGGCGAAAAAGGCGGGAGCCGCGCAGGAAATCGGCCCTTGCCCGATATAGTATGCGTAGGGCGGGAAAgg                                                          >  1:962394/1‑147 (MQ=255)
                                                                                           cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGAAGCCGCGCAGGAAATCGGCCTTAGGCCGATATAGTATGCGTAGGCCGGGTaaaggggt                                                      >  1:722068/1‑145 (MQ=255)
                                                                                                aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgc                                                 <  2:223525/151‑1 (MQ=255)
                                                                                                  ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccg                                               <  2:870458/151‑1 (MQ=255)
                                                                                                  ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGAACCCTTAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGGAGCCGCGCCGTAATTCGGCCTTTGGCAGATATAGTATGAGGAGGCCGGGTAAAGCGGTCAcgcc                                                >  2:329716/1‑150 (MQ=255)
                                                                                                     ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                            <  1:286492/151‑1 (MQ=255)
                                                                                                       cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                          <  1:713393/151‑1 (MQ=255)
                                                                                                          ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                       <  1:568177/151‑1 (MQ=255)
                                                                                                            atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                     <  2:875849/151‑1 (MQ=255)
                                                                                                              atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                                    <  1:316107/150‑1 (MQ=255)
                                                                                                                           gCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                      <  1:610481/151‑1 (MQ=255)
                                                                                                                             aaCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGa                    <  2:64871/151‑1 (MQ=255)
                                                                                                                             aaCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGa                    <  2:827884/151‑1 (MQ=255)
                                                                                                                              aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGa                    <  1:576148/150‑1 (MQ=255)
                                                                                                                               cggcggTTTTTGGGGGATCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCCCATCCGGCATTTTAGCATAGTGCCTGAt                   >  1:800702/1‑150 (MQ=255)
                                                                                                                                ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGc                 <  1:455805/151‑1 (MQ=255)
                                                                                                                                 gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg                >  2:161373/1‑151 (MQ=255)
                                                                                                                                         ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg         >  2:96637/1‑150 (MQ=255)
                                                                                                                                               aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTc  >  2:949956/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
CGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456962‑4457255

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