Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,831 C→G 7.5% intergenic (+30/+109) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,8310CG 7.5% 289.0 / 14.8 93intergenic (+30/+109)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (49/35);  new base G (5/2);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.72e-01

TACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGC  >  NC_000913/2282685‑2282977
                                                                                                                                                  |                                                                                                                                                  
tACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg                                                                                                                                                  >  2:134012/1‑149 (MQ=255)
   aGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcggcg                                                                                                                                               >  7:16130/1‑149 (MQ=255)
   aGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcggcg                                                                                                                                               >  2:16314/1‑149 (MQ=255)
    gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcggcgg                                                                                                                                              <  2:159794/149‑1 (MQ=255)
     cGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCggg                                                                                                                                             <  8:198310/149‑1 (MQ=255)
          ttCTGACCATAGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACCTACAGCGTCACCCGTACACACCGCAGCTGCAGCACATCGACTTCGTTCGCGCTTCATTGCTGAATAAGTTGTAAATAACCCCGCTCCGGCGGGGtttt                                                                                                                                        <  4:231704/149‑1 (MQ=255)
            cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTACGGCGGGGtttttt                                                                                                                                      >  3:167548/1‑149 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCTGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGt                                                                                                                                    <  6:313340/148‑1 (MQ=255)
                 cATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTa                                                                                                                                   <  3:113286/147‑1 (MQ=255)
                 cATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTa                                                                                                                                   >  4:113286/1‑147 (MQ=255)
                 cATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATc                                                                                                                                 >  7:139322/1‑149 (MQ=255)
                 cATCGTTGTTGACGGGAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATc                                                                                                                                 <  6:301571/149‑1 (MQ=255)
                  aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCt                                                                                                                                >  7:130074/1‑149 (MQ=255)
                   tCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCGCAGGACGTACAGCGTCACCCGTACAACCCGAAGCTGCAGCACATCGACTTCGGTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTg                                                                                                                               >  4:49107/1‑149 (MQ=255)
                      ttgttgACGGGAAAGAAATCAAAGTTAAAGCGCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCGCCGGCGGGGTTTTTTGTATCTGCAg                                                                                                                            <  1:30422/149‑1 (MQ=255)
                           gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGttttt                                                                                                                                       >  5:309620/1‑133 (MQ=255)
                           gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGttttt                                                                                                                                       <  6:309620/133‑1 (MQ=255)
                           gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTAt                                                                                                                       >  3:155003/1‑149 (MQ=255)
                           gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACATCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTAt                                                                                                                       >  1:293730/1‑149 (MQ=255)
                              ggTAAAGAAATCAAAGTTAAAGGTCAGGACTTACAGCGTAACCCTTACAAACCGTAGCTGCAGCACATCGACTTCGATCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCGTCTACGGCGGGGTTTTTTGTATCTGCAGATTATGcc                                                                                                                    <  2:293730/149‑1 (MQ=255)
                              ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGAtt                                                                                                                         >  5:15627/1‑144 (MQ=255)
                              ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGAtt                                                                                                                         <  6:15627/144‑1 (MQ=255)
                              ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGGTTTTTGTATCTGCAGATTATGcc                                                                                                                    >  2:130029/1‑149 (MQ=255)
                              ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCACCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGGTCCGGCGGGGGTTTTTGTATCTGCTGATTATggc                                                                                                                    >  2:123425/1‑147 (MQ=255)
                               gTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCTGCGGGGTTTTTTGTATCTGCAGATTATGCCt                                                                                                                   <  8:139322/149‑1 (MQ=255)
                               gTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCt                                                                                                                   >  7:107725/1‑149 (MQ=255)
                                     aaaTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTAc                                                                                                             <  4:167548/149‑1 (MQ=255)
                                          aaaGTTAAAGCTCAGGACGTACAGCGTCACCCGTAAAAACCGAAGCTGCAGCACATCGCCTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGAGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTa                                                                                                         <  3:66911/148‑1 (MQ=255)
                                            aGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTg                                                                                                      >  5:113090/1‑149 (MQ=255)
                                              ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGc                                                                                                     >  3:145344/1‑148 (MQ=255)
                                                    cTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGtttttt                                                                                                                                      >  4:304835/1‑109 (MQ=255)
                                                    cTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGtttttt                                                                                                                                      <  3:304835/109‑1 (MQ=255)
                                                     tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                             <  7:324766/149‑1 (MQ=255)
                                                     tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                             >  7:342571/1‑149 (MQ=255)
                                                       aGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGc                                                                                           <  8:107725/149‑1 (MQ=255)
                                                        ggACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTg                                                                                            >  2:316238/1‑147 (MQ=255)
                                                         gACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                         <  6:113090/149‑1 (MQ=255)
                                                           cGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGc                                                                                       <  4:315500/149‑1 (MQ=255)
                                                            gTACAGCGTCACCCGTACCAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGcc                                                                                      <  2:328652/149‑1 (MQ=255)
                                                             tACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGc                                                                                       >  7:238063/1‑147 (MQ=255)
                                                                        ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                          >  6:60738/1‑149 (MQ=255)
                                                                        ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTGTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                          >  5:54190/1‑149 (MQ=255)
                                                                             cAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCa                                                                                     >  2:105553/1‑133 (MQ=255)
                                                                             cAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCa                                                                                     <  1:105553/133‑1 (MQ=255)
                                                                              aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                         >  7:49465/1‑128 (MQ=255)
                                                                              aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                         <  8:49465/128‑1 (MQ=255)
                                                                               caCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACACCGCTCCGGCGGGGTTATTTGAATATGCAGAATATGCCTAATTACGATAATGCTATATTTTGCAGGCCAGAAAAGTCGTTCAcgccg                                                                   >  2:248874/2‑149 (MQ=255)
                                                                                         gcagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAt                                                         >  4:241030/1‑149 (MQ=255)
                                                                                          cagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAtt                                                        <  1:108404/149‑1 (MQ=255)
                                                                                          cagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAtt                                                        >  5:30491/1‑149 (MQ=255)
                                                                                          cagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAt                                                         <  8:16130/148‑1 (MQ=255)
                                                                                            gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                      <  1:134012/149‑1 (MQ=255)
                                                                                             cacaTCGTCTTCGTTCGCGCATAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGCTATTGCTATTTTTAGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAg                                                     <  1:123425/149‑1 (MQ=255)
                                                                                             cacaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAtt                                                        >  6:79665/1‑146 (MQ=255)
                                                                                               caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                   <  1:316238/149‑1 (MQ=255)
                                                                                                 tCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAaca                                                 <  6:30491/149‑1 (MQ=255)
                                                                                                    aCTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAAcaca                                               >  8:27940/1‑148 (MQ=255)
                                                                                                         gttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCg                                         <  5:60738/149‑1 (MQ=255)
                                                                                                          ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTCTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                        >  2:68775/1‑149 (MQ=255)
                                                                                                              cgcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                    <  4:88279/149‑1 (MQ=255)
                                                                                                              cgcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                    >  6:249040/1‑149 (MQ=255)
                                                                                                               gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTTTATCTGCAGATTATGCCTGATTACGTTATTGCTATTTTTTGAAGGCAAGATAAGGCGTTAACGCCGCATATGGCATTCAGAACACAAAACCGATTAAt                                   >  3:257232/1‑149 (MQ=255)
                                                                                                               gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                   >  3:296604/1‑149 (MQ=255)
                                                                                                                cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                  <  3:241030/149‑1 (MQ=255)
                                                                                                                cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                  >  2:96914/1‑149 (MQ=255)
                                                                                                                 ttAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCTTCTGGCATGAAGACCACAAAGCCTATTAATTg                                 >  6:237045/1‑149 (MQ=255)
                                                                                                                        cTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAg                          <  8:342571/149‑1 (MQ=255)
                                                                                                                           aaTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGGAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                       >  7:182856/1‑149 (MQ=255)
                                                                                                                           aaTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                       >  1:134792/1‑149 (MQ=255)
                                                                                                                                gatacaaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                        <  8:124246/122‑1 (MQ=38)
                                                                                                                                gatacaaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                        >  7:124246/6‑127 (MQ=38)
                                                                                                                                gTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                   <  4:281361/132‑1 (MQ=255)
                                                                                                                                gTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                   >  3:281361/1‑132 (MQ=255)
                                                                                                                                gTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGAGTATGGCTGATTACGGTATTGGTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTTATTGCCGCCAGCCGTCCgg                  >  6:228516/1‑149 (MQ=39)
                                                                                                                                    tAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt              >  1:191267/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg             >  3:188798/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg             >  2:233097/1‑149 (MQ=255)
                                                                                                                                      aaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                  <  2:297102/127‑1 (MQ=255)
                                                                                                                                      aaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                  >  1:297102/1‑127 (MQ=255)
                                                                                                                                      aaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg             >  6:143328/1‑148 (MQ=255)
                                                                                                                                       aaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc            >  2:345895/1‑148 (MQ=255)
                                                                                                                                       aaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGAGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc            >  4:176418/1‑148 (MQ=255)
                                                                                                                                        aaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGACTGATTACTGTATTGCTATTTTTTGAAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCTCTGaaa          >  4:39114/1‑149 (MQ=39)
                                                                                                                                         aaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAt         >  5:73480/1‑149 (MQ=255)
                                                                                                                                         aaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAt         <  3:49107/149‑1 (MQ=255)
                                                                                                                                           ccccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTcc                    >  7:316532/1‑136 (MQ=255)
                                                                                                                                           ccccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTcc                    <  8:316532/136‑1 (MQ=255)
                                                                                                                                           ccccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTg       >  8:148551/1‑149 (MQ=255)
                                                                                                                                            cccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg      >  6:298280/1‑149 (MQ=255)
                                                                                                                                            cccGCTCCGGCGGGGTTTTTTGTATATGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGAAGCATCTGGCATTCAGAAAACAAAACAGATTAATTGCCGACAGACGTCATGCGCTGCAATTgg      >  2:108204/1‑149 (MQ=255)
                                                                                                                                             ccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt     >  7:302499/1‑149 (MQ=255)
                                                                                                                                              cGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTc    <  6:73480/149‑1 (MQ=255)
                                                                                                                                                cTCCGGCGGGGTTTTTTGTATCTGCAGTTTTTGCCTGATTTCGGTTTTGCTATTTTTTGGAGGCCAGTTAAGGCGTTCACGCCGCTTCTGGCATTCAGAACAGAAATGCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTcgg  <  4:151946/149‑2 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGC  >  NC_000913/2282685‑2282977

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