Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 414,832 G→T 5.4% H308N (CAT→AAT)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913414,8320GT 5.4% 386.8 / 7.4 110H308N (CAT→AAT) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (54/50);  new base T (4/2);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 117 total aligned reads displayed.

GGTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCT  >  NC_000913/414685‑414980
                                                                                                                                                   |                                                                                                                                                    
ggTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATgc                                                                                                                                                     >  3:11274/1‑149 (MQ=255)
ggTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATgc                                                                                                                                                     >  8:193434/1‑149 (MQ=255)
ggTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGGATTTCATTCTTCAATCTGCTGGCGAATATgc                                                                                                                                                     >  5:210477/1‑149 (MQ=255)
      gtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAg                                                                                                                                                >  1:126359/1‑148 (MQ=255)
          cTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgcagat                                                                                                                                           <  2:126359/149‑5 (MQ=255)
           tGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGcgccgc                                                                                                                                          >  1:178361/1‑149 (MQ=255)
             tAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCGCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGcgccgc                                                                                                                                          <  5:28852/147‑1 (MQ=255)
               aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                      >  3:175852/1‑149 (MQ=255)
               aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                      >  2:400017/1‑149 (MQ=255)
               aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCt                                                                                                                                       <  1:230183/148‑1 (MQ=255)
               aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCt                                                                                                                                       <  7:193434/148‑1 (MQ=255)
                ccAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgt                                                                                                                                     >  7:184760/1‑149 (MQ=255)
                ccAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTACCGTGGATTTAATTCGTCAATAATCTGGCGAATTTGCGCCAGGGCCGCGCtggg                                                                                                                                    >  1:90013/1‑147 (MQ=255)
                 cAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgtg                                                                                                                                    <  4:274802/149‑1 (MQ=255)
                 cAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgtg                                                                                                                                    <  8:189978/149‑1 (MQ=255)
                    gTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCGGGGTGGTGGCGAATGCTCGCGGGTTGCGGCATTGTGCTCTGTTAGCGAGTATTTTCTTCTTCATTCTGCTGGCGAATATGCGCCAGCGCGGCGCTGTGTTc                                                                                                                                 >  3:197926/1‑149 (MQ=255)
                       ttCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTg                                                                                                                               <  4:175852/148‑1 (MQ=255)
                          aGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCTAATATGCGCCAGCGCCGCGCTGTGTTCTGCGAt                                                                                                                           <  7:165220/149‑1 (MQ=255)
                          aGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGAt                                                                                                                           >  6:400096/1‑149 (MQ=255)
                          aGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGAt                                                                                                                           >  5:117527/1‑149 (MQ=255)
                           ggCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATg                                                                                                                          <  5:400096/149‑1 (MQ=255)
                           ggCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATg                                                                                                                          <  1:400017/149‑1 (MQ=255)
                            gCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGc                                                                                                                         <  8:184760/149‑1 (MQ=255)
                                tgttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGt                                                                                                                       <  7:323426/147‑1 (MQ=255)
                                    gctgctgctGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCa                                                                                                                  <  8:215615/148‑1 (MQ=255)
                                    gctgctgctGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCa                                                                                                                  <  6:160358/148‑1 (MQ=255)
                                      tgctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTc                                                                                                                                 <  7:168731/131‑1 (MQ=255)
                                      tgctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTc                                                                                                                                 >  8:168731/1‑131 (MQ=255)
                                        ctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTAGTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtgt                                                                                                              <  6:100861/148‑1 (MQ=255)
                                          gctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCt                                                                                                                                       <  4:152895/121‑1 (MQ=255)
                                          gctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCt                                                                                                                                       >  3:152895/1‑121 (MQ=255)
                                            tgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGAt                                                                                                                           <  6:84490/131‑1 (MQ=255)
                                            tgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGAt                                                                                                                           >  5:84490/1‑131 (MQ=255)
                                            tgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtgt                                                                                                              <  1:100434/144‑1 (MQ=255)
                                            tgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtgt                                                                                                              >  2:100434/1‑144 (MQ=255)
                                            tgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGAc                                                                                                          >  4:270229/1‑148 (MQ=255)
                                                aTTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACg                                                                                                         >  5:282411/1‑145 (MQ=255)
                                                aTTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGa                                                                                                     <  2:178361/149‑1 (MQ=255)
                                                 ttCTGCTGACTGCTTCTCCGCGTGGTGGCTAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGGGGACGAAGAt                                                                                                    <  1:319756/149‑1 (MQ=255)
                                                 ttCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGa                                                                                                           <  5:231626/142‑1 (MQ=255)
                                                      ctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTc                                                                                                 <  6:117527/147‑1 (MQ=255)
                                                      ctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGt                                                                                               >  1:297873/1‑149 (MQ=255)
                                                       tgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                              >  4:168241/1‑149 (MQ=255)
                                                        gACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATgcgc                                                                                                                                                   <  2:338446/95‑1 (MQ=255)
                                                           tGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgt                                                                                                                <  6:369581/127‑1 (MQ=255)
                                                           tGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCgg                                                                                          >  5:869/1‑149 (MQ=255)
                                                             cTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGtt                                                                                        >  5:1414/1‑149 (MQ=255)
                                                              ttCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCa                                                                                                                  >  2:107588/1‑122 (MQ=255)
                                                              ttCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCa                                                                                                                  <  1:107588/122‑1 (MQ=255)
                                                              ttCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTtg                                                                                       >  5:401052/1‑149 (MQ=255)
                                                                            gCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGa                                                                                                           >  7:61102/1‑115 (MQ=255)
                                                                            gCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGa                                                                                                           <  8:61102/115‑1 (MQ=255)
                                                                            gCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGAc                                                                                                          <  4:269088/116‑1 (MQ=255)
                                                                            gCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGAc                                                                                                          >  3:269088/1‑116 (MQ=255)
                                                                              gAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAgcgc                                                                        <  8:349613/148‑1 (MQ=255)
                                                                              gAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAgcgc                                                                        >  7:349613/1‑148 (MQ=255)
                                                                              gAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGcgcc                                                                       >  2:178139/1‑149 (MQ=255)
                                                                                   cTCGCGCGAAGCGACATTGTGCTCTCTAAGCCAGAATATACTTCTTCAATCTCCTGGCCAATATGCGCCAGCGCCGCTCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAg                                                                  <  4:197926/149‑1 (MQ=255)
                                                                                    tCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAgcgcaacctc                                                                  >  2:254584/1‑142 (MQ=255)
                                                                                    tCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGt                                                                 >  7:171683/1‑149 (MQ=255)
                                                                                      gcgcgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGtt                                                                >  3:70371/1‑148 (MQ=255)
                                                                                      gcgcgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGtt                                                                <  1:98279/148‑1 (MQ=255)
                                                                                             gcgcCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTgc                                                        >  3:178238/1‑149 (MQ=255)
                                                                                                  aTTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttt                                                   <  4:70371/149‑1 (MQ=255)
                                                                                                   ttGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCtttt                                                  <  8:271527/149‑1 (MQ=255)
                                                                                                      tgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttc                                                <  1:390918/148‑1 (MQ=255)
                                                                                                      tgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttc                                                >  2:390918/1‑148 (MQ=255)
                                                                                                        ctctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttttt                                                 >  3:26879/1‑145 (MQ=255)
                                                                                                        ctctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttttt                                                 <  4:26879/145‑1 (MQ=255)
                                                                                                         tctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttct                                            >  7:226601/1‑149 (MQ=255)
                                                                                                          ctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttctt                                           >  4:186176/1‑149 (MQ=255)
                                                                                                          ctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttctt                                           >  4:267452/1‑149 (MQ=255)
                                                                                                           tGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGAGTTATGAGATGCGTTCCCAGTGAGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGATGAGGTTGCGCTTTttcttcttc                                          >  5:305601/1‑149 (MQ=255)
                                                                                                             tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttcttc                                          <  3:318357/147‑1 (MQ=255)
                                                                                                             tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttcttc                                          <  3:168241/147‑1 (MQ=255)
                                                                                                             tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttctt                                           <  3:186176/146‑1 (MQ=255)
                                                                                                             tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCg                                         <  6:1414/148‑1 (MQ=255)
                                                                                                              taGCGAGTATTTACTTCTTCAATCTGCTGGCCAATATGCGCCAGCGCCGCGCTGTTTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                        <  2:90013/147‑1 (MQ=255)
                                                                                                               aGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGc                                       <  6:401052/148‑1 (MQ=255)
                                                                                                                  gAGTTTTTACTTCTGCTATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTACCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAAc                                   <  6:305601/149‑1 (MQ=255)
                                                                                                                  gAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGtt                                                                                                                      >  5:285171/1‑66 (MQ=255)
                                                                                                                  gAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGtt                                                                                                                      <  6:285171/66‑1 (MQ=255)
                                                                                                                  gAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAAc                                   >  5:158813/1‑149 (MQ=255)
                                                                                                                      aTTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCt                               >  1:259831/1‑149 (MQ=255)
                                                                                                                       tttACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCt                               <  2:297873/148‑1 (MQ=255)
                                                                                                                          aCTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGtt                            <  1:178139/148‑1 (MQ=255)
                                                                                                                           cttcttCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGc                          <  6:255973/149‑1 (MQ=255)
                                                                                                                              cttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCaa                        <  3:267452/148‑1 (MQ=255)
                                                                                                                              cttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAg                       >  1:124674/1‑149 (MQ=255)
                                                                                                                                   aTCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctg                  <  3:24336/149‑1 (MQ=255)
                                                                                                                                         tgGCGCATATTCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcg            >  2:173086/1‑149 (MQ=255)
                                                                                                                                         tgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcg            >  2:4885/1‑149 (MQ=255)
                                                                                                                                           gcgcATATTCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGtt                            <  8:203189/131‑1 (MQ=255)
                                                                                                                                             gAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                        >  1:281621/1‑117 (MQ=255)
                                                                                                                                             gAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcggcgg        >  1:251115/1‑149 (MQ=255)
                                                                                                                                             gAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcggcgg        >  7:220543/1‑149 (MQ=255)
                                                                                                                                             gAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcggcgg        >  8:392480/1‑149 (MQ=255)
                                                                                                                                                 atGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctgg    >  2:201240/1‑149 (MQ=255)
                                                                                                                                                  tGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctggc   <  7:81472/149‑1 (MQ=255)
                                                                                                                                                   gcgcCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctggct  >  8:40149/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCT  >  NC_000913/414685‑414980

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