Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,413,542 A→T Y67N (TAT→AAT)  ydaC ← DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,5420AT98.4% 203.0 / NA 62Y67N (TAT→AAT) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (27/33);  minor base C (0/1);  total (27/35)
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

AAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTC  >  NC_000913/1413399‑1413687
                                                                                                                                               |                                                                                                                                                 
aaGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATg                                                                                                                                              <  2:423552/149‑1 (MQ=255)
aaGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATg                                                                                                                                              <  1:278149/149‑1 (MQ=255)
aaGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCAACGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATg                                                                                                                                              <  2:423478/149‑1 (MQ=255)
     tCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAAtt                                                                                                                                         <  1:692140/149‑1 (MQ=255)
        ggCCGTATTCAAGCATGCGCGCTTCATTAAAGATTACGTGTGCCATCGCAGCACCTGTGGCAGGCGCCAATGTTAGACTGAAATTGACTCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTTAATGCAATtctc                                                                                                                                      <  2:403644/149‑1 (MQ=255)
        ggCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATtctc                                                                                                                                      <  2:315060/149‑1 (MQ=255)
         gCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCg                                                                                                                                     >  1:758570/1‑149 (MQ=255)
          ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                    <  1:254024/149‑1 (MQ=255)
          ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                    >  2:152805/1‑149 (MQ=255)
             tATGCAACCTTCCCCTCTCCATTAAAGATTCCGTGAGCCACCGCATCGCCTCGGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGa                                                                                                                                 <  1:22720/149‑1 (MQ=255)
                       tCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtc                                                                                                                       >  1:449943/1‑149 (MQ=255)
                         ccctctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACCGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                     <  2:731118/147‑1 (MQ=255)
                         cACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                     <  1:766002/149‑1 (MQ=255)
                         cACGCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                     <  2:650975/149‑1 (MQ=255)
                           ctctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGc                                                                                                                   >  1:815482/1‑149 (MQ=255)
                            tctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCa                                                                                                                  >  1:369332/1‑149 (MQ=255)
                             cgtCATTCAAGATTACTTGTGCCCTCGCCTCACCTCTCGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCaa                                                                                                                 <  2:319450/147‑1 (MQ=255)
                                 aTTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaa                                                                                                             >  1:640021/1‑149 (MQ=255)
                                  ttAAAGATTACTTCTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGCTTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTATCGAGGACTGAAGTTTCTCGCAATTaaa                                                                                                            >  2:567892/1‑149 (MQ=255)
                                      aGATTACTTGTGCCATCGCATCACCGCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAAttt                                                                                                        <  2:272348/149‑1 (MQ=255)
                                        aTTACTTGTGCCATCGCATCAACTCTGGCCGGCGCCAATGTTAGACTGAAATTGGAGCCCGATTTGCATTATTAATAAGCAGCTATTAAGTTGTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTAt                                                                                                      <  2:477963/149‑1 (MQ=255)
                                            cTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGt                                                                                                  >  1:691067/1‑149 (MQ=255)
                                               gtgCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCCGtttt                                                                                               >  1:307843/1‑149 (MQ=255)
                                               gtgCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATATATCAGtttt                                                                                               >  1:104916/1‑149 (MQ=255)
                                                ggCCATCGCGTCAGCGCGGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTa                                                                                              <  1:738818/148‑1 (MQ=255)
                                                  ccATCGCATCACCTCTGGCAGGCGCCAATGGTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACt                                                                                            <  2:92933/149‑1 (MQ=255)
                                                     tCGCAGCACCTCTGGCAGGCGCCACTGTTAGACTGAAATCGACGCCCGATGTGGATTATTAATAATCAGCTATGAAGTTTTCATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTCATCAGTTTTACTTTc                                                                                         <  2:125100/149‑1 (MQ=255)
                                                         aTCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTTTgat                                                                                     >  2:166580/1‑147 (MQ=255)
                                                             ccTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGGAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctct                                                                                 <  2:582242/149‑1 (MQ=255)
                                                             ccTCTGGCAGGCGCCAATGTTAGACTGAAATCGACGCCAGATGTAGCATATTAAAAAACAGCTATGAAGTTTTAATTTGAATTCAAGGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctct                                                                                 <  1:229839/149‑1 (MQ=255)
                                                               tctGGAAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTgg                                                                               <  2:428509/149‑1 (MQ=255)
                                                                ctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGa                                                                              <  1:82487/149‑1 (MQ=255)
                                                                   gccggCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCCCTCGAAAc                                                                           >  2:22434/4‑149 (MQ=255)
                                                                   gCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAAc                                                                           >  2:18221/1‑149 (MQ=255)
                                                                         gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGAAAACGCCTGc                                                                     >  2:17629/1‑149 (MQ=255)
                                                                         gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGATTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                     >  2:69108/1‑149 (MQ=255)
                                                                           cAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGAACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGctt                                                                   >  2:255401/1‑149 (MQ=255)
                                                                               gTTAGACTTAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGCGGCCTGAAGCTTCTCGCAATTAAAATTTATCAGTTTTACTTCCTGCTCTCTGGAAACGCCTGCTTCttt                                                               <  2:480904/149‑1 (MQ=255)
                                                                                  aGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTCCtttttt                                                            <  1:685460/149‑1 (MQ=255)
                                                                                  aGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTACGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGACACGCCTGCTGCtttgtt                                                            <  2:424568/149‑1 (MQ=255)
                                                                                   gACGGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTCAGAGGTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTa                                                           <  1:229574/149‑1 (MQ=255)
                                                                                      tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCt                                                        >  1:266462/1‑149 (MQ=255)
                                                                                           ttGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGagag                                                   >  1:31002/1‑149 (MQ=255)
                                                                                           ttGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTCagag                                                   >  1:161814/1‑149 (MQ=255)
                                                                                               cccccAATGTTGATTATTAATAATCAGCTACGAAGTTTTAATTTGAATTCAATGAAATTCTAGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACCCCTGCTTCTTTTTTACCTGAGAGCAta                                               <  1:166580/147‑2 (MQ=255)
                                                                                                gCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAttt                                              >  2:146217/1‑149 (MQ=255)
                                                                                                   cGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAtttttt                                           >  1:119735/1‑149 (MQ=255)
                                                                                                          gATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGCGAGCCTTTTTTCTCATTc                                    >  2:155657/1‑149 (MQ=255)
                                                                                                           attattAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCt                                   <  2:691067/149‑1 (MQ=255)
                                                                                                               ttaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGAtt                               >  1:784500/1‑149 (MQ=255)
                                                                                                                 aataatCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAAAGCCTGCTTCGTTTTTTCCTGAGAGCATTTTTTCGCATTCTGTTTTc                             >  2:85821/1‑149 (MQ=255)
                                                                                                                  ataatCAGCTATGAAGTTTTAATTTGAATCCAATGTAATTCTCGAGGACCGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGACACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCg                            <  1:22434/149‑1 (MQ=255)
                                                                                                                   taatCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGt                           <  2:640021/149‑1 (MQ=255)
                                                                                                                     atCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTa                         <  1:446523/149‑1 (MQ=255)
                                                                                                                      tCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAg                        >  1:370704/1‑149 (MQ=255)
                                                                                                                            aTGAAGTTTTCATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGa                  <  1:714346/149‑1 (MQ=255)
                                                                                                                                atttttAATTTGAGTTTTATGTAATTCTGTAGGATGGAAGTTTCGCTCAATTAAAACTTATCTGTTCTACTTTCTTTTCTCTTGAATCGCCTGCTTCTTTTTTACCTTAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAtttt              <  2:229490/147‑1 (MQ=255)
                                                                                                                                 gTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTg             <  1:18221/149‑1 (MQ=255)
                                                                                                                                      aaTTTGAATGAAATGCAATGCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAatat        <  2:689618/149‑1 (MQ=255)
                                                                                                                                         ttGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCtt     >  2:624398/1‑149 (MQ=255)
                                                                                                                                           gAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGt   >  1:698929/1‑149 (MQ=255)
                                                                                                                                            aaTTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTc  <  2:491552/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
AAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTC  >  NC_000913/1413399‑1413687

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