Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 562,609 C→A 7.5% S112Y (TCC→TAC)  rpsD → 30S ribosomal protein S4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947562,6090CA 7.5% 178.2 / 4.1 61S112Y (TCC→TAC) rpsD30S ribosomal protein S4
Reads supporting (aligned to +/- strand):  ref base C (34/22);  new base A (2/3);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CGTATCTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAACT  >  NC_002947/562461‑562756
                                                                                                                                                    |                                                                                                                                                   
cgtATCTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTc                                                                                                                                                     >  2:114613/1‑149 (MQ=255)
 gtATCTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTc                                                                                                                                                     <  7:90195/148‑1 (MQ=255)
 gtATCTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTc                                                                                                                                                     <  3:146470/148‑1 (MQ=255)
    tCTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGa                                                                                                                                                  <  8:304982/148‑1 (MQ=255)
     cTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGt                                                                                                                                                >  8:191028/1‑149 (MQ=255)
        cGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGGAACGAGTcga                                                                                                                                             <  5:125381/149‑3 (MQ=255)
        cGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCCTTTACCGTATTGGCTTCGGCTCGACTCGTTCCTAGTcc                                                                                                                                              >  4:214823/1‑148 (MQ=255)
        cGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACATCGTTTACCGTATGGGCTTCGGCTCGACTCGGAACGAGTcga                                                                                                                                             <  6:51753/149‑3 (MQ=255)
                     cGTCAGTTCCGCGTTTACTACCAAGCTGCTGCCTCGAAAAATGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGttt                                                                                                                                 <  3:214823/148‑1 (MQ=255)
                                      ctacCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGc                                                                                                                <  6:140566/148‑1 (MQ=255)
                                      atacCAAGCTGCTGCCGCGAAAAACGGCGCAAACGGTGAGAACCTGGTGCAACTTCTCGATGGCCGTCTGGATAACGTCGTTTACCGTATTGGCTTCGGCTCGACTCGTTc                                                                                                                                                     <  3:202158/110‑1 (MQ=255)
                                       tacCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGt                                                                                                              >  5:88064/1‑149 (MQ=255)
                                         ccAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTc                                                                                                                                                     <  3:271625/108‑1 (MQ=255)
                                         ccAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTc                                                                                                                                                     >  4:271625/1‑108 (MQ=255)
                                          cAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCaa                                                                                                           >  5:123401/1‑149 (MQ=255)
                                           aaGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAAc                                                                                                          >  1:188870/1‑149 (MQ=255)
                                            aGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACg                                                                                                         <  2:185910/149‑1 (MQ=255)
                                                 ctgcCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGt                                                                                                       >  8:52919/1‑146 (MQ=255)
                                                   gcCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGGAACGAGTcgag                                                                                                                                            <  4:300472/107‑4 (MQ=255)
                                                   gcCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGGAACGAGTcgag                                                                                                                                            >  3:300472/1‑104 (MQ=255)
                                                    ccTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACt                                                                                                 <  6:225484/149‑1 (MQ=255)
                                                      tCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTg                                                                                                >  7:55322/1‑148 (MQ=255)
                                                          aaaaaGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAAc                                                                                           <  3:309689/149‑1 (MQ=255)
                                                          aaaaaGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAAc                                                                                           <  3:314719/149‑1 (MQ=255)
                                                               ggCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTcc                                                                                      <  7:191028/149‑1 (MQ=255)
                                                                  gcAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGa                                                                                                   <  5:304778/133‑1 (MQ=255)
                                                                  gcAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGa                                                                                                   >  6:304778/1‑133 (MQ=255)
                                                                  gcAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATc                                                                                   >  4:45693/1‑149 (MQ=255)
                                                                  gcAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATc                                                                                   >  8:153045/1‑149 (MQ=255)
                                                                   cAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGAAAACATTCCATcc                                                                                  <  7:170281/149‑1 (MQ=255)
                                                                    aaCCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCt                                                                                 <  3:45693/149‑1 (MQ=255)
                                                                       cGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTAcc                                                                              <  1:114613/149‑1 (MQ=255)
                                                                              aaCCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTc                                                                        <  7:52919/148‑1 (MQ=255)
                                                                                ccTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTc                                                                     >  3:24734/1‑149 (MQ=255)
                                                                                 ctgctgCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTAGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTcc                                                                    >  2:40438/1‑149 (MQ=255)
                                                                                       cAACTGCTCGAGTGCCGTCTGGATAACGTCTTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTGGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCATACCAAGTTCGTCCGGGTGa                                                              >  4:246055/1‑149 (MQ=255)
                                                                                       cAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGa                                                              >  5:221016/1‑149 (MQ=255)
                                                                                           tGCTCGAGTGCCGTCTGGATAACGTCGATTACCGTATGGGCTTCGGCTCGACTCGTTACGAGTCCCGTAAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTTAACATTACATCCTACCAAGTTAGTACTGGTTACGTg                                                          >  1:180538/1‑149 (MQ=255)
                                                                                            gCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGCTCAGCGTCAACGGTAAGACTGGAAACATTCCCTCCGACCAACTTCGTCCGGGGGACGTgg                                                         >  6:181692/1‑149 (MQ=255)
                                                                                            gCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTACACTCGTTCCGAGTCCCGTCACCAGGTTTCGCACAAAGCGATCAACGTCAACCGTAAGACTGTTAACATTCCATCCAACCAAATTAGTCCTGGAGACGAgg                                                         >  2:294230/1‑149 (MQ=255)
                                                                                              tCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTc                                                       >  3:18157/1‑149 (MQ=255)
                                                                                                     ccGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTAGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAAGATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTcc                                                <  6:88064/149‑1 (MQ=255)
                                                                                                          tGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCcgc                                              >  3:158788/1‑146 (MQ=255)
                                                                                                            gATAACGTCGTTTACCGTATGGGATTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCga                                            <  2:188870/146‑1 (MQ=255)
                                                                                                                  gtcgtTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCt                                   >  1:308792/1‑149 (MQ=255)
                                                                                                                  gtcgtTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCt                                   <  7:153045/149‑1 (MQ=255)
                                                                                                                  gtcgtTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCt                                   >  2:278366/1‑149 (MQ=255)
                                                                                                                      tttACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGc                               >  2:62600/1‑149 (MQ=255)
                                                                                                                           cGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTTCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGtt                          <  1:62600/149‑1 (MQ=255)
                                                                                                                             tATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTgc                        >  5:97545/1‑149 (MQ=255)
                                                                                                                                 ggCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTggg                                <  2:155895/137‑1 (MQ=255)
                                                                                                                                 ggCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTggg                                >  1:155895/1‑137 (MQ=255)
                                                                                                                                 ggCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCAt                    <  6:221016/149‑1 (MQ=255)
                                                                                                                                  gCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCTACGTCAACGGTAAGACTGTTAACATTCCATGCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCAtt                   >  1:85190/1‑149 (MQ=255)
                                                                                                                                      cGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGtt                >  4:187675/1‑148 (MQ=255)
                                                                                                                                      cGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGtt                <  6:123401/148‑1 (MQ=255)
                                                                                                                                        gCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCGGTTGCGCATTGTTCaa             <  8:55322/149‑1 (MQ=255)
                                                                                                                                        gCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTc               >  1:15013/1‑147 (MQ=255)
                                                                                                                                              cTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTTTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCtt       >  7:186336/1‑149 (MQ=255)
                                                                                                                                               tCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAGGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGccc         <  7:986/146‑1 (MQ=255)
                                                                                                                                               tCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAGGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGccc         >  8:986/1‑146 (MQ=255)
                                                                                                                                               tCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCtt       >  2:285566/1‑148 (MQ=255)
                                                                                                                                                  ttCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAAc   >  8:275318/1‑149 (MQ=255)
                                                                                                                                                  ttCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAAc   >  1:188380/1‑149 (MQ=255)
                                                                                                                                                  ttCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATCCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAAc   >  8:14033/1‑149 (MQ=255)
                                                                                                                                                   tCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAACt  >  8:207009/1‑149 (MQ=255)
                                                                                                                                                   tCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAACt  >  4:135809/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTATCTACGGTGTTCTGGAGCGTCAGTTCCGCGGTTACTACCAAGCTGCTGCCTCGAAAAAAGGCGCAACCGGTGAGAACCTGCTGCAACTGCTCGAGTGCCGTCTGGATAACGTCGTTTACCGTATGGGCTTCGGCTCGACTCGTTCCGAGTCCCGTCAGCTGGTTTCGCACAAAGCGATCAGCGTCAACGGTAAGACTGTAAACATTCCATCCTACCAAGTTCGTCCGGGTGACGTGGTCGCAGTCCGCGAGAAATCGCTGGGCCAGTTGCGCATTGTTCAAGCCCTTGAACT  >  NC_002947/562461‑562756

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: ‑