Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,810,294 T→G 6.5% intergenic (+55/‑114) eno → / → ftsB enolase/cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,810,2940TG 6.5% 192.5 / 3.3 61intergenic (+55/‑114)eno/ftsBenolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base T (26/31);  new base G (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTT  >  NC_002947/1810150‑1810437
                                                                                                                                                  |                                                                                                                                               
cgcgTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTcc                                                                                                                                             >  3:291494/1‑149 (MQ=255)
    tttCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACg                                                                                                                                         >  1:217342/1‑149 (MQ=255)
    tttCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTTCAGTACGAACGTTCCTAca                                                                                                                                         <  5:285503/149‑2 (MQ=255)
          aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCCGTTTACCCCCGTCGTGCCGAGTTTCGCGCCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGATCGTACTGAACTTTCCTACGTACTaa                                                                                                                                   >  3:78698/1‑149 (MQ=255)
            tACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACg                                                                                                                                 >  5:33392/1‑149 (MQ=255)
             aCAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACgg                                                                                                                                >  3:83668/1‑149 (MQ=255)
                  cAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCACGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                            <  1:241338/149‑1 (MQ=255)
                      tgctgcGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGATTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTc                                                                                                                       >  6:282158/1‑149 (MQ=255)
                             cATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACg                                                                                                                                         <  6:308160/124‑1 (MQ=255)
                             cATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACg                                                                                                                                         >  5:308160/1‑124 (MQ=255)
                              aTCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGc                                                                                                               >  6:134088/1‑149 (MQ=255)
                                cGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTTCAGTACGAACGTTCCTAc                                                                                                                                          <  5:50719/120‑1 (MQ=255)
                                cGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTTCAGTACGAACGTTCCTAc                                                                                                                                          >  6:50719/1‑120 (MQ=255)
                                 gAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTg                                                                                                            <  8:132866/149‑1 (MQ=255)
                                   ccagCCACTCGGCGCCAAACCCGTCTACCGCGGTCGCGCCGAGTTCCGCGGCTAAGCAAGAGATGGTAAAAACGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGc                                                                                                          <  4:230437/147‑1 (MQ=255)
                                     agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                                        >  7:60556/1‑149 (MQ=255)
                                      gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACgg                                                                                                                                >  2:77524/1‑124 (MQ=255)
                                      gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACgg                                                                                                                                <  1:77524/124‑1 (MQ=255)
                                      gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTc                                                                                                       <  1:221685/149‑1 (MQ=255)
                                      gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTc                                                                                                       <  2:107426/149‑1 (MQ=255)
                                       cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTATGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTAc                                                                                                                                          <  1:283550/113‑1 (MQ=255)
                                       cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTAc                                                                                                                                          >  2:283550/1‑113 (MQ=255)
                                       cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCg                                                                                                      >  6:66607/1‑149 (MQ=255)
                                        aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCgg                                                                                                     >  5:56186/1‑149 (MQ=255)
                                          cTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGcc                                                                                                   <  3:56620/149‑1 (MQ=255)
                                            gggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTcc                                                                                                                                             >  6:164817/1‑105 (MQ=255)
                                            gggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTcc                                                                                                                                             <  5:164817/105‑1 (MQ=255)
                                             ggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAgg                                                                                                <  6:56186/149‑1 (MQ=255)
                                              gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCACGc                                                                                               >  4:236863/1‑149 (MQ=255)
                                               cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCt                                                                                              <  5:257830/149‑1 (MQ=255)
                                               cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                              <  4:291494/149‑1 (MQ=255)
                                                gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCtt                                                                                             <  2:217342/149‑1 (MQ=255)
                                                    aaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTg                                                                                         <  3:81836/149‑1 (MQ=255)
                                                     aGCGGTTTACCGCGGTCGTGCCGAGTTTCTCTGCTAAGCAAGTGTTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGTAGGCTGGCGTCGGCCAGGCTTCGt                                                                                          <  4:78698/147‑1 (MQ=255)
                                                     aGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGc                                                                                        <  6:33392/149‑1 (MQ=255)
                                                          tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                            <  2:21479/108‑1 (MQ=255)
                                                          tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                            >  1:21479/1‑108 (MQ=255)
                                                            tACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCg                                                                                  <  3:189860/148‑1 (MQ=255)
                                                                      gTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGaa                                                                        <  4:83668/148‑1 (MQ=255)
                                                                       tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                      <  5:66607/149‑1 (MQ=255)
                                                                         ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                                    >  7:99720/1‑149 (MQ=255)
                                                                         ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCCCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTa                                                                     <  5:134088/148‑1 (MQ=255)
                                                                          cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                                    >  3:220791/1‑148 (MQ=255)
                                                                                tCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                             >  2:194122/1‑149 (MQ=255)
                                                                                   cgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAcgcg                                                          <  4:152821/149‑1 (MQ=255)
                                                                                                 aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAc                                            >  5:274507/1‑149 (MQ=255)
                                                                                                     tAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTgg                                        <  1:194122/149‑1 (MQ=255)
                                                                                                       aaaaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGa                                       <  4:220791/148‑1 (MQ=255)
                                                                                                         aaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCCGCCAGGCTTCGTGCTATTCGAGGACCCTAAGTAACAATTACGCGGCAGACTTTTTAACCGGGAGaa                                    >  2:201021/1‑149 (MQ=255)
                                                                                                              agcagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                               <  8:60556/149‑1 (MQ=255)
                                                                                                              agcagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                               <  8:99720/149‑1 (MQ=255)
                                                                                                                 agcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc                            <  6:274507/149‑1 (MQ=255)
                                                                                                                      ggAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAg                        >  8:88996/1‑149 (MQ=255)
                                                                                                                         gTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTcc                     <  4:222575/149‑1 (MQ=255)
                                                                                                                             tAGGAAAGTTCAGTACGAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                 >  5:277532/1‑149 (MQ=255)
                                                                                                                               ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                             >  6:234296/1‑103 (MQ=255)
                                                                                                                               ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                             <  5:234296/103‑1 (MQ=255)
                                                                                                                               ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg               >  6:76988/1‑149 (MQ=255)
                                                                                                                                            cTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGt   >  5:108259/1‑149 (MQ=255)
                                                                                                                                             tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGt   >  6:135429/1‑148 (MQ=255)
                                                                                                                                               aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGtt  <  3:236863/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTT  >  NC_002947/1810150‑1810437

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