Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG43.2% 108.3 / 28.8 81I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (5/41);  new base G (35/0);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.45e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAT  >  NC_000913/4456955‑4457250
                                                                                                                                                      |                                                                                                                                                 
cATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGa                                                                                                                                                   <  2:806499/151‑1 (MQ=255)
   atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATcc                                                                                                                                                <  2:180890/151‑1 (MQ=255)
            caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGaaa                                                                                                                                        >  2:321899/1‑150 (MQ=255)
             aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCCTTAGaaa                                                                                                                                        >  1:311387/1‑149 (MQ=255)
               gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCATCAGaaa                                                                                                                                        >  2:114252/1‑147 (MQ=255)
                tggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCATCAGaaa                                                                                                                                        >  1:1098848/1‑146 (MQ=255)
                 ggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgct                                                                                                                                   <  2:1167722/150‑1 (MQ=255)
                  gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCCCCAGaaaaatccgg                                                                                                                                 >  1:398156/1‑144 (MQ=255)
                    agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTTAGAAAatgctg                                                                                                                                  >  1:1116857/1‑148 (MQ=255)
                     gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                              <  1:560470/151‑1 (MQ=255)
                     gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCATCAGAAAATGCTgg                                                                                                                                 >  2:622573/1‑148 (MQ=255)
                      aagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                              <  2:5665/150‑1 (MQ=255)
                                 gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAACTGCTGggtggggtgggg                                                                                                                      >  1:78913/1‑145 (MQ=255)
                                   cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCCGAAAACGCTGGGGGGGGTGGCg                                                                                                                      >  1:853977/1‑145 (MQ=255)
                                   cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGGTgg                                                                                                                        >  2:959727/1‑143 (MQ=255)
                                    tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTACCTCTGAAAAAGCTGGGTGGGGCGGTGAAAAAgg                                                                                                               >  2:759325/1‑151 (MQ=255)
                                     gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGc                                                                                                              <  1:107899/151‑1 (MQ=255)
                                      cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg                                                                                                             <  1:928504/151‑1 (MQ=255)
                                       gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGGGAAGAAGGCgg                                                                                                            >  1:902160/1‑151 (MQ=255)
                                            gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCgg                                                                                                            >  2:25300/1‑146 (MQ=255)
                                                gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                     <  1:1037970/149‑1 (MQ=255)
                                                gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                     <  2:263892/149‑1 (MQ=255)
                                                  acacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGGGAAGAAAGCGGGAGCcg                                                                                                      >  2:816726/1‑146 (MQ=255)
                                                     cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGGGGAGAAGGCGGGAGCcgc                                                                                                     >  1:355789/1‑144 (MQ=255)
                                                     cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTa                                                                                               <  2:1069806/150‑1 (MQ=255)
                                                      tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGTaa                                                                                              >  2:212932/1‑150 (MQ=255)
                                                       ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAt                                                                                             <  2:1177247/150‑1 (MQ=255)
                                                         aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCg                                                                                          <  1:485964/151‑1 (MQ=255)
                                                          gTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCCGGAATTCgg                                                                                         >  2:65857/1‑151 (MQ=255)
                                                          gTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCgg                                                                                         <  2:1005475/151‑1 (MQ=255)
                                                          gTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCgg                                                                                         <  1:659081/151‑1 (MQ=255)
                                                          gTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCg                                                                                          <  1:381213/150‑1 (MQ=255)
                                                           tGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAAAAAATGATGGATGGGGTTGCGGAGAAAGCGGAAGCCGCGCAGGaaat                                                                                            >  1:439284/1‑145 (MQ=255)
                                                             gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGGGGGAGCCGCGCAGGAATTCGGcc                                                                                       >  2:124756/1‑150 (MQ=255)
                                                             gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCt                                                                                      <  1:65852/151‑1 (MQ=255)
                                                             gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGGTGGGGGAGAAGGGGGGAGCCGCGGAGGAAATTGGcc                                                                                       >  1:10020/1‑150 (MQ=255)
                                                              cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGcc                                                                                       <  1:857988/149‑1 (MQ=255)
                                                               gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                    <  1:873361/151‑1 (MQ=255)
                                                                 tgtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTTGGCCtt                                                                                     >  2:1217082/1‑148 (MQ=255)
                                                                        aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATTCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAAATTGGccgtttcccgat                                                                             >  1:440804/1‑139 (MQ=255)
                                                                        aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGatat                                                                           <  1:279253/151‑1 (MQ=255)
                                                                          gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGGGGAGAAGGCGGGAGGCGCGCAGTAATTTGGCCTTAGCCAGAtata                                                                          >  1:186718/1‑150 (MQ=255)
                                                                          gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAg                                                                         <  2:795018/151‑1 (MQ=255)
                                                                             ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGt                                                                        <  2:1125609/149‑1 (MQ=255)
                                                                             ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                      <  2:1101798/151‑1 (MQ=255)
                                                                               cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGTCCTTAGCCAGATATAGTATGc                                                                    >  1:757776/1‑151 (MQ=255)
                                                                                    atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                               <  2:1027306/151‑1 (MQ=255)
                                                                                      ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGc                                                              <  1:622516/150‑1 (MQ=255)
                                                                                         ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCgg                                                           >  2:1123212/1‑150 (MQ=255)
                                                                                          gtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGAt                                                         <  2:10020/151‑1 (MQ=255)
                                                                                             aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                      <  1:664855/151‑1 (MQ=255)
                                                                                               aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                    <  2:1116968/151‑1 (MQ=255)
                                                                                               aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGGGGGGACCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGGAGGCGCGCAGGAATTCGGCCTTAGCCAGGTATAGTATGCGTAGGCCGGGTAAgg                                                     >  1:54085/1‑150 (MQ=255)
                                                                                                aCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCg                                                   <  2:349479/151‑1 (MQ=255)
                                                                                                aCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCg                                                   <  2:94601/151‑1 (MQ=255)
                                                                                                 ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                  <  1:884732/151‑1 (MQ=255)
                                                                                                  cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGtt                                                 <  1:1157874/151‑1 (MQ=255)
                                                                                                   gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGtt                                                 >  1:999579/1‑150 (MQ=255)
                                                                                                    cTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                >  1:1236408/1‑150 (MQ=255)
                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAACTGCTGGGTGGGGTGGGGAAGAAGGCGGGAGCCGCGCAGTAATTCGGGCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGTTTAc                                              >  1:427235/1‑150 (MQ=255)
                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCa                                               <  2:78919/149‑1 (MQ=255)
                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcg                                             <  1:285380/151‑1 (MQ=255)
                                                                                                        gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAAATGCTGGGTGGGCTGGCGAAGAAAGCGGGAGCCGCGCCGTAATTCCGCCTTAGCCAGATATAGTATTCGTAGGCCGGGTAAGGCGGTCAcgcc                                           >  1:859102/1‑151 (MQ=255)
                                                                                                        gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                           <  1:476005/151‑1 (MQ=255)
                                                                                                        gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                           <  1:441883/151‑1 (MQ=255)
                                                                                                         ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGGCGGCGGAGAAGGCGGGAGCCGCGCAGTAAATCGGGCTTAGCCAGGTATAGTATGCGTAGGGCGGGTAAGGCGGTCAc                                              >  1:1064072/1‑147 (MQ=255)
                                                                                                         ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGCTGGGGGAGAAGGCGGGAGCCGCGCAGGAATTTGGCCTTTGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgcc                                           >  1:1144435/1‑150 (MQ=255)
                                                                                                          cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTTGGCCTTAGCCCGGTATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccgc                                         >  1:240176/1‑151 (MQ=255)
                                                                                                          cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCCGTAAATCGGCCCTAGCCCGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcccc                                         >  1:1065576/1‑149 (MQ=255)
                                                                                                            ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAAGCGTTCACGCCCCa                                        >  1:648149/1‑150 (MQ=255)
                                                                                                            ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCa                                        <  1:595754/150‑1 (MQ=255)
                                                                                                               gACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGATCCATCAGAAAATGCTGGAAGGGGTGGCGAAGAAAGCGGGAGACGCGCCGTAAATTGGCCCTAACCAGATATAGGATGCGGAGGCCGGGGAAAGGGGTCACGCCGCATcc                                     >  1:126233/1‑150 (MQ=255)
                                                                                                               gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                     <  2:838920/150‑1 (MQ=255)
                                                                                                                aCCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                   <  1:25297/151‑1 (MQ=255)
                                                                                                                        cTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGGGAAGAAGGCGGAAGCCGCGCCGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCCttt                              >  1:699309/1‑148 (MQ=255)
                                                                                                                             gatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAg                      <  2:609604/151‑1 (MQ=255)
                                                                                                                                gaGCAACGGCGGTTGTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTg                    <  1:759258/150‑1 (MQ=255)
                                                                                                                                  gCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                 <  1:114246/151‑1 (MQ=255)
                                                                                                                                   cAACGGGGGTTTTGATGGGAGCCATCAGAAACGGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTACCCGGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGt                     <  1:495124/146‑1 (MQ=39)
                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc       >  1:233521/1‑151 (MQ=255)
                                                                                                                                              ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCtt     >  2:1175422/1‑151 (MQ=255)
                                                                                                                                                 gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt  >  1:1194593/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAT  >  NC_000913/4456955‑4457250

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