Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,044 A→G Q50Q (CAA→CAG serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,0440AG96.1% 152.3 / ‑2.5 52Q50Q (CAA→CAGserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (25/24);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAA  >  NC_000913/4624810‑4625281
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
gaggcTTAATGGTGCCAGAAGGAATCGATACAAGGTAGCCTCATGCGTTATTATCCCTGCTTCGAACGATTTTTCAGGAGCCTTAATGCCTAACATTTCCTGGTGCGACCTGCCTGAGGATGTCTCTTTATGGCCGGGGCTGCCTCTTTCAGTAAGTGGTGATGAAGTGATGCCACCGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATAc                                                                                                                                                                                                                                <  2:3854/246‑1 (MQ=255)
      tAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGACCGATTTTACAGAAGCCTTAATGCCTAACATAACCTGGTGCGCCCTGCCTGAAGATGTCTCTTTATGGCCGGGCCTGCCTCTTTCATTAAGGGTTGATGAAGGGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCGGGATAAAAAGCGTTGGACCAAAAACCAGAGc                                                                                                                                                                                                                          >  1:345610/1‑250 (MQ=255)
               cAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTg                                                                                                                                                                                                                <  1:68209/251‑1 (MQ=255)
                   aGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCCGTATGGTCGGGGGCAGGATAAACAGCGGCTGACCCAATACCAGAGCAAACTGGGggcggc                                                                                                                                                                                                            >  2:15049/1‑251 (MQ=255)
                       aTGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTGCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATg                                                                                                                                                                                                        <  1:97224/251‑1 (MQ=255)
                                 ggtggctttagtcGTTATTTTCCCTGCGTCGAACGATTTTACAGGAGCCTTAAGGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGAttgtgg                                                                                                                                                                                              <  2:424883/239‑3 (MQ=255)
                                        tCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTg                                                                                                                                                                                       >  1:343223/1‑251 (MQ=255)
                                         cATGCTTTTTTTTCCCTTCGTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGCCCTGCCTGAAGATGTTTCTTTATGGCCGGGTCTGCCTCTTTCAGTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGTCGATGGTGATTGTTGCCGCCTg                                                                                                                                                                                       <  2:345610/250‑1 (MQ=255)
                                         cATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGACTTAATGCCTAACATTACCTGGTGCGACCGGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATGAAGGGGTGATGAAGTGATGCCCCTGGATTACCTCGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCTTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCTCCTgg                                                                                                                                                                                      <  2:358684/251‑1 (MQ=255)
                                                      cccTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCACACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGAtt                                                                                                                                                                         >  1:240607/1‑251 (MQ=255)
                                                               gAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAgg                                                                                                                                                                   >  1:8361/1‑248 (MQ=255)
                                                                            ggAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGt                                                                                                                                                                                 <  2:169735/221‑1 (MQ=255)
                                                                            ggAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGt                                                                                                                                                                                 >  1:169735/1‑221 (MQ=255)
                                                                             gAGCCTTAATGCCTAACATTACCGGGTGCGACCTGCCTGAAGATGTCTCTTTATGGACGGGTCTGTCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCTGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGt                                                                                                                                                  <  2:35358/251‑1 (MQ=255)
                                                                                  ttAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGGGATTGTTGCCGCCTGGTGCGTGGACGATTATCCGGTGATTCGTCGGGCAGGTtcact                                                                                                                                             >  1:468331/1‑251 (MQ=255)
                                                                                            aCATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAAc                                                                                                                                                                                                                      >  2:138842/1‑168 (MQ=255)
                                                                                            aCATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAAc                                                                                                                                                                                                                      <  1:138842/168‑1 (MQ=255)
                                                                                            aCATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGGTTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACgg                                                                                                                                   >  1:380620/1‑251 (MQ=255)
                                                                                              aTTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGATGGATAAACAGCGTCTGAACCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCTGGTTCACTCAACGCACGGGc                                                                                                                                 >  1:354239/1‑251 (MQ=255)
                                                                                                 aCCTGGTGCGTCCTACCTGAAGATGTCTCTATATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGCCGTGGGCTGGAACACCAGCGTATGACCCAATACCAGAGACAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGCGCGTGGAAGATTATCAGGTGATTTGTCTGGcgagttaacctcaagcccggtcgcc                                                                                                                              >  2:468337/1‑227 (MQ=255)
                                                                                                  ccTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGGGCGTGGAAGATTATCAGGTGATTCGTCTGGCCGGTTCACTACCCGCCCGGGCtccg                                                                                                                             >  2:157606/1‑248 (MQ=255)
                                                                                                     ggTGCGCCCTGCCTGAAGATTTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTacgcg                                                                                                                           <  2:287289/250‑4 (MQ=255)
                                                                                                                   gAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGc                                                                                                                                       >  1:147803/1‑224 (MQ=255)
                                                                                                                   gAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGc                                                                                                                                       <  2:147803/224‑1 (MQ=255)
                                                                                                                   gAAGATGTCGCTTTATGGCCGGGTCTGCCTCTTTCATGAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAgc                                                                                                            <  2:247783/251‑1 (MQ=255)
                                                                                                                    aaGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAgcg                                                                                                           >  1:377296/1‑251 (MQ=255)
                                                                                                                      gATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAgcgc                                                                                                          >  1:112013/1‑250 (MQ=255)
                                                                                                                         gTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCaac                                                                                                       <  2:112013/250‑3 (MQ=255)
                                                                                                                           ctctTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCCGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCCGGTTCACTCACCGCACGGGCTACGCGTCTGGCCCATGAATCGCAACTgg                                                                                                    >  1:224205/1‑251 (MQ=255)
                                                                                                                           ctctTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCGGGTGCGTGGAAGATTATCAGGTGATTc                                                                                                                                                             <  1:182076/194‑1 (MQ=255)
                                                                                                                           ctctTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGAGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTc                                                                                                                                                             >  2:182076/1‑194 (MQ=255)
                                                                                                                               tgaTGGCCGGGTCTGCCTCTTTCAGTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCCGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGTCCCATGAAGCGCAACTGGATg                                                                                                 <  2:439434/248‑1 (MQ=255)
                                                                                                                               gtATGGCCTGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATTCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCTCATGAAGCGCAACTGGAt                                                                                                  <  1:72363/248‑1 (MQ=255)
                                                                                                                                 aTGGCCGAGTCTGCCTCTTTCATTCAGTGGTGATGAAGTGCTGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGGGATTGTTGCCGCCTGGTGCGTGGCAGGTTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCCCGGGCTACGCGCCTGGCCCATGAAgccacactggctgtcc                                                                                              >  2:189948/1‑237 (MQ=255)
                                                                                                                                                tctTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGAGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGGCCCGCTGGGGaaa                                                                                >  1:268837/1‑250 (MQ=255)
                                                                                                                                                   ttCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAAt                                                                              <  1:157606/249‑1 (MQ=255)
                                                                                                                                                     cATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATccc                                                                           <  1:82185/250‑1 (MQ=255)
                                                                                                                                                        tAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCAc                                                                       >  2:441244/1‑251 (MQ=255)
                                                                                                                                                                 atgaAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACg                                                              >  2:166650/1‑251 (MQ=255)
                                                                                                                                                                     aGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCgg                                                          <  1:263466/251‑1 (MQ=255)
                                                                                                                                                                      gTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCggg                                                         >  1:240153/1‑251 (MQ=255)
                                                                                                                                                                       tGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGt                                                        >  1:3722/1‑251 (MQ=255)
                                                                                                                                                                            ccACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGCTTGCt                                                   >  2:96114/1‑251 (MQ=255)
                                                                                                                                                                             cACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTg                                                  <  1:166650/251‑1 (MQ=255)
                                                                                                                                                                                     tACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATgg                                           >  1:59888/1‑250 (MQ=255)
                                                                                                                                                                                        cACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTTATGGatat                                       <  2:240153/251‑1 (MQ=255)
                                                                                                                                                                                              ggTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGattc                                 <  1:325873/251‑4 (MQ=255)
                                                                                                                                                                                                  gTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCGGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                             <  1:441244/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAt         >  1:352433/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTCGTTGGCTTGATAAGCATCGTCTGACACAATACTATAGCATACTGGGTGCGGCTATGGTGATTGTTGCCGCTTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTCCACGCACCGCACGGGCTACGTGCCTGGCCCACGAAGCGCAACTGTGTGTCGCCCCGCTGGGGAAAAGCTCGCACCTACGCACGCCGTGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAtt        <  2:7914/251‑1 (MQ=255)
                                                                                                                                                                                                                            gggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGAATGAATGTATtgatga   <  1:374444/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGaa  <  2:3722/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
GAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAA  >  NC_000913/4624810‑4625281

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: