Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,589 G→A S60S (TCC→TCT)  pyrE ← orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,5890GA89.6% 144.5 / ‑6.5 49S60S (TCC→TCT) pyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/30);  minor base C (5/0);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 5.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGT‑TGTGGTGGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCGGAATCCACCAACGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAA  >  NC_000913/3815344‑3815832
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATc                                                                                                                                                                                                                                                 <  1:231518/251‑1 (MQ=255)
            cccATCGATTCACGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCCCCGTGGTCTTTTGCCGCTTTGCGGGTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAACCCACCAATGCTTc                                                                                                                                                                                                                                     <  2:61416/250‑1 (MQ=255)
             ccATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTc                                                                                                                                                                                                                                     <  1:104180/250‑1 (MQ=255)
               aTCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTcggc                                                                                                                                                                                                                                  <  1:154842/251‑4 (MQ=255)
                      cgcgAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑Taaaaa                                                                                                                                                                                                                          <  1:149115/251‑1 (MQ=255)
                               ccGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAAc                                                                                                                                                                                                                 <  1:434120/251‑1 (MQ=255)
                                                         aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTa                                                                                                                                                                                         <  1:168885/249‑1 (MQ=255)
                                                               aTTACGCGTCCTTGTAACGCGCTACCAACCAGATGGCCGCCGTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTGCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGccc                                                                                                                                                                                 <  2:286113/251‑1 (MQ=255)
                                                                 tACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGG‑AGTGGTCGCAATCGGGATACCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATTCCCGAATCCCCCCATGCTTCCACG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCCGGATTAAACAGCCCgg                                                                                                                                                                               >  1:440958/1‑251 (MQ=255)
                                                                      gTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTg                                                                                                                                                                          <  1:347966/251‑1 (MQ=255)
                                                                           tGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaaga                                                                                                                                                                      <  1:365302/250‑1 (MQ=255)
                                                                                   gcTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGAGCCACGGCGGG‑AGGGGGCGGAATCCGGGTCCCTTTGTAACGC‑GGGGCAAAAAGCAGATCGGACTCA‑ATGCCCGAATCCAACAATGCTTCAGCG‑GAAAAACGGCCCAAAAGGGCCAGGGGGCGCCCGGGGTTAAAAAGCCCGGCGTTGAAAAAATAggggg                                                                                                                                                             >  1:426287/1‑250 (MQ=255)
                                                                                     tACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGG‑AGGGGTAGCCCTCGGGATACCTTTTTAA‑GCCGGGCCCAACACCAGATCGAA‑ACACATGCCAGAATCCACCCACGCTTCAACG‑TAAAAAACGCCTAACAGTGACACGGCGCGCCCGGTATTAAACAACCCCGCGTTTAAGAAAAAgggggct                                                                                                                                                           >  2:260458/1‑248 (MQ=255)
                                                                                              gATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCg                                                                                                                                                  <  1:232273/251‑1 (MQ=255)
                                                                                               aTTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGG‑AGTGGGCGCCATCGGGATCCCCTTGTAA‑GCAGGGCCAAACAGCAGATCGAACCCC‑ATGCCCAAATCCACCAATGCTTCAGCG‑TAAAAAAGGGCTAACAGAGCCCGGGCGCGCCCGGGATTAAAAAACCCGGCGTTGAAAAAAAAGGGGGTTTTGGGCCCgg                                                                                                                                                 >  2:331783/1‑251 (MQ=255)
                                                                                                     gccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCATGATGCTCCGCCCGCGCCACGGCGGT‑AGTGGTGGCACACCGGATCCCTTTTTAA‑CCCGGGCCCAACAGCCGATCCAACCCC‑ATGCCCGAATCCACCCATGCTTCACCGGAAAAAACGCCCTA‑CAGTGCCCGGGCACCCCCGGGATTACACAACCCGGCGTTGGAGAAAAAAGGGCCTTTGCGCCCGGGttttt                                                                                                                                           >  2:392618/1‑249 (MQ=255)
                                                                                                         tCACCGTGGTCTTTTGACTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGACACGGCGGT‑AGGGGTCGCAAAAGGGGTCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATACACCAATGCTTCAGCG‑TAAAAACGGCCTAAAAGTGCCAGGGCGCCCCCGGCATTAAACAGCCCCGCGTTGAAGAAATAAGGGCTTTTGCGCCCGGATTTCAGCgg                                                                                                                                       >  1:89320/1‑250 (MQ=255)
                                                                                                                     tttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCa                                                                                                                             <  1:314314/249‑1 (MQ=255)
                                                                                                                                gCGGTTCAAGCAGTCCGGCAGGTCAAGGTAATGATGCGCCGCCAGCGCCGCGGCGGT‑AGTGGTCGCAAGCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCCCCAATGCTTCAGCG‑TAAAAACGGCCTAACATTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGTCAAACTTTAACCCCt                                                                                                                <  2:266723/251‑1 (MQ=255)
                                                                                                                                gCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGTCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGGTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACAcc                                                                                                                 <  1:438345/250‑1 (MQ=255)
                                                                                                                                        cccAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAATCTTTAACACCtgcttgct                                                                                                         <  2:415283/248‑1 (MQ=255)
                                                                                                                                        aGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTa                                                                                                        <  1:331783/251‑1 (MQ=255)
                                                                                                                                                ggCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCaa                                                                                                 >  1:486948/1‑250 (MQ=255)
                                                                                                                                                  caggtcaaggtcATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCCGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAtt                                                                                              >  2:413284/1‑251 (MQ=255)
                                                                                                                                                        aaggtcATGATGCTCCGCCAGCGCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAAAGGCCTAACAGTGCCAGGTCACGCCCGGAATTAAACAGCCCGGCGTTGAAGAAATAGGGGCGTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTCAACACCTGCTTGCTAAGCGCAAATTCAATa                                                                                         >  2:406507/2‑250 (MQ=255)
                                                                                                                                                                      ccGCCAGCGCCACGGCGGTATG‑GGTCGCAATCGGGATCCCGTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATGTCAGCGTAAACTCGCGAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGat                                                                            <  2:435929/249‑1 (MQ=255)
                                                                                                                                                                      ccGCCAGCGCCACGGCGGG‑AGTGGTCGCAATCGGGATCCCTATGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ACCCCAGAAGCCCCCAACGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCTGATTTCAGCGTAAACTCGCCAAACTTGAACACCTGCTTGCTGAGCGCAAATTCAATAAACTTGCGCTGat                                                                            <  2:21061/249‑1 (MQ=255)
                                                                                                                                                                           agcgCCACGGCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGttt                                                                     >  1:495371/4‑251 (MQ=255)
                                                                                                                                                                                   ggCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACt                                                                                                                         <  2:286010/190‑1 (MQ=255)
                                                                                                                                                                                   ggCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACt                                                                                                                         >  1:286010/2‑191 (MQ=255)
                                                                                                                                                                                    gCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCaat                                                                                                                                                                                                                                          >  2:325287/1‑77 (MQ=255)
                                                                                                                                                                                    gCGGT‑AGTGGTCGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCaat                                                                                                                                                                                                                                          <  1:325287/78‑2 (MQ=255)
                                                                                                                                                                                                cgCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAtt                                                                                              >  2:326953/2‑206 (MQ=255)
                                                                                                                                                                                                cgCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAtt                                                                                              <  1:326953/205‑1 (MQ=255)
                                                                                                                                                                                                cgCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATc                                                  <  1:111471/249‑1 (MQ=255)
                                                                                                                                                                                                cgCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATc                                                  <  1:190797/249‑1 (MQ=255)
                                                                                                                                                                                                 gCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCtt                                                <  1:219928/251‑1 (MQ=255)
                                                                                                                                                                                                      cGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACtt                                            <  2:378564/250‑1 (MQ=255)
                                                                                                                                                                                                      cGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACtt                                            <  1:392618/250‑1 (MQ=255)
                                                                                                                                                                                                       gggATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACATCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACtttt                                          >  1:347002/1‑251 (MQ=255)
                                                                                                                                                                                                         gATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCTGAATCCACCAGTTCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGcc                                                                                                                                                                                  <  1:473680/113‑1 (MQ=255)
                                                                                                                                                                                                         gATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCt                                        <  1:413284/251‑1 (MQ=255)
                                                                                                                                                                                                                ttGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAc                                  >  1:478871/1‑250 (MQ=255)
                                                                                                                                                                                                                         aGGGCCAAACAGCAGATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTGCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACGGGCTCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAg                         <  1:414834/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aCAGCAGATCTAACGCA‑ATGCCAGAATCCACNAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                 <  2:244264/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gATCGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCATCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc           >  1:251182/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGAACTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc        <  2:440958/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCTCA‑ATGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCtt      >  1:83844/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 a‑aTGCCAGAATCCACCAATGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa  >  1:389351/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGT‑TGTGGTGGCAATCGGGATCCCTTTGTAA‑GCAGGGCCAAACAGCAGATCGAACTCA‑ATGCCGGAATCCACCAACGCTTCAGCG‑TAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAA  >  NC_000913/3815344‑3815832

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