Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,688 A→T 16.3% I287N (ATT→AAT)  dacD ← D‑alanyl‑D‑alanine carboxypeptidase; penicillin‑binding protein 6b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,081,6880AT16.3% 205.0 / 32.8 86I287N (ATT→AAT) dacDD‑alanyl‑D‑alanine carboxypeptidase; penicillin‑binding protein 6b
Reads supporting (aligned to +/- strand):  ref base A (31/41);  new base T (6/8);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

GGTAGCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATG  >  NC_000913/2081589‑2081786
                                                                                                   |                                                                                                  
ggTAGCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCaaaa                                                                                                    <  2:1280391/100‑1 (MQ=255)
ggTAGCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCAACGGTGCaaaa                                                                                                    >  1:1787911/1‑100 (MQ=255)
   aGCACCATCCAGATCTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAAttt                                                                                                 <  2:2168920/100‑1 (MQ=255)
   aGCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAtttt                                                                                                 <  2:972825/100‑1 (MQ=255)
   aGCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAAttt                                                                                                 <  1:30560/100‑1 (MQ=255)
    gCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTg                                                                                                <  1:1531628/100‑1 (MQ=255)
     cACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGc                                                                                               <  2:146624/100‑1 (MQ=255)
      aCCATCCAGAACTCTTGTTCCGTTCCCAGGTCTATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCa                                                                                              <  1:747559/100‑1 (MQ=255)
      aCCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCa                                                                                              <  1:1113294/100‑1 (MQ=255)
      aCCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACATTTTTCCCACGGTGCAAAATTTGCa                                                                                              <  2:1806447/100‑1 (MQ=255)
         aTCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCg                                                                                           >  2:619048/1‑100 (MQ=255)
          tCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgt                                                                                          <  2:109061/100‑1 (MQ=255)
          tCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCg                                                                                           <  2:1954230/99‑1 (MQ=255)
          tCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTACACCgt                                                                                          >  1:759032/1‑100 (MQ=255)
           ccAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtg                                                                                         <  2:1045554/100‑2 (MQ=255)
           ccAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtg                                                                                         <  1:2260454/100‑2 (MQ=255)
           ccAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtg                                                                                         <  1:1441310/100‑2 (MQ=255)
           ccAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCgta                                                                                         <  2:354677/100‑1 (MQ=255)
            cAGAACTCTTGTTCTGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtgg                                                                                        <  2:99808/100‑3 (MQ=255)
            cAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtgg                                                                                        <  2:430841/100‑3 (MQ=255)
            cAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtgg                                                                                        <  1:2412356/100‑3 (MQ=255)
            cAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAATTTTGCACCgtgg                                                                                        <  2:2230639/100‑3 (MQ=255)
            cAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCgtag                                                                                        <  2:2341857/100‑1 (MQ=255)
            cAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCgtag                                                                                        <  1:619045/100‑1 (MQ=255)
                 ctctTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATCCCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAg                                                                                   >  1:353286/1‑100 (MQ=255)
                 ctctTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAg                                                                                   <  1:1526678/100‑1 (MQ=255)
                    ttGTTCCGTTCCTAGGTCGTTATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGttt                                                                                <  1:1692196/100‑1 (MQ=255)
                       ttccgtaccCAGTTCGATATTTTCTTTTTCGCCAGACCAGTTGCGTTACGTTCCGACCTTCTTCCCACAGTTCAAAATTTGCTCCGTAGTAAAGTTttg                                                                              <  1:1210769/99‑1 (MQ=255)
                         ccgttccCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTttgtt                                                                            >  1:560077/1‑99 (MQ=255)
                           gttccCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGcc                                                                         <  1:1486206/100‑1 (MQ=255)
                              cccAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGccccc                                                                      >  2:1391434/1‑100 (MQ=255)
                              cccAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGccccc                                                                      >  2:630586/1‑100 (MQ=255)
                              cccAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGccccc                                                                      <  2:832335/100‑1 (MQ=255)
                               ccAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCa                                                                     >  2:1951719/1‑100 (MQ=255)
                                  ggTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCCCGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACg                                                                  >  1:2330803/1‑100 (MQ=255)
                                    tCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCa                                                                >  2:314055/1‑100 (MQ=255)
                                    tCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCa                                                                >  2:1589206/1‑100 (MQ=255)
                                     cGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAg                                                               >  2:1799610/1‑100 (MQ=255)
                                      gATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGt                                                              <  1:76290/100‑1 (MQ=255)
                                       atatTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTa                                                             >  1:388201/1‑100 (MQ=255)
                                        tgttttCTTTTTCGCCATACCAGGTGCCTTCCCTTCCTACCTTTTTCCCACGGTTCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGt                                                              <  1:2172059/96‑1 (MQ=255)
                                        tatTTTCTTTATCGCCATACCAGATGCGTTACGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCAACGTAGTATAGTTTTGTTGCCCCCAACGCAGt                                                              >  1:416476/1‑98 (MQ=255)
                                          ttttCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAAtt                                                          >  2:1220551/1‑100 (MQ=255)
                                           tttCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAAttt                                                         >  2:140807/1‑100 (MQ=255)
                                           tttCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAAttt                                                         <  2:1073601/100‑1 (MQ=255)
                                             tCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTc                                                       <  1:1955785/100‑1 (MQ=255)
                                                  aTCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGcc                                                  <  2:1536405/100‑1 (MQ=255)
                                                   tCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCct                                                 >  2:161344/1‑100 (MQ=255)
                                                   tCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCct                                                 >  1:40498/1‑100 (MQ=255)
                                                     gCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCctct                                               >  1:469645/1‑100 (MQ=255)
                                                      ccATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCtt                                              <  2:2563653/100‑1 (MQ=255)
                                                        aTACCAGATGCGTTCCGTTCCGACCTTTTGCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTcc                                            <  1:227679/100‑1 (MQ=255)
                                                         tACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCt                                           >  1:107171/1‑100 (MQ=255)
                                                         tACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCt                                           >  1:1360580/1‑100 (MQ=255)
                                                            cggATGCGTTCCGGTCCGACCTGTGTCCCATGGTGCAATATGTGCGCGGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGGCTCTTCCtacc                                        <  1:1135981/98‑4 (MQ=255)
                                                            cAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCAc                                        <  2:2047686/100‑1 (MQ=255)
                                                            cAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCAc                                        >  2:2536142/1‑100 (MQ=255)
                                                            cAGATGCGTTCCGTTCCGACCTGTTTACCACCGTGGTAAATTTGCAACGTAGTAAAGTTT‑‑TTGCCCCCTACGCAGTAATTTTCTTGCCTCTTCCTCTCGa                                      <  1:370859/100‑1 (MQ=255)
                                                             aGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACg                                       <  2:469647/100‑1 (MQ=255)
                                                              gATGCGTTCCGTTCCGACCTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGAc                                     >  2:1960051/1‑100 (MQ=255)
                                                               aTGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGAc                                     <  2:198124/100‑1 (MQ=255)
                                                                     tccgttccgACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTg                               <  1:744941/100‑1 (MQ=255)
                                                                        gttccgACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCAc                            <  2:1249265/100‑1 (MQ=255)
                                                                         ttccgACCTTTTTCCCTCAGTGCAAAATTTGCACCGTAGTACAGTTTTGTTGCCTCCAACTCATTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCAc                            <  1:1165356/99‑1 (MQ=255)
                                                                         ttccgACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACt                           <  2:1709351/100‑1 (MQ=255)
                                                                          tccgACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTg                          >  1:833338/1‑100 (MQ=255)
                                                                            cgACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTc                        >  1:1945965/1‑100 (MQ=255)
                                                                            cgACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTCGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTc                        <  1:599131/100‑1 (MQ=255)
                                                                              aCCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAg                      <  2:525143/100‑1 (MQ=255)
                                                                               aaactttactaCGGTGCAAATTTTGCACCGTGGGAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGc                     >  1:2490346/11‑100 (MQ=25)
                                                                                ccttttCCCACGGAGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcc                    >  2:1757566/3‑100 (MQ=255)
                                                                                cTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTTTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcc                    <  1:297919/100‑1 (MQ=255)
                                                                                cTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcc                    >  2:1901027/1‑100 (MQ=255)
                                                                                aactttactaCGGTGCAAATTTTGCACCGTGGGAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCTTCTTCCTCACGACCTTTTGCACTGTCTGcc                    >  1:147825/10‑100 (MQ=25)
                                                                                 tttttCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGccc                   <  2:1984517/100‑1 (MQ=255)
                                                                                  ttttCCCACGGTGCAAAATTTGAACCGTAGTAAAGTTTTGTTGGCCGCAACGCCGTAAGTTTGTTGCTTCTTCCTCGCGACCTTTTGCACTGTCAGcccc                  >  1:612912/1‑100 (MQ=255)
                                                                                   tttCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGccccc                 <  1:1901017/100‑1 (MQ=255)
                                                                                   tttCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGccccc                 <  1:1757557/100‑1 (MQ=255)
                                                                                      cccACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAtt              <  2:1978706/100‑1 (MQ=255)
                                                                                      actaCGGTGCAAATTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCATCCCCCAtt              >  1:611686/4‑100 (MQ=255)
                                                                                       ctaCGGTGCAAATTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTa             >  2:1843167/3‑100 (MQ=255)
                                                                                       ctaCGGTGCAAATTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTATCAGCCCCCATTa             >  1:1472317/3‑100 (MQ=255)
                                                                                       ccACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTa             >  1:1074264/1‑100 (MQ=255)
                                                                                         aCGGTGCAAATTTTGCAACGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTTCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTAcc           >  2:122935/1‑100 (MQ=255)
                                                                                         aCGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTAcc           <  1:641090/100‑1 (MQ=255)
                                                                                           ggTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCAc         <  2:107172/100‑1 (MQ=255)
                                                                                             tGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTCATTTTCTTGCCTCTTCCTCCCGCCCTTTTGCACTGTCCGCCCCCATTccc           >  1:1158900/1‑96 (MQ=255)
                                                                                             tGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTg       >  1:773784/1‑100 (MQ=255)
                                                                                              gCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGc      >  2:907589/1‑100 (MQ=255)
                                                                                                aaaaTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCaa    >  2:2345194/1‑100 (MQ=255)
                                                                                                aaaaTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCaa    >  1:938040/1‑100 (MQ=255)
                                                                                                 aaaTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAAt   >  2:1363297/1‑100 (MQ=255)
                                                                                                  aaTTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATg  <  1:2071602/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GGTAGCACCATCCAGAACTCTTGTTCCGTTCCCAGGTCGATATTTTCTTTATCGCCATACCAGATGCGTTCCGTTCCGACCTTTTTCCCACGGTGCAAAATTTGCACCGTAGTAAAGTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATG  >  NC_000913/2081589‑2081786

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