Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,601 C→T A56A (GCG→GCA)  pyrE ← orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,6010CT98.0% 162.3 / ‑5.6 51A56A (GCG→GCA) pyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (14/34);  minor base G (1/0);  total (17/34)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCG  >  NC_000913/3815362‑3815847
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
gattCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCCTGATGCTCCGCCAGCGCCCCGGCGGGAGGGGTCGCCAACGGGGTCCCCTTT‑TAAGCAGGG‑CCCAAACA‑CAAATCGGACTTAACGCCCGAATCCA‑CCAAAGCTTccccgg                                                                                                                                                                                                                                               >  2:97570/4‑245 (MQ=255)
gattCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGGGTGGGCGCCAACGGGATCCCTTTG‑TTAGGAGGG‑GC‑CAACAGCAGAACGCACTCCATTCCAGAAACCC‑CCAATGCTTCAGCgg                                                                                                                                                                                                                                               >  1:434993/4‑249 (MQ=255)
                        gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAggtc                                                                                                                                                                                                                     <  1:203251/251‑4 (MQ=255)
                        gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGGGGTCGCAATCGGGATCCCTTTG‑TAAACAGGG‑CC‑AAACAGGCGATAGGACTCAAAGCCCGAATCCA‑CCCATGCTTCAGCGTAAAAACGGACTAACaaggccagggc                                                                                                                                                                                                                      >  1:363690/1‑241 (MQ=255)
                               catcatCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGGGGTCGCAATCGGGATCCCTTTGGTAACCGGGG‑CC‑AC‑CAGCAGGTCGGACTCCATTCCAGAATCCAACCATGGCT‑CAGCGTAAAAAAGGCCCAACAC‑TGCCAGGTCGGGcccc                                                                                                                                                                                                              >  1:253718/1‑250 (MQ=255)
                               catcatCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGTGGGCGCAATCGGGGTCCCCTTG‑TAAGGCGGG‑CC‑CAACAGCCGATCGGACACCATGCCAAAATCCA‑CCAATGCTTCAACGGAAAAACGGCCTAACAG‑GGCCCGGGCGCGCCCg                                                                                                                                                                                                              >  1:101675/1‑251 (MQ=255)
                                  catcTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGt                                                                                                                                                                                                            <  1:291383/250‑1 (MQ=255)
                                   atcTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATTCTCCGCCAGCGCCACGGCGGTTGTGGTCGCAATACGGATCCCCTTG‑TAAGCAGGGGCC‑AACCA‑CCGATCGAAATCCATGCCAGAATCCA‑CCAATGCTTCAGCGGAAAAACCGCCTAACAG‑TGCCAggtccccccccggat                                                                                                                                                                                                          >  1:292961/1‑237 (MQ=255)
                                       aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGTCCTAACAG‑TGCCAGGTCGCGCCCGGTAt                                                                                                                                                                                                          <  1:373636/247‑1 (MQ=255)
                                       aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAtt                                                                                                                                                                                                         <  2:434993/248‑1 (MQ=255)
                                       aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAtt                                                                                                                                                                                                         <  1:46016/248‑1 (MQ=255)
                                       aCCAGCATTACGCGTCCTGGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCATAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAt                                                                                                                                                                                                          <  1:88163/247‑1 (MQ=255)
                                           gCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGc                                                                                                                                                                                                  <  1:282872/251‑1 (MQ=255)
                                               tACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCg                                                                                                                                                                                               <  1:243618/250‑1 (MQ=255)
                                                                    accaaccaGATTGCCGCCTTCTCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTCCGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCtt                                                                                                                                                                          <  2:296541/250‑1 (MQ=255)
                                                                    accaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAATCAGTCCGGCAGTTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGAGCCCTTTG‑TAAGCAGGG‑CC‑AAACATCAGAGCGAACTCAATGCCAGAATCCA‑CCAATGCTCCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGTCGTTGAAGAAATAGGGGCt                                                                                                                                                                           <  2:432589/249‑1 (MQ=255)
                                                                    accaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCt                                                                                                                                                                           <  1:335850/249‑1 (MQ=255)
                                                                         ccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTgcgc                                                                                                                                                                    <  1:432612/251‑1 (MQ=255)
                                                                                       gtaCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCATGGTCATTATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGAGCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGAGCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGTCTCTTGCGCCCGGATTTCAGc                                                                                                                                                        <  2:382223/247‑1 (MQ=255)
                                                                                              ggatctttgcttatgggcgtttAAAGCAGTTCGGCAGGTCAAGGTAATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCATAACAG‑TGCCAGGTCGCGACCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAAcgc                                                                                                                                                <  2:177027/230‑3 (MQ=255)
                                                                                                   ttggcctcttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGcc                                                                                                                                             <  1:447506/242‑1 (MQ=255)
                                                                                                    ttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCATGATGCTCCGCCCGCGCCACGGCGGGAGTGGGCGCCATCGGGGTCCCCTTG‑TAAGCCGGG‑CC‑AAACCACCGATCGGAATCCATGGCAGAATCCA‑CCCACGCGTCCTCGTAAAAAACGCCTAACAGCTGCC‑GGCCGCGCCCGGGATTAAACAGCCCCGCGTTGCCGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGGAAAcccacccaac                                                                                                                                         >  2:240599/1‑242 (MQ=255)
                                                                                                                ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                                                                                                                             <  1:367969/251‑1 (MQ=255)
                                                                                                                ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                                                                                                                             <  1:327145/251‑1 (MQ=255)
                                                                                                                   tAAAGCAGTACGGCCGGTCAATGTCAGGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgctt                                                                                                                           <  2:26589/250‑1 (MQ=255)
                                                                                                                   tAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCt                                                                                                                                                                           <  2:424512/202‑1 (MQ=255)
                                                                                                                   tAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgctt                                                                                                                           <  1:223825/250‑1 (MQ=255)
                                                                                                                   tAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGGCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAAAAG‑TGCCAGGGCGCGCCCGGGATTAAACAGCCCGGCGGTGAAGAAATAGGGGCt                                                                                                                                                                           >  1:424512/1‑202 (MQ=255)
                                                                                                                     aaGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgcttgct                                                                                                                        <  1:423298/251‑1 (MQ=255)
                                                                                                                               gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt                                                                                                              <  1:97570/251‑1 (MQ=255)
                                                                                                                                 aggtcaaggtcaTGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCACTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTc                                                                                                            <  1:240599/251‑1 (MQ=255)
                                                                                                                                  ggtcaaggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGCAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTc                                                                                                            >  1:436868/1‑250 (MQ=255)
                                                                                                                                               aTGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGGACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGGAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCt                                                                                              >  1:230880/1‑251 (MQ=255)
                                                                                                                                                      gCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTcg                                                                                                                                                                                                                    >  1:237594/1‑126 (MQ=255)
                                                                                                                                                      gCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTcg                                                                                                                                                                                                                    <  2:237594/126‑1 (MQ=255)
                                                                                                                                                             ccACGGCGGTAGTGGTCGCAATCGAGATCCCTTGG‑TCAGCAGGG‑CC‑AAACAGCAGATCGACCTCAATGCCAGAATCCC‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCagg                                                                                <  2:56136/251‑3 (MQ=255)
                                                                                                                                                                 ggCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCtt                                                                            <  2:422325/250‑1 (MQ=255)
                                                                                                                                                                    ggtagtggtcGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCAGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGCTATGGTTTCATGCCTTCGc                                                                         >  1:403232/5‑251 (MQ=255)
                                                                                                                                                                     gtagtggtcGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGct                                                                        <  1:322680/248‑1 (MQ=255)
                                                                                                                                                                             cgCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATc                                                                 <  1:240550/249‑1 (MQ=255)
                                                                                                                                                                             cgCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACGCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAATCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCatcatc                                                                 <  2:138404/249‑1 (MQ=255)
                                                                                                                                                                                       aTCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCt                                                       >  1:382114/1‑250 (MQ=255)
                                                                                                                                                                                                  aaGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAAGCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCACTCCTCATCTTACTTTTCTACAGACaaaa                                             <  2:292961/250‑1 (MQ=255)
                                                                                                                                                                                                      aGGG‑CC‑AAACAGCAGATCGAACTCAATGCCANAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAg                                        <  2:230880/251‑1 (MQ=255)
                                                                                                                                                                                                      aGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTAAGCGTAAAAACGGCCTACCAG‑TGCCAGGTCGCCCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCCGCGTAAACTCGCCAAACTTTAAAACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCaaaaaaa                                          >  2:365291/1‑249 (MQ=255)
                                                                                                                                                                                                                          gAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCt                      <  2:298604/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGAATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGGTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaa                  >  2:102178/1‑250 (MQ=255)
                                                                                                                                                                                                                                aaTGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                <  2:363690/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gccagAATCCA‑CCAATGCTTCAGCGTAAAAANGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTTGCGCTGATATGGGTTCATGCCTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa             <  2:204825/247‑1 (MQ=255)
                                                                                                                                                                                                                                   gccagAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc              >  2:378956/5‑250 (MQ=255)
                                                                                                                                                                                                                                               ccaATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg  >  1:169696/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
GACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCG  >  NC_000913/3815362‑3815847

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