Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,830 T→C 7.6% intergenic (+29/+110) 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,8300TC 7.6% 272.1 / 10.7 93intergenic (+29/+110)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base T (50/35);  new base C (5/2);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

TTCTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGC  >  NC_000913/2282682‑2282977
                                                                                                                                                    |                                                                                                                                                   
ttCTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAACGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCt                                                                                                                                                     >  2:108404/1‑149 (MQ=255)
   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                                                                                                                                   >  4:113286/1‑147 (MQ=255)
                    cATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTa                                                                                                                                   <  3:113286/147‑1 (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                                                                                                                                       <  6:309620/133‑1 (MQ=255)
                              gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGttttt                                                                                                                                       >  5:309620/1‑133 (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:342571/1‑149 (MQ=255)
                                                        tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                             <  7:324766/149‑1 (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                                                                                         <  8:49465/128‑1 (MQ=255)
                                                                                 aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                         >  7:49465/1‑128 (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)
                                                                                                                                                    |                                                                                                                                                   
TTCTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGC  >  NC_000913/2282682‑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: