Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 232,959 C→T A333T (GCC→ACC)  mltD ← putative membrane‑bound lytic murein transglycosylase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913232,9590CT98.8% 317.8 / NA 86A333T (GCC→ACC) mltDputative membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (47/37);  minor base A (1/0);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTGGCCTGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGCAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAG  >  NC_000913/232813‑233103
                                                                                                                                                  |                                                                                                                                                
aCTTTGGCCTGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCTCTGCTGCAAATGTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAAcc                                                                                                                                              <  2:971730/151‑1 (MQ=255)
aCTTTGGCCTGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAAcc                                                                                                                                              <  1:855157/151‑1 (MQ=255)
aCTTTGGCCTGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAAcc                                                                                                                                              <  1:745375/151‑1 (MQ=255)
aCTTTGGCCTGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAAc                                                                                                                                               >  1:58840/1‑150 (MQ=255)
         tGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGAc                                                                                                                                     <  1:212360/151‑1 (MQ=255)
            cTTCATCTTAGATCCGCGCTGTTTGTTCCACTGCTGCAAATCTTTGGTTCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCGAGAGCGTACGGTGTAAACACGGCTGATAAGCGGCGGATTGTCGGTAACCAGCGTCGACttt                                                                                                                                  >  1:56504/1‑149 (MQ=255)
            cTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGa                                                                                                                                      >  2:613843/1‑147 (MQ=255)
            cTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGt                                                                                                                                  <  2:606836/151‑1 (MQ=255)
              tCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAc                                                                                                                                <  1:775071/151‑1 (MQ=255)
                 gCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcag                                                                                                                              >  1:75146/1‑150 (MQ=255)
                 gCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcag                                                                                                                              >  2:705475/1‑150 (MQ=255)
                 gCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcag                                                                                                                              >  1:660341/1‑150 (MQ=255)
                  cTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcagc                                                                                                                             <  1:124633/150‑1 (MQ=255)
                  cTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcag                                                                                                                              <  2:909228/149‑1 (MQ=255)
                  cTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcag                                                                                                                              <  1:1013832/149‑1 (MQ=255)
                   ttAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTAcagcag                                                                                                                           <  2:85146/151‑1 (MQ=255)
                      gATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCaa                                                                                                                        >  1:908338/1‑151 (MQ=255)
                      gATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGNNNNNNNGGTAACCANNGTCGACTGTACNGCAGCaa                                                                                                                        >  1:232081/1‑151 (MQ=255)
                              cAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGcc                                                                                                                 <  1:868697/150‑1 (MQ=255)
                                gTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGa                                                                                                              >  2:926438/1‑151 (MQ=255)
                                gTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGa                                                                                                              >  1:264558/1‑151 (MQ=255)
                                gTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGa                                                                                                              >  2:249600/1‑151 (MQ=255)
                                  ttGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAg                                                                                                            >  1:989879/1‑151 (MQ=255)
                                     ttCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGcc                                                                                                          <  1:329658/150‑1 (MQ=255)
                                             ctgcAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCa                                                                                                 >  2:316966/1‑151 (MQ=255)
                                              tgcAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCa                                                                                                 >  1:348314/1‑150 (MQ=255)
                                                cAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGc                                                                                              <  2:568651/151‑1 (MQ=255)
                                                    tCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGtt                                                                                          >  1:380054/1‑151 (MQ=255)
                                                        tGGTGCTTACGCCGGGACGTGAAGCGATACTTGAAAGCGGGTCGCCAGAGCGGACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGGAACCAGCGTCGACTGTGCAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATc                                                                                      <  2:957289/151‑1 (MQ=255)
                                                           tGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTg                                                                                    <  2:58349/150‑1 (MQ=255)
                                                           tGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTg                                                                                    <  2:622487/150‑1 (MQ=255)
                                                           tGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTg                                                                                    <  2:442158/150‑1 (MQ=255)
                                                             cTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCAt                                                                                 <  1:597922/151‑1 (MQ=255)
                                                                aCGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATg                                                                                >  2:799165/1‑149 (MQ=255)
                                                                aCGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGct                                                                              <  2:616404/151‑1 (MQ=255)
                                                                   ccGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGcttc                                                                            <  1:510368/150‑1 (MQ=255)
                                                                             aaGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCAt                                                                 >  2:290681/1‑151 (MQ=255)
                                                                                cGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATc                                                                <  2:748120/149‑1 (MQ=255)
                                                                                cGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCa                                                               >  1:85825/1‑150 (MQ=255)
                                                                                cGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCa                                                               <  2:396915/150‑1 (MQ=255)
                                                                                  aTACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGt                                                            >  2:529824/1‑151 (MQ=255)
                                                                                  aTACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGt                                                            <  1:1008242/151‑1 (MQ=255)
                                                                                   tACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGt                                                            <  2:517301/150‑1 (MQ=255)
                                                                                    aCTTGAAGGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTAc                                                          <  1:655268/151‑1 (MQ=255)
                                                                                    aCTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTa                                                           >  1:244561/1‑150 (MQ=255)
                                                                                    aCTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTAc                                                          <  1:280104/151‑1 (MQ=255)
                                                                                    aCTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTAc                                                          <  1:288556/151‑1 (MQ=255)
                                                                                      ttGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTg                                                        <  1:232622/151‑1 (MQ=255)
                                                                                            gCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGc                                                   >  1:923690/1‑150 (MQ=255)
                                                                                             cGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGc                                                   <  2:621072/149‑1 (MQ=255)
                                                                                                 tCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTg                                             <  1:316871/151‑1 (MQ=255)
                                                                                                 tCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTg                                             >  2:563218/1‑151 (MQ=255)
                                                                                                          cGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGc                                     <  2:359720/150‑1 (MQ=255)
                                                                                                           gTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGt                                   >  1:435973/1‑151 (MQ=255)
                                                                                                               ggTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:484960/1‑151 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:1052408/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  2:711595/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:286321/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  2:613221/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  2:908584/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:334426/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  2:486791/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:394938/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  2:411168/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:432501/1‑150 (MQ=255)
                                                                                                                gtgtAAACACGGCTGTTAAGCGGCGTATTGACGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCCGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTggcg                               >  1:255970/1‑148 (MQ=255)
                                                                                                                gtgtAAACACAGCTGTTAAGCAGCGCATTGGCGGAAACCAGCGTCGACTGTACACGAGCAACTTCGCCTGAAGCCAGAGATCCACGCCATTGATCTGCAATCTTCTTTGGCACCATCACGTACCGCGGGCCACTTGCGCCCAGCGTGGAg                               >  1:115439/1‑150 (MQ=255)
                                                                                                                         cGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGc                     <  1:332392/151‑1 (MQ=255)
                                                                                                                         cGGCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGc                     <  2:384794/151‑1 (MQ=255)
                                                                                                                           gCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGcc                    <  1:153249/150‑1 (MQ=255)
                                                                                                                           gCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGcc                    <  1:553147/150‑1 (MQ=255)
                                                                                                                           gCTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCa                   >  2:562452/1‑151 (MQ=255)
                                                                                                                            cTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAg                  >  2:829924/1‑151 (MQ=255)
                                                                                                                            cTGTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAg                  >  2:68078/1‑151 (MQ=255)
                                                                                                                              gTTAAGCGGCGTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCg                >  2:192839/1‑151 (MQ=255)
                                                                                                                                    cggcgTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAAt          <  1:1057944/151‑1 (MQ=255)
                                                                                                                                        gTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGt        >  1:321333/1‑149 (MQ=255)
                                                                                                                                        gTATTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGGGGAGCCTTTTACGCCAGCGTTGAATGt        >  2:17025/1‑149 (MQ=255)
                                                                                                                                          aTTGTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTc    <  2:565691/151‑1 (MQ=255)
                                                                                                                                             gTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTAc                       >  2:76799/1‑129 (MQ=255)
                                                                                                                                             gTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTAc                       <  1:76791/129‑1 (MQ=255)
                                                                                                                                             gTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAg  >  2:717057/1‑150 (MQ=255)
                                                                                                                                             gTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAg  >  2:589193/1‑150 (MQ=255)
                                                                                                                                             gTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAg  >  1:810985/1‑150 (MQ=255)
                                                                                                                                             gTCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGCGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAg  >  2:710347/1‑150 (MQ=255)
                                                                                                                                              tCGGTAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAg  >  2:535481/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACTTTGGCCTGGCTTCAGCTTAGATCCGCGCAGTTTGTTCCACTGCTGCAAATCTTTGGTGCTTACGCCGAGACGTGAAGCGATACTTGAAAGCGTGTCGCCAGAGCGTACGGTGTAAACACGGCTGTTAAGCGGCGTATTGTCGGCAACCAGCGTCGACTGTACAGCAGCAATTTCGCCTGAAGCCAGAGATTCACGCAGTTGATCTGCATGCTTCTTTGGCACCATCACGTACTGCGGGCCACTTGCGCCCAGCGTGGAGCCTTTTACGCCAGCGTTGAATGTCTTCAG  >  NC_000913/232813‑233103

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