Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,210 A→G 19.7% K287E (AAA→GAA)  amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,2100AG19.7% 155.6 / 37.4 86K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (11/6);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

GTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGC  >  NC_000913/2553113‑2553308
                                                                                                 |                                                                                                  
gTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTACACTGGTTCGACAACCAGAAAGCACATTCGaaa                                                                                                  >  2:547802/1‑100 (MQ=255)
 ttAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGaaaa                                                                                                 <  2:802960/100‑1 (MQ=255)
  tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGaaaaa                                                                                                >  1:1239682/1‑100 (MQ=255)
    ggCACGGCTGAGTTTTGTCATAAACTCGCGCCATAGATTGCTGATGGCGTGATCAGTTATTTACTCTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                              <  2:1913045/100‑1 (MQ=255)
    ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCTATTGCTGAAGGCGTGATCAGTTATTTCCTCTGGTTTGACATCTAGAAAGCAAATTAGAAAAAGc                                                                                              <  2:242366/100‑1 (MQ=255)
    ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAg                                                                                               <  1:455655/99‑1 (MQ=255)
    ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                              >  1:473349/1‑100 (MQ=255)
    ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                              >  1:1096263/1‑100 (MQ=255)
     gCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCg                                                                                             >  1:2166915/1‑100 (MQ=255)
      catggcggcGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGCCAACCAGAAAGCACATTCGAAAAAGc                                                                                              <  1:1363555/95‑1 (MQ=255)
         ggcggcGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATaa                                                                                         >  1:872306/1‑100 (MQ=255)
           cggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                                                       <  1:180180/100‑1 (MQ=255)
             gcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCTACAACCAGAAATCACATTCGGAAAAGCGATAAGTTa                                                                                     <  2:1753038/100‑1 (MQ=255)
                   cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAACGCGATAAGTTATGAAAc                                                                               >  1:1941556/1‑100 (MQ=255)
                     tCATAAAATCGCCACATCGATTGCGGAAGGCGTGATCAGTTATTTCCACTTTTTCGACAACCAGAAAGCACATTCGATAAAGCGATAAGTTATGAAAcc                                                                              <  2:1926117/99‑1 (MQ=255)
                     tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                             >  2:705622/1‑100 (MQ=255)
                      cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcac                                                                             >  1:184990/1‑97 (MQ=255)
                        gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGAc                                                                          >  1:1665226/1‑100 (MQ=255)
                        gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGAc                                                                          <  2:935901/100‑1 (MQ=255)
                         aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                                                         <  2:480737/100‑1 (MQ=255)
                          aaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAACGCGATAAGTTATGAAACCCGACGc                                                                        >  2:488547/1‑100 (MQ=255)
                             tCGCCACAGCGATTTCTGAAGCCTTGCTCAGTTAGTGCCACTCGTTCGACCACCATAAAGCACATCCGAAAAAGCGCTAAGTTATCAAACCCGAcaca                                                                       <  2:937610/98‑1 (MQ=255)
                              cGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac                                                                    <  1:1704642/100‑1 (MQ=255)
                               gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                   >  1:576802/1‑100 (MQ=255)
                               gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATCAAACCCGACGCACAcc                                                                   >  2:846359/1‑100 (MQ=255)
                                ccACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCa                                                                  >  2:1492492/1‑100 (MQ=255)
                                      cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                            <  2:2166900/100‑1 (MQ=255)
                                      cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                            <  1:1334142/100‑1 (MQ=255)
                                      cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                            <  1:1764826/100‑1 (MQ=255)
                                         ttGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACa                                                         <  1:1545366/100‑1 (MQ=255)
                                         ttGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACa                                                         <  1:1729497/100‑1 (MQ=255)
                                           gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAg                                                        >  1:1247488/1‑99 (MQ=255)
                                              gAAGGCGTGGTCAGTTATTTCCACTGGTTCGTCAACCAGAAATCACGTTCGAAAAGGCGGTAAGGTAGGAAACCCGACGCACAGGAGGTTAAGCAGTTTc                                                    >  1:2134552/1‑100 (MQ=255)
                                              gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACAAGGTTAAACAGTTTc                                                    >  2:509911/1‑100 (MQ=255)
                                              gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                    <  1:1045539/100‑1 (MQ=255)
                                               aaGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCt                                                   <  1:1091714/100‑1 (MQ=255)
                                                 ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                                 >  2:42796/1‑100 (MQ=255)
                                                  gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                                <  1:813956/100‑1 (MQ=255)
                                                   cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                               >  2:111845/1‑100 (MQ=255)
                                                   cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                               >  1:1377551/1‑100 (MQ=255)
                                                      gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                            >  1:469296/1‑100 (MQ=255)
                                                      gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                               >  2:1401946/1‑97 (MQ=255)
                                                        tCAGTGATTTCCACTGGTTCGACAACCAGTAAGCATATTCGAAAAAGCGATAAGTGATGAAACACGACGCACTCCAGGTTAAACAGTTTCTGc                                                 <  2:1727463/93‑1 (MQ=255)
                                                           gTTATTTCCACTGGTTCAACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATAAAACACGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                       >  2:911814/1‑100 (MQ=255)
                                                             tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                     <  1:1785352/100‑1 (MQ=255)
                                                              aTTTCCACTGGTTCGACAACTAGAAAGCACATTCGAAAAAGCGATAAGTTAGGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                    <  1:439100/100‑1 (MQ=255)
                                                              aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                    <  1:1002157/100‑1 (MQ=255)
                                                              aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                    >  2:1785340/1‑100 (MQ=255)
                                                              aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                    <  2:1604430/100‑1 (MQ=255)
                                                              aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTCCTCAACCTTCAg                                      >  1:2141636/1‑98 (MQ=255)
                                                                ttCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                                  >  2:688944/1‑100 (MQ=255)
                                                                ttCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                                  >  2:388459/1‑100 (MQ=255)
                                                                ttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                                  <  2:1257536/100‑1 (MQ=255)
                                                                   cACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGa                               >  1:971608/1‑100 (MQ=255)
                                                                   cACTGGTTCGACAACCAAAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGTATACGa                               >  2:454535/1‑100 (MQ=255)
                                                                     cTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAtt                             >  2:959941/1‑100 (MQ=255)
                                                                       ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                           >  1:292202/1‑100 (MQ=255)
                                                                        gTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                          <  1:1777472/100‑1 (MQ=255)
                                                                         ttCGACAACCTGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAATAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                         <  2:1474001/100‑1 (MQ=255)
                                                                         ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                         <  2:1914556/100‑1 (MQ=255)
                                                                         ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                         <  1:1063903/100‑1 (MQ=255)
                                                                         ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                         <  1:804670/100‑1 (MQ=255)
                                                                          tCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTca                        <  1:1500609/100‑1 (MQ=255)
                                                                           cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGCTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                       <  2:1223387/100‑1 (MQ=255)
                                                                           cGACAACCAGAAAGCACATTCGAAAAATCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                         <  2:2163193/98‑1 (MQ=255)
                                                                            gACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                      >  1:977583/1‑100 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCGGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    <  2:962304/100‑1 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    >  1:778636/1‑100 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    >  2:29195/1‑100 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    >  1:986375/1‑100 (MQ=255)
                                                                               aaCCAGAAAGCACATTCGAAATAGCGATAAGTTATGAAAACCGACGCACACCAGGTTAAACAGTTTCTGCTCAACTTTCAGGATACTATTTGTcagc                      <  2:1546415/97‑1 (MQ=255)
                                                                                aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                  >  1:356689/1‑100 (MQ=255)
                                                                                aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                  >  1:2061693/1‑100 (MQ=255)
                                                                                aCCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                  <  2:1496702/100‑1 (MQ=255)
                                                                                 ccAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                  >  2:711070/1‑99 (MQ=255)
                                                                                 ccAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTg                 >  1:2184095/1‑100 (MQ=255)
                                                                                 ccAGAAAGCACAATCGAAAAAGCGTTAAGTTATGAAACCCGATGCACACCAGGTTAAACTGTTTCTGCTCAACCTTCAGTATACTATTTGTCAgctgct                  <  2:1721785/99‑1 (MQ=255)
                                                                                  caaaaaGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGa                <  2:2141621/97‑1 (MQ=255)
                                                                                    gAAAGCACATTCGAAAAAGCGATAAGTTATGGAACCCCACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAcc              <  2:1516261/100‑1 (MQ=255)
                                                                                     aaaGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccg             <  2:1049249/100‑1 (MQ=255)
                                                                                        gCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgccg          >  1:631316/1‑100 (MQ=255)
                                                                                        gCACATTCGAAAAAGCGAAAAGATATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCACCAGCTGAccgc            >  2:1908332/1‑98 (MQ=255)
                                                                                         cacaTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTCCTCAACCTTCAGGATACGATTTGTCAGCTGCTGAcaaatcg         <  1:1565674/100‑7 (MQ=255)
                                                                                           caTTCGAAAAAGCGATAAGTTATGACACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc        <  2:1007025/99‑1 (MQ=255)
                                                                                            aTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGa      >  2:46336/1‑100 (MQ=255)
                                                                                              tCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATg    <  2:1946592/100‑1 (MQ=255)
                                                                                               cGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAACTGACCGCCGTCGATgg   >  2:981846/1‑100 (MQ=255)
                                                                                                gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGgc  >  2:167577/1‑100 (MQ=255)
                                                                                                gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGgc  >  2:1415460/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGC  >  NC_000913/2553113‑2553308

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