Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,145 C→G 13.5% intergenic (+128/‑70) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1450CG13.5% 1023.4 / 53.1 342intergenic (+128/‑70)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (224/70);  new base G (3/45);  total (227/115)
Fisher's exact test for biased strand distribution p-value = 3.69e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

CGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCT  >  NC_000913/3780998‑3781293
                                                                                                                                                     |                                                                                                                                                     
cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACACTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCGCGCCCCTCCGGGGAGAGGGCCg                                                                                                                                                      >  4:535682/1‑149 (MQ=255)
 gAAAGGGAATGCGGGATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCCCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     >  3:63874/1‑149 (MQ=255)
 gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     >  2:566856/1‑149 (MQ=255)
 gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     <  4:246829/149‑1 (MQ=255)
  aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                                                                    <  2:474916/149‑1 (MQ=255)
  aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                                                                    <  1:272843/149‑1 (MQ=255)
  aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTCCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                                                                    <  3:59050/149‑1 (MQ=255)
  aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACGCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                                                                    <  4:696780/149‑1 (MQ=255)
     gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                                                 >  3:131762/1‑149 (MQ=255)
       gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTggg                                                                                                                                               >  3:693469/1‑147 (MQ=255)
       gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                                                               <  1:285012/149‑1 (MQ=255)
        aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                                              <  1:122775/149‑1 (MQ=255)
        aaTGCGGCATAGCCGTTTGCGCCCCCTCCCCCTAAACCTCCCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                                              <  1:173278/149‑1 (MQ=255)
              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                                                                       >  2:539646/1‑149 (MQ=255)
              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                                                                       <  2:231820/149‑1 (MQ=255)
                aTAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                                              >  4:570619/1‑141 (MQ=255)
                aTAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGAGAGGGGAAAAgg                                                                                                                                     >  4:784220/1‑149 (MQ=255)
                aTAGTCGTTTGCCCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGCTCGGCGCTGTATGTACTCCCTCACTCTGACACGACACCGCACTCTTTTTTTCCCCATCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                                                                     <  2:439820/149‑1 (MQ=255)
                 tAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                                                                                    >  2:553322/1‑149 (MQ=255)
                   gTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGGGGGGAAAAGGCCg                                                                                                                                  >  2:304301/1‑149 (MQ=255)
                   gTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                                                                                                                  >  1:583/1‑149 (MQ=255)
                    tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                                                             <  1:803293/138‑1 (MQ=255)
                    tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                                                                                                                 >  2:608709/1‑149 (MQ=255)
                    tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                                                                                                                 >  2:362676/1‑149 (MQ=255)
                    tCGTTTGCGCCCCCTCACCATAACCCTCTCCCTCAGGGAGAGGGGACCGCTTGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGACATGGCCGc                                                                                                                                 >  2:839139/1‑149 (MQ=255)
                      gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                                                                               >  3:435876/1‑149 (MQ=255)
                         tGCGCCCCCTCACCCAACCCCCCCCCCTCGGGGGGGGGGGCCCGTCCGGTGCTTTATCTCCCCCCCCCCTCTAAAACGACCCCGCTCTTTTTTTTCCCCCCCCGCCCCCCGGGGGAGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGc                                                                                                                            <  4:261717/149‑1 (MQ=255)
                                 cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                                                                    >  4:679143/1‑149 (MQ=255)
                                 cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                                                                    >  4:320014/1‑149 (MQ=255)
                                  tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                                                                                   >  4:722819/1‑149 (MQ=255)
                                      ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                               >  4:383235/1‑149 (MQ=255)
                                      ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                               >  3:501673/1‑149 (MQ=255)
                                      ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                               >  1:719137/1‑149 (MQ=255)
                                      ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                               >  3:21934/1‑149 (MQ=255)
                                       cTACCCCTCCCCCCCGGGAGGGGGGGACGTTCGGGCCCTTTTTGTCCTCCCCCCCCCGGAAAAGACCCCGCCTTTTTTTTTCCCCCCTCCCCCCCCGGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                                                              <  1:890256/149‑1 (MQ=255)
                                       cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                                                              >  2:890358/1‑149 (MQ=255)
                                        tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                                                             >  3:875111/1‑149 (MQ=255)
                                         aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                                            >  4:60586/1‑149 (MQ=255)
                                           cccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                                                             <  3:326842/115‑1 (MQ=255)
                                            ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  3:532464/1‑149 (MQ=255)
                                            ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  1:610344/1‑149 (MQ=255)
                                            ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  1:523526/1‑149 (MQ=255)
                                            ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                                                 >  1:220358/1‑110 (MQ=255)
                                            ccTCTCCCTCAGGGAGAGGAGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  3:729919/1‑149 (MQ=255)
                                             ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                                                                        >  1:642549/1‑149 (MQ=255)
                                             ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGATAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                                                                        >  4:776344/1‑149 (MQ=255)
                                              tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                                                                       >  1:824222/1‑149 (MQ=255)
                                               cccccTCGGGGAAGGGGGCCCTTTCGGCGTTTTATTTCCCCCCCCCCTCGAAACCGCCACCGCATTTTTTTTTTCTCCCTCGCCCCCCGGGGGAGAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                                                                      <  4:710156/147‑1 (MQ=255)
                                                  ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTATGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                   >  2:162260/1‑149 (MQ=255)
                                                   cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                                                                  >  1:185480/1‑149 (MQ=255)
                                                    tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                                                                 >  1:775776/1‑149 (MQ=255)
                                                      aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                                                               >  3:252545/1‑149 (MQ=255)
                                                            gggggCCCGTCCGGCGTTGTATTTCCCCCCTCCCTCGAAACCGCCCCCCCCTTTTTTTTTTCCCCCCCGCCCCTGGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                                                                         <  1:65657/148‑1 (MQ=255)
                                                              gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                                                               <  2:275558/94‑1 (MQ=255)
                                                                gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                                                                     >  3:361426/1‑149 (MQ=255)
                                                                gACCGTTCGGCGCTGTATCTACTCCCTCACTCTGAAACGAGACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                                                                     >  1:260917/1‑149 (MQ=255)
                                                                  ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAGAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                                                                                   <  1:228100/149‑1 (MQ=255)
                                                                   cgctCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTCTCCCTCGCCCCCCCGGGGGGAGGGGCCGGGGGAGGGGAAAAGGGCGCCCTTGTTCCCAATTTTCTTCTTTTTTGCCCCCGGCCTAAggggggcgc                                                                                 >  3:417601/4‑145 (MQ=255)
                                                                            gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                                                                          >  4:418079/1‑149 (MQ=255)
                                                                             ttgtttgtcCCCCCTCCCTCTGAAACGCCACCCCATTTTTTTTTTCCCCCCCCCCCCGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                                                                         <  1:30288/148‑1 (MQ=255)
                                                                                 ttgtaCCCCCCCCCTCTGAAACGCCCCCCCTCTTTTTTTTCCCCCCTCGCCCCCCGGGGGGGAGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                                                                     <  1:401464/148‑1 (MQ=255)
                                                                                     aCTCCCTCACTCTGAAACGCCCCCGCTCTTTTTTTTTCCCCCTCGCCCCCCGGGGGGGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                                                                 <  4:175011/149‑1 (MQ=255)
                                                                                     aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                                                                 >  1:784777/1‑149 (MQ=255)
                                                                                      cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTTCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                                                                >  1:365149/1‑149 (MQ=255)
                                                                                       cccccCACTCGGAAACGACCCCGCACTTTTTTTTCCCCCCCCGCCCCTCGGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt                                                               <  4:76105/148‑1 (MQ=255)
                                                                                         ctccACTCTGAAACACCCCCAAACTTTTTTTTCTTCCCCGCCCCCCCCGGGGGGGGGGTCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                                                             <  2:695543/146‑1 (MQ=255)
                                                                                         cctccctcAAAAACGCCCCCCCTCTTTTTTTTCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                                                             <  3:409705/149‑1 (MQ=255)
                                                                                          cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTTAGTGTGCAGTATGATAAACATCGTACTTTAc                                                            >  3:301939/1‑149 (MQ=255)
                                                                                           tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg                                                           >  3:154283/1‑149 (MQ=255)
                                                                                            cACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGTTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa                                                          >  2:196609/1‑149 (MQ=255)
                                                                                             aCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa                                                         >  3:859178/1‑149 (MQ=255)
                                                                                              cgctGAAACGACCCCGCCCTTTTTTTTTCTCCCTCCCCCCCCGGGGGAGAGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc                                                        <  2:780665/147‑1 (MQ=255)
                                                                                                 tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg                                                     >  4:295018/1‑149 (MQ=255)
                                                                                                    aaCGACACCGCACTCTTTTTTTCTCCATCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaaa                                                  <  3:814293/149‑1 (MQ=255)
                                                                                                              cACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaaa                                        >  2:648094/1‑149 (MQ=255)
                                                                                                                 tttttttttCCCCCCCCCCCCCCGGGGGGGAGGGCGGGGGTGAGGGGAAAGGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATAc                                     <  2:595482/147‑1 (MQ=255)
                                                                                                                   tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTg                                   >  3:46349/1‑149 (MQ=255)
                                                                                                                   ttttttgCCCCCTCGCCCCGCGGGGGGGGGGCCGGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTg                                   <  4:154283/149‑1 (MQ=255)
                                                                                                                    tttttcCTCCCCCCCCCCCCGGGGGGGAGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTgg                                  <  4:343803/149‑1 (MQ=255)
                                                                                                                    tttttcCCCCCCCCCCCCCCGGGGGGGGGGGCCGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTgg                                  <  1:435790/149‑1 (MQ=255)
                                                                                                                       cccggcccTCTCCCCGGAGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCaa                               >  4:457441/6‑149 (MQ=255)
                                                                                                                           tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTTCGCCAAATACTGGCAAcatc                           >  2:393629/1‑149 (MQ=255)
                                                                                                                           tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatc                           >  3:790719/1‑149 (MQ=255)
                                                                                                                            cccTCCCCCCCCCGGGGAGGGGGCCGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatca                          <  4:131762/149‑1 (MQ=255)
                                                                                                                             ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                                                                       >  1:570137/1‑39 (MQ=255)
                                                                                                                              cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACatcatc                        >  1:201581/1‑149 (MQ=255)
                                                                                                                              cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACatcatc                        >  2:723977/1‑149 (MQ=255)
                                                                                                                              cTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACatcatc                        >  2:91843/1‑149 (MQ=255)
                                                                                                                               tCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc                       >  2:720082/1‑149 (MQ=255)
                                                                                                                                    cccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCt                   >  1:725783/1‑149 (MQ=255)
                                                                                                                                      ttcgggggggggggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCttt                 <  2:765200/145‑1 (MQ=255)
                                                                                                                                      cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCttt                 >  4:590923/1‑149 (MQ=255)
                                                                                                                                       tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTg                >  1:555225/1‑149 (MQ=255)
                                                                                                                                          ggggggggggCCGGGGGAGGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcgc             <  2:256766/149‑1 (MQ=255)
                                                                                                                                             gagaGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCa                                                      <  4:682283/105‑1 (MQ=255)
                                                                                                                                             gagaGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCa                                                      >  3:682283/1‑105 (MQ=255)
                                                                                                                                             gagaGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTaa          <  4:427947/149‑1 (MQ=255)
                                                                                                                                                ggggCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATAc       <  4:4844/148‑1 (MQ=255)
                                                                                                                                                    ccGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGc   <  2:582601/149‑1 (MQ=255)
                                                                                                                                                     cGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCt  >  1:456206/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCT  >  NC_000913/3780998‑3781293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: