Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,490,095 G→A intergenic (+424/‑807) aldA → / → cybB aldehyde dehydrogenase A/superoxide oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,490,0950GA93.5% 123.3 / ‑6.1 48intergenic (+424/‑807)aldA/cybBaldehyde dehydrogenase A/superoxide oxidase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (22/21);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTG‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTTTTCAT  >  NC_000913/1489848‑1490338
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
cgaaccgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTa‑ggg                                                                                                                                                                                                                                                   <  1:501312/251‑1 (MQ=255)
              tGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTg                                                                                                                                                                                                                                     >  1:361117/1‑251 (MQ=255)
                                 gaCTATTTCGTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCTTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACg                                                                                                                                                                                                                  <  1:523807/251‑1 (MQ=255)
                                  actctttCTTCATCAGTATAACCAAATGACTCGTTGTTGGGGGTTGCTTGTTTAATAGCGTTATTCACCTATTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGACCCTGTTTTCACCGGCACGAGTTGCGCCTGATCTTTCGGTTGGCCGCTCAGTTCCGGGATCACCAGACCACTGGCTTTTGACGCCCCCGTAGTGTA‑GGGAATACTAGTTTCTGCAGATGCGCGGGAAGCActt                                                                                                                                                                                                                 >  2:497098/5‑249 (MQ=255)
                                  aCTATTTCGTCATCGGTATAAACAAATGACTCGGTATTGGTGGTTGCTTGTTTAAGTGCGTTAGTCACCTCTTCGGCAGTCACTTTTTTCCGAGAAACGATTCCAGTTCAGTGACCGTACCGGTTTTCCCCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTGCCGGGATCACCAGATCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGt                                                                                                                                                                                                                 <  2:41267/251‑1 (MQ=255)
                                               cGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCCGACCAATGGCTTTTGCCGCCCCAGTAGTGTA‑GGGAATGCTATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCAc                                                                                                                                                                                                    >  1:107608/1‑251 (MQ=255)
                                                ggaataACCAAATGACTCGTTATTGGTGGGTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGTTATTTTCTGACGCTGCGCGTCAAGCACGTAAATCTTTACCAcc                                                                                                                                                                                                   >  2:535918/4‑250 (MQ=255)
                                                ggTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTG‑GGGAATGATATTTTCTg                                                                                                                                                                                                                                     >  2:175654/1‑217 (MQ=255)
                                                   ataACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATCCCAGTTCAGTGGCCGAACCTGTTTTCATCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCATACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCGGCGCGTGAAGCACGTAAATCTTTACCACGCgg                                                                                                                                                                                                <  1:586402/251‑1 (MQ=255)
                                                   ataACCAAATGAATAGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTTACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCAcgcg                                                                                                                                                                                                 <  1:660760/250‑1 (MQ=255)
                                                     aaCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGc                                                                                                                                                                                              <  1:551689/251‑1 (MQ=255)
                                                                  gTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGa                                                                                                                                                                                 <  1:314447/251‑1 (MQ=255)
                                                                     aTTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTg                                                                                                                                                                              <  1:294246/251‑1 (MQ=255)
                                                                       tggtggTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTggg                                                                                                                                                                            >  1:224287/1‑251 (MQ=255)
                                                                       tggtgaTTGCTTGTTTAAGTGCGTTATTCAAATCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTAAATGACCGACCCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGTATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTC‑GGGACGGATATTTTCTTCCGCTGCGCGTGAAGCACGTAAATCTTTAACAGGCGGGcccaccaaccatgaatgcg                                                                                                                                                                            >  2:437554/1‑234 (MQ=255)
                                                                           ggtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGcc                                                                                                                                                                        >  1:567799/1‑251 (MQ=255)
                                                                                  tGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGtat                                                                                                                                                                   >  1:211000/1‑249 (MQ=255)
                                                                                       aaGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATg                                                                                                                                                           <  1:328849/251‑1 (MQ=255)
                                                                                            cGTTATTCACCTCTTCGGCAGTAACCTTTTTCCGAGAATCGATCCCAGTTCAGTGCCCGAACCTGTTTTCCCCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATgg                                                                                                                                                      <  1:125936/251‑1 (MQ=255)
                                                                                              ttATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTACTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGgtc                                                                                                                                                    >  1:521316/1‑251 (MQ=255)
                                                                                                   ccccTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCAt                                                                                                                                               >  1:265872/3‑251 (MQ=255)
                                                                                                        ctTCGGCAGTCACTTTTTTCCGAGAAGCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGt                                                                                                                                            <  1:527242/249‑1 (MQ=255)
                                                                                                             gCAGTCACTTTTTTCCTCGAACCGATACCAGTTCAGGGACCTAACCTGTTTTCACCGGCACGGGTTGCGCATGACCTTTCAGGTTGCCGCTCAGTGCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTAGAGG‑CATGAATGGTCGTCATCGTGCCGAcg                                                                                                                                     <  2:658321/251‑2 (MQ=255)
                                                                                                               aGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CTTGAATGGTCGTCATCGTGCCGACTTc                                                                                                                                   >  1:310550/1‑251 (MQ=255)
                                                                                                                gTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCt                                                                                                                                  <  1:535918/251‑1 (MQ=255)
                                                                                                                           ccGAGAATCGATACCAGGTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGCTATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCCGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGGGCCGACTTCTATCCCAAAACt                                                                                                                       >  2:50295/1‑251 (MQ=255)
                                                                                                                               gaATCGATTCCAGTTCAGTGTCCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTGGCCGCTTAGGTCCTGGATCACCAGACCAATGGCTAATGCCGCCCCCTTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTAAAGCACGTAAATCTTTTCCACGCGGGCGATCCTCCATTGTCTGGATGCGAGTATAGG‑CATGAATGGTCGTCATTGTGCCGACTTCTAACCCGGAACTGTCa                                                                                                                   <  2:331372/251‑1 (MQ=255)
                                                                                                                                   cGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATTACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCa                                                                                                               <  2:361117/251‑1 (MQ=255)
                                                                                                                                     aTACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAGGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGGGCCAGTATAGGACA‑GATTGGGCGGCATCGAGCCCACTTCTATCCCGAAACTATCATGCa                                                                                                               >  2:317897/1‑249 (MQ=255)
                                                                                                                                                gTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTCTAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCAt                                                                                                  >  1:550515/1‑251 (MQ=255)
                                                                                                                                                  gtccGAACCTGTTTCCACCGGCTCTCGTGGCGCATGACCTTTCAGTTTGCCGCTCAGTGCCGGGATCCCCAGACCAATTGCTTTTTCCGCCCCCGTAGTGTA‑GGGAATGATAGTTTCTGCCTCTGCGCGTGAAGCACGTAAATCTTTACCCCGCGGGCCATCCACCAGTGCCTTGGTGCCAGTATAGG‑CAGGAATGGTCGGCATCGTGCCGATGTCGATTCCGAAACTGTTATGCAATGCTTTGGCCAta                                                                                                 <  2:401433/248‑2 (MQ=255)
                                                                                                                                                   aCCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGACCACCAGACCAATGGCTTTTGCCGCCCCCGTCGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGGGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCAGCGTGCCGGCTTCTCTCCCGAACCTGTCATACAAGGCATTGGCCATCgg                                                                                               >  2:41133/1‑251 (MQ=255)
                                                                                                                                                    ccGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCC‑GTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCCTTACCACGCGGGCCCACCAACAGTGACTGGGTGCCAGAATAAG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGgcg                                                                                             >  1:397429/1‑251 (MQ=255)
                                                                                                                                                     cGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGgcg                                                                                             >  1:521870/1‑251 (MQ=255)
                                                                                                                                                              tttCACCGGCCCGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAg                                                                                     <  2:567799/250‑1 (MQ=255)
                                                                                                                                                                tCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGtt                                                                                   <  1:672878/250‑1 (MQ=255)
                                                                                                                                                                   ccGGCACGAGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGGCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGGGACTGGGGGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCACGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAgtg                                                                               >  2:519768/1‑251 (MQ=255)
                                                                                                                                                                           cgTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCAAGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCTACTTCTATCCAGAAACTGTCATGCAAGGCTTTGGACATCGGCGCAAGAACGTTAGTGGTGcctg                                                                        >  2:80056/1‑247 (MQ=255)
                                                                                                                                                                                                      cgtatgtCCGGGATCACCAGGCCAATGGCGTTTGCCCCCCCCGTAGTTTC‑TGGAATGATATTTTCTGCCCCGGCGCGTGAAGCACGTAAATCTTTTCCTCGCGGGACATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTTGTCAACGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTGAGGGGTGCATGACGCCACGTAAACAATGGTGTCTTTGcc                                            <  2:265872/245‑1 (MQ=255)
                                                                                                                                                                                                        caggcccggaataacaaGACCATTGGCTTTTACCGCCCCCGTAG‑GTGAGGGAATGATATTTTCTGCCGCGGCGCGTGAAGCTAGTAAATCTTTGCCACGCGTGCCATCCACCAGGGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGGCTTCTATCCCGAAACTGACATGCAAGGCTTTGGCCATCTGCGTAATACAGTTAGTGGTGAATGACGCCTCGGAAACAATGGTGTCGTTGcc                                            <  2:70295/233‑1 (MQ=255)
                                                                                                                                                                                                           ttCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGAGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTAATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGGGGGGCATGACGCCACGGAAAAAATGGTGTCGTTGcaaaca                                        >  2:361678/1‑245 (MQ=255)
                                                                                                                                                                                                            tCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATcc                                        >  1:321399/1‑249 (MQ=255)
                                                                                                                                                                                                                      ccAGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTgtcgtc                             >  1:598670/1‑250 (MQ=255)
                                                                                                                                                                                                                        cgaccatgggCTTTTTCCGCCCCCTTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCACCCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGGCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCAtt                          <  2:521316/243‑1 (MQ=255)
                                                                                                                                                                                                                        aGACCAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGGCTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGAc                                                                                                                                      <  1:592387/143‑1 (MQ=255)
                                                                                                                                                                                                                           ccAATGGCTTTTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGAc                       <  1:636687/251‑1 (MQ=255)
                                                                                                                                                                                                                           ccAATGGCTTGTGCCGCCCCCGTAGTGTA‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGAc                       <  1:252829/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggCTGTTGCCGCCCCCGTAGGATG‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTta                   <  2:521870/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAGTGTA‑GGGAATGATATTTTCTGACGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGTCCATCCACCAGTGACTGGGGGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGCAACCATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTtttcat  >  2:223474/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTG‑GGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGG‑CATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTTTTCAT  >  NC_000913/1489848‑1490338

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