Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 365,783 T→C 25.8% S175G (AGC→GGC)  lacZ ← beta‑D‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913365,7830TC25.8% 96.6 / 5.6 66S175G (AGC→GGC) lacZbeta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base T (17/29);  new base C (16/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 8.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GAACTTCAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGG  >  NC_000913/365570‑366011
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
taaCTTCCGCCTCCAGTACAGGGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGa                                                                                                                                                                                                                         <  1:267979/228‑1 (MQ=255)
agatcgccgcttCCAGTACACCGCGGATGAAATCAGCATGAAAGCGAGTGG‑AACACTGGAAATCGCTGATTTGGGTAGGCGGTTTAGGCAGCAACGCGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATATCTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGa                                                                                                                                                                                                                         <  2:413250/218‑1 (MQ=255)
      cAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTGGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCCTAAAAACGCACCCAGGCCAAATTCaacagggaa                                                                                                                                                                                                                   >  2:320230/1‑221 (MQ=255)
      cAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGCTAACTGCCGTCACTCCAGCGCACCACCATCACCGCGAGGCGGTTTTCCCC‑GGCGCGTGAAAATTCCCCCAGGGCAAAATCAAACGCCaa                                                                                                                                                                                                                   >  2:413154/1‑229 (MQ=255)
        gCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCGACACCATCACCGCCAGGCGGTTTTCTCC‑GGGCCGCAAAAATAAGCCCAGGCCAAATTAACACGACAAAc                                                                                                                                                                                                                 >  2:280528/1‑229 (MQ=255)
        gCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCCTGGCGCTTAAAAA‑GCCCCCCGGTCAACTTCAAACGGCaaag                                                                                                                                                                                                                 >  1:735592/1‑228 (MQ=255)
         ccTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCCCTCCAGCGCAGCACCCTCACCGCCAGGCGGTTTTCTCC‑GGCCCGTCAACATGCGCTCAGGTCAAATTTCGCCGGCAAAc                                                                                                                                                                                                                 >  1:89892/1‑228 (MQ=255)
          tcccAGTCCAGCCCGGCTCAAATCATCCATAAAGCGAGTGGCCACA‑TGGAAACCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGa                                                                                                                                                                                                               <  1:179554/227‑1 (MQ=255)
          gtccagtccagCGCGGTGTAAATCCTCATTAAAGCGAGTGGCAACA‑GGGAAATCGCGGATTTGTGTAGTCGGTTTAAGCAGCAACGTGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGa                                                                                                                                                                                                               <  2:721385/228‑1 (MQ=255)
               gcacggctcgGCTTAAATCATCAGTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCCCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑Tc                                                                                                                                                                                                         <  1:714874/221‑1 (MQ=255)
               atacccgcgggttGAAATCCTCTATAAGGGGGGGGGCAAAA‑TGGAAATCCGTGATTTGTGTAGTCGGTTTATGCAGCACCGAGACGTCACGGAAAATGCCGCTCATCCCCCACATATCCTGATCTGCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑Tc                                                                                                                                                                                                         <  2:298018/217‑1 (MQ=255)
                     cgcgGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGc                                                                                                                                                                                                    >  1:18581/1‑228 (MQ=255)
                     cgcgGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCCCGTAAAAATCGGCCCAGGCCAAATTCGCACGGCAAAC‑ACTGTTCCTGGc                                                                                                                                                                                                    >  1:460508/1‑228 (MQ=255)
                       cgGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCCCGTAAAAATGCGCTCAGGTCAAATTCCGACGGCAAACGACCT‑TCCTGGCCg                                                                                                                                                                                                  >  1:203621/1‑228 (MQ=255)
                          gttaaaTCTTAATGAAACCGAGGGGAAAGA‑GGGAAATCGCTAATTTGTGTATTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAAc                                                                                                                                                                                              <  1:327848/226‑1 (MQ=255)
                          cTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATCCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTaa                                                                                                                                                                                               >  1:695910/1‑228 (MQ=255)
                          cTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAAAGCGCTCAGGCCAAATTCAGACGGCAAACGGCGT‑TCCTGGCCGGaa                                                                                                                                                                                               >  1:680241/1‑228 (MQ=255)
                          cTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGGAAAAATGGGCTCAGGTCAAATTCCGACGGCAAACTACTG‑TCCTGGCCGTaa                                                                                                                                                                                               >  1:557647/1‑228 (MQ=255)
                          cTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑CGCGCCTAAAAATGCGCCCAGGTCAAATTCCGACGGCAAACGGCCG‑TCCTGAGCGTaa                                                                                                                                                                                               >  1:367473/1‑228 (MQ=255)
                          cTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCACCACCATCACCGCGAGCCGGTTTTTTCC‑GGCACGTAAAAATCCGCCCAGGCCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGtca                                                                                                                                                                                               >  1:323501/1‑226 (MQ=255)
                                            cGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGt                                                                                                                                                                                                                                   <  1:603972/175‑1 (MQ=255)
                                            cGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGt                                                                                                                                                                                                                                   >  2:603971/1‑175 (MQ=255)
                                             gggTGGCAACA‑TGGAAAACGATGATTGGGGGGGTCGGTGTAGTAACGAACGAGACGTCCCGGAAAATGGCGCACACCCGCCCAAAATCCGGATCTGCCAGATAACCGCCGTCACTACAGCGCAGCCCCATCACCGCGAGGCGGTTTTCGCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGca                                                                                                                                                                          <  2:411910/227‑1 (MQ=255)
                                              aGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑AGCCCGTAAAAATGCGCTCAGGTCAAATTCAGACGGGAAAAGCCTT‑TCCCGGCCGTAACCGACCCA‑CGCACCGTTGccc                                                                                                                                                                         >  1:646813/1‑227 (MQ=255)
                                                              aTCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCCCAGGTCAAATTCAGCCGGCAAACGACTG‑TCCTGGCCGTCACCGACCCACCGC‑CCCTTGcacccacaatcaaacccc                                                                                                                                                          >  1:126880/1‑215 (MQ=255)
                                                                   ttaTTTGTGTAGTCGGTTTTTGCGGAAACGAAACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTACAGCTAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGTCGGTTTTATCC‑AGCGCGTAAAAATGCGCTCAGGTCAACTTCAGACGGCAAACGCCTG‑TCCTGGCCGCCACCAACCCAGCGC‑CCGTCGACCCACACATAAAACAcccaagg                                                                                                                                                     >  2:64175/3‑224 (MQ=255)
                                                                     aTTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCCCATATCCTGAGCTTCCAGATAACTGCCGTCGCTCCAACGCACAACCATCACTCCGAGGCGGTTTTTGCC‑GGCGAGTCAAAATGCGCCCAGGTCCAATTCAGACGGCAAACCCGTT‑TCCTGGCAGTAACCCACCCAGATC‑CCGGTGCCCaacaaagcaaacacgcgttgaa                                                                                                                                                   >  2:316067/1‑211 (MQ=255)
                                                                     aTTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGCGGCGGTTTTCTCC‑GTCGCGTCAAAATGCGCCCAGGGCAGATTCCGACGGCCAACAACCT‑TCCCGGCCGTAACCCACCCCCCGC‑CCCTTGGCCCCCAGaaaaaaacccccaataa                                                                                                                                                   >  2:419089/1‑213 (MQ=255)
                                                                               tCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGCAAAAATTCGCTCAGGGCAAATTCAGACGGCCAACGCCTG‑TCCTGGCCGTAAACGACCCAGCCC‑CCGTTGCACCACAGGTGAAACGCCGCGTTAACGcactcca                                                                                                                                          >  1:242830/1‑222 (MQ=255)
                                                                                 ggTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCCCCTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCACCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACACCCACaaaa                                                                                                                                        >  2:59215/1‑228 (MQ=255)
                                                                                    ttATGCAGCAACGAGACGTCACGCAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCCCCATCACCGCGGGGCGGATTTCTCC‑GGCACGGCACAACGCGCCCACGCCACATTTCGACGGCAAAAGGCTT‑CCCGGGCCGGAACACACCCACCCA‑CCGCTGTCCCACACACGAACCGCCCAGGTAACCCCCCCAAAaataa                                                                                                                                    >  2:482213/1‑229 (MQ=255)
                                                                                         cagcaACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACg                                                                                                                                                 <  1:88019/211‑1 (MQ=255)
                                                                                         cagcaACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGGAAAAATGCGCTCAGGGCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGGCCCAGCAC‑CCGTTGCACCACAGATGAAAACCCGCGTTAAACCCCACaaaaaacttgcc                                                                                                                                >  2:88019/1‑221 (MQ=255)
                                                                                           gcaACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCCCAGGTCAAATTCAGACGGCAAACCACTT‑TCCTGGCCGTAACCCACCCAGCGC‑CCGTTGCCCCCCCGACGAAACGCCGGGTTAACCCCCACAAAAATAAATCGCGt                                                                                                                             >  2:225362/1‑229 (MQ=255)
                                                                                             aaCGAGACGTCACTGAAACTGCCGCGCATCGCCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTc                                                                                                                            <  1:570995/228‑1 (MQ=255)
                                                                                                 gtccctcccggacacgtccccTCATCCGCCACATCTCCCGTACTTCCAGCTAACTGCCGTCACTCCAGGGCAGCCCCATCACCGCGAGGCGGTTTTCCCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑CCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGcc                                                                                                                       <  2:121858/209‑1 (MQ=255)
                                                                                                 agccGTCACGGAAAATGCAGCTCATCCGCCACATCTCCTGATCTTCCAGATACCTGCCGTCACTCCAGCGCAGCACCATCACCGCGCGGCGGTTTTCTCC‑GGCGCGTAACAATCCGCACCGGGCAACATCAGACGGCAAACGACCG‑GCCTCGCCGTATACGCCGCACAGC‑CCGTTGCACCCCAGATGAGACGCCGAGTTACAGCCATCCAAAACAAATCCCGTCTggac                                                                                                                       >  2:225276/4‑227 (MQ=255)
                                                                                                         cGGAAAATGCCGCTCATCCGCCACATATCATGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTAC‑GGCGCGGCAAAATGCGCCCAGGTCCAATTCAGCCGGCCACCCACTG‑GCCCGGCCCAAACCGCCCCCCCCC‑CCCATTCCCCACAAATGAAACCCCGAGTTACCCCCCTCCAAAATCAtttccgccccgccctcctcc                                                                                                                >  2:588083/1‑210 (MQ=255)
                                                                                                               aTGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCCCCGATGAAACCCCGGGGTAACGCCCTCAAAAATAATTCCCGCCTGGCCTTCCCGTAGCCAGc                                                                                                         >  1:366442/1‑229 (MQ=255)
                                                                                                                      tCATCCGCCACATAGCCTGATCGTCCAAATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTtcatc                                                                                                  <  1:459838/229‑1 (MQ=255)
                                                                                                                      tCACCCCCCCCATACCCTAATCTTCCAGATACCCGCCGTCACGCCAGCGCAGCCCCATCCCCGCGAGGCGGTTTTCCCC‑GGCGCGTAAAAATGCGCTCAGGGCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑ACGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTtcatc                                                                                                  <  2:494033/229‑1 (MQ=255)
                                                                                                                       tagccGCCACATATCCTGACCGTCCAGATAACTGCCGTCACTCCAGCGCAGCCACAACACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATACAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTtcatc                                                                                                  <  1:740995/225‑1 (MQ=255)
                                                                                                                       cctCCGCCACATAGCCTGGTCTTCCAGAAAGCTGCCGTCACTCCAGCGCAGCACCATCCCCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTtcatc                                                                                                  <  1:280528/226‑1 (MQ=255)
                                                                                                                       cATCCGCCCCATATCCTGATCTTCCAGATAACTGCCGTCCCTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTtcatc                                                                                                  <  1:225362/228‑1 (MQ=255)
                                                                                                                       cATCCGCCACATACCCTGATCTCCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTtcatc                                                                                                  <  1:413154/228‑1 (MQ=255)
                                                                                                                          cccccccatctcccggcctcccAGTTAACTGCCGTCCCACCAGCGCGGCACCATCACCGCGAGGGGGTTTCCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACa                                                                                              <  2:680240/209‑1 (MQ=255)
                                                                                                                                      cTGATCTTCCAGATAACTGCCGTCCCTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGc                                                                                  <  1:59215/229‑1 (MQ=255)
                                                                                                                                        gATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATTCGCCCAGGGCAAATTCAGCCGGCAAACGACTG‑TCCTGGCCTTAACCTACCCAACCC‑CCGCTGCACCACAGATTACACCCCGGATTAACACCCTCAAAACAAATTCACCCCTGGCCTTCCCGTTCCCACATTTCATCAACATTAAATTGGAGCg                                                                                 >  2:513767/1‑228 (MQ=255)
                                                                                                                                        gATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑GCCTGGCCGTAACCGACCCAGCGC‑CCGTTGAACCACAGATGAAACGCCCAGTTAACGCCACCAAAAATAATTCCCGTCCGGCCTTCCCTTAGCCACCCTTCCTCACAATTAAATTTGAGCGa                                                                                >  2:10244/1‑229 (MQ=255)
                                                                                                                                              ccAGATACCCGCCGTCCCACCAGCGCACCACCAACACCGCGCGCCGGTTTTCCCG‑GGCGCGAAAAACTGCGCTCAGGTAAAATTCAGACGGCAAACGCCTT‑GCCTGGCAGTAACCGACCCAGGCC‑CCGTTGCCCCACAGATGAAACGTCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTaaca                                                                          <  2:581221/229‑1 (MQ=255)
                                                                                                                                               cAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGGCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAACAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGGGTaacaa                                                                         >  1:441162/1‑229 (MQ=255)
                                                                                                                                                                          caCCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCCAATTCAGACGGGAAACGACTG‑GCCTGGCCGTAACCGACCCAGCGC‑CCCTTGCACCCCAGATGAAACGCCGAGTTAACGCCATCAAAAATACTTCGCCTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTACCAACCCGCCGGATTCTCCGCGGGAACaaa                                               >  2:553501/1‑228 (MQ=255)
                                                                                                                                                                          caCCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACCCCATCAAAAATAATTCCCGTCTGGCCTTCCTTTAGCCAGCTTTCATCAAAATTAAATGTGAGCGAGTAACAAACCGTCGGATTCTCCGTGGGGACaaa                                               >  1:321274/1‑228 (MQ=255)
                                                                                                                                                                          caCCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCCCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGAAGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACAATAAAAGTGAGCGCGAAAAAACCCCTCGAATTGCCCGGGGGaaaaaa                                               >  2:724888/1‑228 (MQ=255)
                                                                                                                                                                            ccATCGCCCAGAGGCGGTTTTCACC‑GGCGCGTAAAAATGCGCTCAGGTCAAATGCAGACGGCAAACGTCTG‑TCCTGGCCGTAACCGACCCAGCCC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAAcgg                                            <  1:284793/229‑1 (MQ=255)
                                                                                                                                                                                ccacgtgacggggTTTTCACC‑GGCGCGTAAAAATGCGCTTAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGa                                        <  1:748658/218‑1 (MQ=255)
                                                                                                                                                                                  ccGCGGGGCGGTTTTCTCC‑GGCGCTAAAAACGGCGCTCAGGTCAAATTCAGACGACAAACGATTG‑GCCGGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGGGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGAtt                                      <  2:367473/229‑1 (MQ=255)
                                                                                                                                                                                     cgAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGAc                                   <  1:10244/229‑1 (MQ=255)
                                                                                                                                                                                       aGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCCCAGGGCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCCC‑CCGTTGCAACACAGATGAAACGCCGAGTTAACGCCATCAAAAATCATTCGCGTCCGGCCTTCCTGCAGCCAGCTTTCATCAACATTAAATGTGAGCGGGGAAAAACCCCGCGGATTCCCCCGGGGAACaaaagcgggtatgaacc                                 >  2:489672/1‑215 (MQ=255)
                                                                                                                                                                                          gggTTTTCCCC‑GGGGCCGAAAAAGGCCGTCAGTTCAAATTAAGCCGGCAAACGACTG‑TCCTGGCCGCAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGTGAACAAACGGCGGATTGACCGTaa                              <  2:427460/228‑1 (MQ=255)
                                                                                                                                                                                                ttccccgcgcgatAAAAATGCCCCACGGGCAAATTCAGACGGCAAACGACGG‑TACTGGCCGTAACCGCCCCAGCGC‑CCGTTGCCCCACAGATGAAACGCCGAGTTAACGCCAACAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGTTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGa                         <  2:366442/217‑1 (MQ=255)
                                                                                                                                                                                                      ggCGGGTAAAAATGGGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGt                    <  1:588084/228‑1 (MQ=255)
                                                                                                                                                                                                      ggCGCTTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑GCCTGGACGTAACCGACCCAGCGC‑CCGTTGCACCACAGATTAAACGTCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGt                    <  1:320230/228‑1 (MQ=255)
                                                                                                                                                                                                       gccggTAAAAAGGCCCTCAGGTCAAATTCAAACGGCCACCGACTG‑TCCTGGCCGTAACCGGCCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCa                  <  2:246591/225‑1 (MQ=255)
                                                                                                                                                                                                        cgcgTAAAAATGCGACCAGGTCAAATTCAGACGGCAAACGACTG‑ACCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCa                  <  1:246409/228‑1 (MQ=255)
                                                                                                                                                                                                         gcttAAAAAGGCGGCAAGGCCAAATTCAGACGGCAAACGACTG‑TCCGGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGGCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACg                <  2:18581/226‑1 (MQ=255)
                                                                                                                                                                                                                       tCAGGTCAAATTCCGACGGCAAACGACTG‑TCCGGGCCATCACCGACCCCGCGC‑CCGTTGCACCACAGATGAACCGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACACGTCGGATTCTCCGTGGGAACAACCGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATggg  <  1:648181/229‑1 (MQ=255)
                                                                                                                                                                                                                       tCAGGGCAAATTCAGCCGGCAAACGACTG‑TCCGGGCCGTAACCGACCCAGCGC‑CCGTTGCACCAAAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATggg  <  1:419089/229‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
GAACTTCAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACA‑TGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCC‑GGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTG‑TCCTGGCCGTAACCGACCCAGCGC‑CCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGCGTCTGGCCTTCCTGTAGCCAGCTTTCATCAACATTAAATGTGAGCGAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGG  >  NC_000913/365570‑366011

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