Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,168 Δ1 bp 5.7% coding (2227/2334 nt) bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,1680A. 5.7% 779.4 / 37.2 229coding (2227/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base A (103/108);  new base . (6/7);  total (112/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-02

Only 100 of 231 total aligned reads displayed.

GCGCTGATGCAGTATGTGCAGCAGCGCCATCCCCACCTGATGCTGGAGGTTTATCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAAC  >  NC_000913/3713913‑3714434
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                          
gcgcTGATGCAGTATGTGCAGCAGCGCCATCCCCACCTGATGCTGGAGGTTTATCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCCCGCGCTGGGGAGATCTTTGGTCAgaa                                                                                                                                                                                                                                                     >  1:269847/1‑277 (MQ=255)
 cgcTGATGCAGTATGTGCAGCAGCGCCATCCCCACCTGATGCTGGAGGTTTATCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                                                                                                                                         <  2:521614/274‑1 (MQ=255)
                     cagcGCCATCCCCACCTGATGCTGGAGGTTTATCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                                                                                                                                                     <  2:602658/274‑1 (MQ=255)
                                                  ttATCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCCAGCTGGGGAGATCTTTGGTTAGCACGTTCCCTGCGCTGTTTTTACAACTACCGTTCGCGGTTAAATCCa                                                                                                                                                                                                        >  2:479176/1‑274 (MQ=255)
                                                     tCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCCAGCTGGGGAGATCTTTGGTCAGCCCGTCCCCGGCCCCGTTTTTCCAAATCCCGCCAGCGGTTAAACCCatct                                                                                                                                                                                                     >  4:381871/1‑271 (MQ=255)
                                                                   cGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCAGGTGGTTAAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCGGCCCCTTATTACCCGCACACCCATTACACAGCCGCGAGCTGGGGAGATCTTTGGTCCGCCCGTTCCCTTCACCTTTTTTAAAATTACCCACCGCGGTTAAACCCAAGTCCGGCCATGCCCCt                                                                                                                                                                                       >  4:561165/1‑274 (MQ=255)
                                                                           aaTTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                                                                                                                                                              <  2:381875/275‑1 (MQ=255)
                                                                                                                gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCg                                                                                                                                         >  2:218365/1‑275 (MQ=255)
                                                                                                                 cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTggg                                                                                                                                                                                                                                                                      <  5:475182/149‑1 (MQ=255)
                                                                                                                 cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTggg                                                                                                                                                                                                                                                                      <  7:1393742/149‑1 (MQ=255)
                                                                                                                  gcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgggg                                                                                                                                                                                                                                                                     <  5:123454/149‑1 (MQ=255)
                                                                                                                   cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                                                                                                                                                    <  5:1005209/149‑1 (MQ=255)
                                                                                                                    aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGgag                                                                                                                                                                                                                                                                   >  5:1573776/1‑149 (MQ=255)
                                                                                                                    aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGGAAGTTTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGgag                                                                                                                                                                                                                                                                   <  8:1500047/149‑1 (MQ=255)
                                                                                                                     tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaat                                                                                                                                                                                                                                                                  <  5:818989/149‑3 (MQ=255)
                                                                                                                     tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaat                                                                                                                                                                                                                                                                  <  7:601340/149‑3 (MQ=255)
                                                                                                                          ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                                                                                                                                                             <  7:667442/149‑1 (MQ=255)
                                                                                                                           gATGAAACCCAACTACCCACATGGCTTATGAGCTGGCCGCTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGAGGAGATCTTTg                                                                                                                                                                                                                                                            <  7:1110666/149‑1 (MQ=255)
                                                                                                                           gATGAAACCCAACTACCCACATGGATTATGTGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTg                                                                                                                                                                                                                                                            <  5:257598/149‑1 (MQ=255)
                                                                                                                           gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTg                                                                                                                                                                                                                                                            <  7:1056798/149‑1 (MQ=255)
                                                                                                                           gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTg                                                                                                                                                                                                                                                            <  5:182378/149‑1 (MQ=255)
                                                                                                                            aTGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTgg                                                                                                                                                                                                                                                           >  8:1351823/1‑149 (MQ=255)
                                                                                                                             tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                                                                                                                                          >  8:1650814/1‑149 (MQ=255)
                                                                                                                             tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                                                                                                                                          >  5:1637934/1‑149 (MQ=255)
                                                                                                                             tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                                                                                                                                          <  5:565266/149‑1 (MQ=255)
                                                                                                                             tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                                                                                                                                          >  5:1336364/1‑149 (MQ=255)
                                                                                                                               aaaCCCAACTGCCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGCCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCTCAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCa                                                                                                                                                                                                                                                        >  8:1096258/1‑149 (MQ=255)
                                                                                                                                 aCCCAACGACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                                                                                                                                      >  6:932686/1‑149 (MQ=255)
                                                                                                                                    cAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                                                                                                                                                   >  6:1356100/1‑149 (MQ=255)
                                                                                                                                      aCTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAAAATTTGCCCgcgtgc                                                                                                                >  1:113296/1‑278 (MQ=255)
                                                                                                                                         aCCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCAc                                                                                                                                                                                                                                              <  7:761079/149‑1 (MQ=255)
                                                                                                                                          cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                                                                                                                                      <  5:1119450/140‑1 (MQ=255)
                                                                                                                                          cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACt                                                                                                                                                                                                                                             <  6:474461/149‑1 (MQ=255)
                                                                                                                                          cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGCAGATCTTTGGTCAGc                                                                                                                                                                                                                                                      >  6:1119450/1‑140 (MQ=255)
                                                                                                                                           ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCtt                                                                                                                                                                                                                                                              >  8:1203868/1‑131 (MQ=255)
                                                                                                                                           ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTg                                                                                                                                                                                                                                            >  6:275143/1‑149 (MQ=255)
                                                                                                                                             acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                                                                                                                                                          >  8:1517719/1‑149 (MQ=255)
                                                                                                                                              caTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                                                                                                                                                         >  8:147899/1‑149 (MQ=255)
                                                                                                                                                    ttATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                                                                                                                                                   >  7:1308938/1‑149 (MQ=255)
                                                                                                                                                          gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                                                                                                                                                             <  7:1552578/149‑1 (MQ=255)
                                                                                                                                                            tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATa                                                                                                                                                                                                                           >  8:811637/1‑149 (MQ=255)
                                                                                                                                                            tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCAGCCCAGGCCAGCGCGGTATTACCGCACCAGCAATCCCCCAGCCGCGAGAGGGGGATATCTTTGGTCACCACGTTCACGGCCCCGTTTTCACaattc                                                                                                                                                                                                                           >  6:1545536/1‑146 (MQ=255)
                                                                                                                                                             ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                                                                                                                                                    >  3:258750/1‑107 (MQ=255)
                                                                                                                                                             ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAc                                                                                                                                                                                                                          >  7:1672325/1‑149 (MQ=255)
                                                                                                                                                             ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTACTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCCTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCCTGAATGCAAATCACGCCAGTGTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGGCTTCCGTCCACAAACGAAccaccg                                                                                            >  4:350756/1‑270 (MQ=255)
                                                                                                                                                              gCCGGTGGTTCAAACGCTGTAAGTGTCCGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGCGCTGGGGAGATCTTTGGTCAGCACGTTCACTTCACCGTTTTTCCAAATCCCGTCAGCGGTTAAATCCATATCCGGCCATGCCCCTTCGTGCATTCAAATCa                                                                                                                                                                       >  4:153348/1‑199 (MQ=255)
                                                                                                                                                                cGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGt                                                                                                                                                                                                                       <  7:228434/149‑1 (MQ=255)
                                                                                                                                                                 ggtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTACTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCca                                                                                                          <  2:350760/257‑1 (MQ=255)
                                                                                                                                                                       tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                                                                                                                                                                                                      >  8:969933/1‑127 (MQ=255)
                                                                                                                                                                        cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGtt                                                                                                                                                                                                               >  5:1620793/1‑149 (MQ=255)
                                                                                                                                                                               tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                                                                                                                                        >  7:486806/1‑149 (MQ=255)
                                                                                                                                                                                 tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                                                                                                                                          <  5:1384997/97‑1 (MQ=255)
                                                                                                                                                                                 tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                                                                                                                                                      <  6:1330383/149‑1 (MQ=255)
                                                                                                                                                                                  aaGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                                                                                                                                                     <  7:947752/149‑1 (MQ=255)
                                                                                                                                                                                  aaGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                                                                                                                                                     >  5:1359765/1‑149 (MQ=255)
                                                                                                                                                                                    gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                                                                                                                                                     <  6:643373/147‑1 (MQ=255)
                                                                                                                                                                                    gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATcc                                                                                                                                                                                                   >  7:62082/1‑149 (MQ=255)
                                                                                                                                                                                      gtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCgg                                                                                                                                                                                                 <  5:1329982/149‑1 (MQ=255)
                                                                                                                                                                                       tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                                                                                                                                          >  8:1008183/1‑91 (MQ=255)
                                                                                                                                                                                         aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                                                                                                                                                            <  8:712242/135‑1 (MQ=255)
                                                                                                                                                                                         aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                                                                                                                                              <  5:468293/149‑1 (MQ=255)
                                                                                                                                                                                         aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCCCTCCTGGGc                                                                >  2:160071/1‑275 (MQ=255)
                                                                                                                                                                                          gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                                                                                                                                                             >  7:530994/1‑149 (MQ=255)
                                                                                                                                                                                           ttCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                                                                                                                                                            >  7:163200/1‑149 (MQ=255)
                                                                                                                                                                                                  ccGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc                                                                                                                                                                                     <  6:335312/149‑1 (MQ=255)
                                                                                                                                                                                                   cGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAAt                                                                                                                                                                         >  4:64205/1‑160 (MQ=255)
                                                                                                                                                                                                    gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCCTTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGCTTTTACAAATACCGTCAGCGGCTAAATCCAGATCCGGCCCTGCCCCTTCGt                                                                                                                                                                                   >  8:946100/1‑149 (MQ=255)
                                                                                                                                                                                                    gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                                                                                                                                                   <  8:1308938/149‑1 (MQ=255)
                                                                                                                                                                                                     ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTAGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                                                                                                                                                              >  1:77103/1‑153 (MQ=255)
                                                                                                                                                                                                       gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACg                                                                                                                                                                     <  2:64205/160‑1 (MQ=255)
                                                                                                                                                                                                         aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                                                                                                                                                              >  6:1172934/1‑149 (MQ=255)
                                                                                                                                                                                                          tttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                                                                                                                                                   >  5:679208/1‑143 (MQ=255)
                                                                                                                                                                                                          tttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                                                                                                                                                   <  6:679208/143‑1 (MQ=255)
                                                                                                                                                                                                              tCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTCATGACCGCTCCGGCAAGAATTTCCCCGCGTTCCTTCCACAACCGAACAGTACCCCCATCTGTTATTCCGCCCTCCTGGGCACCGTCAGGATGAATGGtta                                            >  4:471564/1‑272 (MQ=255)
                                                                                                                                                                                                                   ccATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGc                                                                                                                                                                    <  7:924784/149‑1 (MQ=255)
                                                                                                                                                                                                                     aTGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTTGGGAGATCTTTGGTCAGCACGTTCCCTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCAGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCa                                    <  2:293735/275‑1 (MQ=255)
                                                                                                                                                                                                                        ccAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                                                                                                                                               <  7:877628/149‑1 (MQ=255)
                                                                                                                                                                                                                              ccGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                                                                                                                                               >  8:446492/1‑143 (MQ=255)
                                                                                                                                                                                                                              ccGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTATTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAtt                                                                                                                                                         <  7:466427/149‑1 (MQ=255)
                                                                                                                                                                                                                                         gCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                                                                                                                                              >  6:779770/1‑149 (MQ=255)
                                                                                                                                                                                                                                         gCACAGCCATTCCCCAGCAGTGAGCTGGGGAGATCTTTGGTCTGTACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                                                                                                                                              <  6:1328710/149‑1 (MQ=255)
                                                                                                                                                                                                                                           acaGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCa                                                                                                                                                                       >  5:925554/1‑122 (MQ=255)
                                                                                                                                                                                                                                                ccATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                                                                                                                                              <  6:1046241/95‑1 (MQ=255)
                                                                                                                                                                                                                                                ccATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                                                                                                                                              >  5:1046241/1‑95 (MQ=255)
                                                                                                                                                                                                                                                 cATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACt                                                                                                                                                                                                                                             <  6:1421732/47‑1 (MQ=255)
                                                                                                                                                                                                                                                 cATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCt                                                                                                                                       >  8:488888/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 atctCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                                                                                                                                                 >  6:296077/4‑75 (MQ=37)
                                                                                                                                                                                                                                                 atctCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                                                                                                                                                 <  5:296077/72‑1 (MQ=37)
                                                                                                                                                                                                                                                  tctCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCt                                                                                                                                                                                       >  5:644439/3‑100 (MQ=255)
                                                                                                                                                                                                                                                    tCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg                                                                                                                                    >  6:911308/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ccccAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCgg                                                                                                                                   >  7:1513784/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ccccAGCCGCGAGGTGGGGAGATCTTTGTTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGACAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCagag     <  2:568722/275‑1 (MQ=255)
                                                                                                                                                                                                                                                     ccccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCgg                                                                                                                                   >  6:1489030/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      cccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGc                                                                                                                                  <  7:1142977/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGc                                                                                                                                  <  8:225119/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa                                                                                                                                 >  8:1021957/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gaGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCCCTGCGCCGTTCTTACAAGTACCGTGAGAGGTTAAATGCAGATCCGGCCGTGCCGCTTCGGGAATGCCAGTCACGCCAGGTTTAATTCCCTCGCTCATGACCGCTCCGGCaa                                                                                                                                >  6:979064/2‑149 (MQ=255)
                                                                                                                                                                                                                                                             cgcgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAAt                                                                                                                            <  5:1190767/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cgcgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAAt                                                                                                                            <  8:935455/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               cgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGTATTTGCCCgcgtgcgt                                                                                                              >  4:100137/1‑161 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                          
GCGCTGATGCAGTATGTGCAGCAGCGCCATCCCCACCTGATGCTGGAGGTTTATCAAAAAAATCAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAAC  >  NC_000913/3713913‑3714434

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