Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,396,233 C→T 6.3% R57C (CGC→TGC)  amiB → N‑acetylmuramoyl‑l‑alanine amidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,2330CT 6.3% 196.8 / 2.6 63R57C (CGC→TGC) amiBN‑acetylmuramoyl‑l‑alanine amidase II
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base T (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCTGCTGCTGCTGTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAAC  >  NC_000913/4396100‑4396379
                                                                                                                                     |                                                                                                                                                  
gctgctgctgctGTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTc                                                                                                                                      >  1:209728/1‑148 (MQ=255)
  tgctgctgctgTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGa                                                                                                                                    >  4:189987/1‑148 (MQ=255)
  tgctgctgctgTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGa                                                                                                                                    >  8:106538/1‑148 (MQ=255)
   gctgctgctgTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGata                                                                                                                                  >  4:115133/1‑149 (MQ=255)
       ctgctgTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATcaa                                                                                                                              <  3:189987/149‑1 (MQ=255)
       ctgctgTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATcaa                                                                                                                              <  3:115133/149‑1 (MQ=255)
               cACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACgg                                                                                                                      <  1:117500/149‑1 (MQ=255)
               cACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACgg                                                                                                                      <  6:223304/149‑1 (MQ=255)
                         ggtgCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCa                                                                                                            <  1:103205/149‑1 (MQ=255)
                               gcgACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATAAAACAAACGGGCGTGATTCAGGGACt                                                                                                      >  3:73695/1‑149 (MQ=255)
                               gcgACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATAACAAACAAACGGGCGTGATTCAGGGACt                                                                                                      >  1:218701/1‑149 (MQ=255)
                                    gCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCgt                                                                                                 <  2:209446/149‑1 (MQ=255)
                                        tctGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCgttgtt                                                                                             >  1:99322/1‑149 (MQ=255)
                                          tGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCgttgttgt                                                                                           >  4:86257/1‑149 (MQ=255)
                                          tGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCgttgttgt                                                                                           >  4:101962/1‑149 (MQ=255)
                                            atatTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCgttgttgtt                                                                                          <  1:35777/148‑1 (MQ=255)
                                               ttCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGc                                                                                      >  7:202787/1‑149 (MQ=255)
                                               ttCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGc                                                                                      <  5:1006/149‑1 (MQ=255)
                                               ttCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGc                                                                                      >  6:1006/1‑149 (MQ=255)
                                                tCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACgg                                                                                                                      <  7:119010/116‑1 (MQ=255)
                                                tCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACgg                                                                                                                      >  8:119010/1‑116 (MQ=255)
                                                     tttCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAATCAAATGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCaa                                                                                 >  4:5174/1‑148 (MQ=255)
                                                     tttCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATACAGGGACTGCCGTTGTTGTTCAGCGGCaa                                                                                 >  6:17009/1‑148 (MQ=255)
                                                      ttCTAAAGGTAATCAACATGCTCGGATTACGTTGAGTTTTATTGGCGATCCTGAATATCCGATTAGCCATCAAAGCAAACGAACAGTGGCGCTGGATATCAAACAAACGGCCGTGATTCAGGGACTGCCGTTCGTGTTAAGCGGCaat                                                                                >  4:35607/1‑148 (MQ=255)
                                                           aCGGTAATCAACATGCGCGGATAACGTTTATTTTTATTGTCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACAGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTg                                                                          <  2:155326/149‑1 (MQ=255)
                                                           aCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCTAACGCACCGTGGCGCTCGATATCAAACACACGGGCGTGATTCAGGGACGGCCGTTGTTGTTCAGCGGCaataat                                                                             <  2:226202/146‑1 (MQ=255)
                                                           aCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCt                                                                           <  5:219846/148‑1 (MQ=255)
                                                           aCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGAGTCAGGGACTGCCGTTGTTGTTCAGCGGCATTAATCTg                                                                          >  4:163859/1‑149 (MQ=255)
                                                           aCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGAGTTTAGCCATCAAAGCAAACGCACCGTTGCGCTCGATATCAAACAAACGGGAGTGATTCAGGGACTGCCGTTGTTGTTCATCGGCAATAATATg                                                                          >  5:4045/1‑149 (MQ=255)
                                                             ggTTATCAAGAGTCGCGGATAACTTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAATGCAAACGCACCGTGGCGCTCGATATCAAACTTACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGt                                                                        <  3:35607/149‑1 (MQ=255)
                                                                 aTCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGaa                                                                     <  5:47922/148‑1 (MQ=255)
                                                                 aTCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAg                                                                    >  2:121661/1‑149 (MQ=255)
                                                                      cAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGAt                                                               >  4:91052/1‑149 (MQ=255)
                                                                       aGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACTGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGAtt                                                              <  8:87079/149‑1 (MQ=255)
                                                                       aGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGAtt                                                              >  7:87079/1‑149 (MQ=255)
                                                                       aGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGAtt                                                              >  7:133523/1‑149 (MQ=255)
                                                                         gcgcGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCg                                                            <  3:101962/149‑1 (MQ=255)
                                                                            cgGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAgg                                                                                                          <  4:42716/100‑1 (MQ=255)
                                                                            cgGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAgg                                                                                                          >  3:42716/1‑100 (MQ=255)
                                                                                    gTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAACCGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGAGTCTGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGc                                                 <  5:17009/149‑1 (MQ=255)
                                                                                       gAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCt                                               <  1:121661/148‑1 (MQ=255)
                                                                                              aTTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATg                                        >  8:116349/1‑148 (MQ=255)
                                                                                              aTTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATg                                        >  6:178445/1‑148 (MQ=255)
                                                                                               ttGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGca                                      <  3:86257/149‑1 (MQ=255)
                                                                                                 ggCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGcac                                     >  4:75713/1‑148 (MQ=255)
                                                                                                        cTGATTATGCGTTTAGCCATCAAAGCAAATGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTa                             >  3:195506/1‑149 (MQ=255)
                                                                                                           aTTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACgg                          >  3:74534/1‑149 (MQ=255)
                                                                                                              aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACg                                                  >  6:19885/1‑122 (MQ=255)
                                                                                                              aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACg                                                  <  5:19885/122‑1 (MQ=255)
                                                                                                              aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGc                         <  6:159583/147‑1 (MQ=255)
                                                                                                              aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTGGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGc                         >  5:159583/1‑147 (MQ=255)
                                                                                                               tGCGTTAAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACTGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCtgg                      <  5:178557/149‑1 (MQ=255)
                                                                                                                        ccATCAAAGCAAATGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATc             <  3:5174/149‑1 (MQ=255)
                                                                                                                        ccATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATc             <  7:106538/149‑1 (MQ=255)
                                                                                                                        ccATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATc             >  1:99071/1‑149 (MQ=255)
                                                                                                                        ccATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATc             >  8:159504/1‑149 (MQ=255)
                                                                                                                          aTCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTg                                                                          <  6:135833/86‑1 (MQ=255)
                                                                                                                          aTCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTg                                                                          >  5:135833/1‑86 (MQ=255)
                                                                                                                          aTCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCCAAAGATGCACAAACGCTACGGCTGGTGGTCGAt              <  2:231531/146‑1 (MQ=255)
                                                                                                                          aTCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCCAAAGATGCACAAACGCTACGGCTGGTGGTCGAt              >  1:231531/1‑146 (MQ=255)
                                                                                                                             aaaGCAAATGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTAcc        <  4:195506/149‑1 (MQ=255)
                                                                                                                             aaaGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTAc         <  4:133400/148‑1 (MQ=255)
                                                                                                                                 cccccGCCCCGTGGCGCTCGATACCCACCAAACGGGCGTGATTCAGTGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGGGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGaaa    <  1:206139/145‑1 (MQ=255)
                                                                                                                                   aaCGCACCGGGGCGCTCGATATCAAACAAACGGGCGTGATTTAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAAc  <  6:186420/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                  
GCTGCTGCTGCTGTGCACGCCGGTGGGTGCCGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAAC  >  NC_000913/4396100‑4396379

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑