Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 641,366 T→C 5.3% intergenic (+48/+73) ahpF → / ← uspG alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding/universal stress protein UP12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913641,3660TC 5.3% 483.2 / 10.8 150intergenic (+48/+73)ahpF/uspGalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding/universal stress protein UP12
Reads supporting (aligned to +/- strand):  ref base T (73/69);  new base C (5/3);  total (78/72)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 150 total aligned reads displayed.

TACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTC  >  NC_000913/641219‑641511
                                                                                                                                                   |                                                                                                                                                 
tACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                                  >  4:63468/1‑149 (MQ=255)
tACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                                  >  7:349060/1‑149 (MQ=255)
            gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                                      >  5:216675/1‑149 (MQ=255)
            gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                                      <  8:158689/149‑1 (MQ=255)
            gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                                      <  1:278579/149‑1 (MQ=255)
            gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                                      >  6:253395/1‑149 (MQ=255)
             tACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTAt                                                                                                                                     >  3:282103/1‑149 (MQ=255)
                   cAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                                <  1:286906/132‑1 (MQ=255)
                      atcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaa                                                                                                                             >  2:310424/1‑148 (MQ=255)
                       tcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaaaa                                                                                                                           <  6:315090/149‑1 (MQ=255)
                            atcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGccc                                                                                                                      >  7:51923/1‑149 (MQ=255)
                            atcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGccc                                                                                                                      >  1:374362/1‑149 (MQ=255)
                             tcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGcccc                                                                                                                     >  2:251762/1‑149 (MQ=255)
                                ccACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgcc                                                                                                                  <  6:404161/149‑1 (MQ=255)
                                 cACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGACGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgccg                                                                                                                 <  1:169965/149‑1 (MQ=255)
                                 cACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgccg                                                                                                                 >  7:179188/1‑149 (MQ=255)
                                    tGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGa                                                                                                              >  2:279396/1‑149 (MQ=255)
                                      gCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAt                                                                                                             <  3:70510/148‑1 (MQ=255)
                                        gAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGc                                                                                                          >  7:35073/1‑149 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                                >  6:265856/1‑102 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                                <  5:265856/102‑1 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCa                                                                                                                                 <  7:341104/117‑1 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCa                                                                                                                                 >  8:341104/1‑117 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                                 >  4:146160/1‑149 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                                 >  1:106578/1‑149 (MQ=255)
                                                 aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                                 <  5:314055/149‑1 (MQ=255)
                                                 aaaGCCGCTCTGAGTGCCTTTGACTACCTGATTCCCACCAAAACTGCATAAGAAGAACTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGGCAAACCTGTCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCGCCGGTCAg                                                                                                 <  6:367143/149‑1 (MQ=255)
                                                   aGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAggg                                                                                               <  7:420640/149‑1 (MQ=255)
                                                   aGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAggg                                                                                               >  1:302883/1‑149 (MQ=255)
                                                     ccTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCg                                                                                             >  3:336491/1‑149 (MQ=255)
                                                      ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCgg                                                                                            <  8:179188/149‑1 (MQ=255)
                                                        ctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGa                                                                                          >  4:147259/1‑149 (MQ=255)
                                                             aGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAc                                                                                       <  5:161440/147‑1 (MQ=255)
                                                               tGCTTTTGACTACCTGATTCGCACCAAAAATGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTAt                                                                                   >  4:316069/1‑149 (MQ=255)
                                                                  ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGa                                                                                          >  5:448216/1‑139 (MQ=255)
                                                                  ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGa                                                                                          <  6:448216/139‑1 (MQ=255)
                                                                  ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCt                                                                                >  7:135437/1‑149 (MQ=255)
                                                                  ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCt                                                                                >  2:390800/1‑149 (MQ=255)
                                                                    ttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTg                                                                               <  5:253395/148‑1 (MQ=255)
                                                                      gACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGc                                                                             <  4:282103/148‑1 (MQ=255)
                                                                         tACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTgg                                                                              <  6:201390/144‑1 (MQ=255)
                                                                           ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                                       <  4:336491/149‑1 (MQ=255)
                                                                           ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                                       >  4:220416/1‑149 (MQ=255)
                                                                           ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                                       >  1:447631/1‑149 (MQ=255)
                                                                            cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                                      >  8:383924/1‑149 (MQ=255)
                                                                            cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                                      <  1:279396/149‑1 (MQ=255)
                                                                              gATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAg                                                                                                                         <  1:19621/96‑1 (MQ=255)
                                                                                    cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAATGCACAAc                                                               >  4:244067/1‑148 (MQ=255)
                                                                                    cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAcc                                                              <  3:428952/149‑1 (MQ=255)
                                                                                    cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGAAATTAACGCACAAc                                                               >  4:90683/1‑148 (MQ=255)
                                                                                     aCCAAAACTGCGTAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAGAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                               >  3:332042/1‑147 (MQ=255)
                                                                                     aCCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                               >  4:316729/1‑147 (MQ=255)
                                                                                         aaaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAc                                                         <  8:51923/149‑1 (MQ=255)
                                                                                         aaaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAc                                                         >  1:401255/1‑149 (MQ=255)
                                                                                             tGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCgg                                                      >  2:7904/1‑148 (MQ=255)
                                                                                             tGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCgg                                                      >  2:448860/1‑148 (MQ=255)
                                                                                               cATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAg                                                   <  3:316069/149‑1 (MQ=255)
                                                                                                aTAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCgg                                                      >  1:442016/1‑145 (MQ=255)
                                                                                                 tAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGAt                                                 >  6:341418/1‑149 (MQ=255)
                                                                                                     aagaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggc                                             <  3:98714/149‑1 (MQ=255)
                                                                                                        aagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtg                                          <  7:414754/149‑1 (MQ=255)
                                                                                                         agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtg                                          <  2:442016/148‑1 (MQ=255)
                                                                                                          gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                        <  1:7904/149‑1 (MQ=255)
                                                                                                          gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                        >  3:198258/1‑149 (MQ=255)
                                                                                                           tAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACaa                                                                >  2:283934/1‑124 (MQ=255)
                                                                                                           tAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                       >  2:464369/1‑149 (MQ=255)
                                                                                                              gATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAc                                                                                       >  3:356176/1‑98 (MQ=255)
                                                                                                              gATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAc                                                                                       <  4:356176/98‑1 (MQ=255)
                                                                                                               aTTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGAt                                    >  6:68069/1‑148 (MQ=255)
                                                                                                               aTTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                                   <  3:220416/149‑1 (MQ=255)
                                                                                                               aTTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCACCACCGGCAGATTGGCGTGGCGGATc                                   <  5:341418/149‑1 (MQ=255)
                                                                                                                ttCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCa                                  >  6:200055/1‑149 (MQ=255)
                                                                                                                ttCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCa                                  <  2:302883/149‑1 (MQ=255)
                                                                                                                   aCCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                               >  7:170757/1‑149 (MQ=255)
                                                                                                                   aCCTGCAATAGATTAGCCGCTGGGGTCAAACATGCCGGCTTTTTTATCGCAATAAAAAGCCCCGCCGGGATGGCTCCCGGCAGGGCGGAATAATTAATCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGTCGGATAACGc                               >  8:313007/1‑149 (MQ=255)
                                                                                                                    ccTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCg                                                       >  4:226365/1‑124 (MQ=255)
                                                                                                                     cTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTc                             <  8:349060/149‑1 (MQ=255)
                                                                                                                        cAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAg                          >  8:411835/1‑149 (MQ=255)
                                                                                                                          ataaaaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGc                        >  1:24804/8‑149 (MQ=255)
                                                                                                                          aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                                   <  7:50095/138‑1 (MQ=255)
                                                                                                                          aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                                   >  8:50095/1‑138 (MQ=255)
                                                                                                                          aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGc                        <  1:71406/149‑1 (MQ=255)
                                                                                                                           ttGCTTAGCCGCCGGTGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                       >  4:78896/1‑149 (MQ=255)
                                                                                                                           ttGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                       >  4:51264/1‑149 (MQ=255)
                                                                                                                              cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt                      >  7:393734/1‑147 (MQ=255)
                                                                                                                              cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa                    <  5:68069/149‑1 (MQ=255)
                                                                                                                              aaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt                      >  2:127454/4‑147 (MQ=255)
                                                                                                                              aaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa                    <  2:24804/146‑1 (MQ=255)
                                                                                                                               ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCTCCGGGATGGCTCCGGGCAGGGAGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATAGGCGTGGCGGATCACGCTCGAGGCGt                      >  6:232991/1‑146 (MQ=255)
                                                                                                                               ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt                      >  1:223627/1‑146 (MQ=255)
                                                                                                                               ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAAGGCTCCGGGCAGGGCGGGATACTTAATCTGGCAATTAACGCAAAACCATCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt                      >  8:421354/1‑146 (MQ=255)
                                                                                                                                   ccgccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAcc               >  8:439337/1‑149 (MQ=255)
                                                                                                                                     gccgGGGTCAGACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                               >  5:334297/1‑131 (MQ=255)
                                                                                                                                     gccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTa             <  2:106578/149‑1 (MQ=255)
                                                                                                                                      ccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCTGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa            <  3:146160/149‑1 (MQ=255)
                                                                                                                                      ccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa            <  3:147259/149‑1 (MQ=255)
                                                                                                                                            tCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                                   >  8:153358/1‑120 (MQ=255)
                                                                                                                                            tCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATg      <  1:251762/149‑1 (MQ=255)
                                                                                                                                              aaaCCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAcc               <  5:15517/138‑1 (MQ=255)
                                                                                                                                              aaaCCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCGCGAGGCGTTAGAACCTAACAGATgg     <  2:86935/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTC  >  NC_000913/641219‑641511

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