Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,637 T→G 6.2% Y27D (TAT→GAT)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,6370TG 6.2% 206.8 / 2.9 64Y27D (TAT→GAT) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base G (2/2);  total (34/30)
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.98e-01

AGTTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTG  >  NC_002947/980492‑980783
                                                                                                                                                  |                                                                                                                                                  
aGTTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATc                                                                                                                                                 <  5:184011/149‑1 (MQ=255)
    caacaaCCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATagc                                                                                                                                               >  7:71702/1‑148 (MQ=255)
    caacaaCCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATcgcg                                                                                                                                              <  5:58568/148‑1 (MQ=255)
        aacCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTg                                                                                                                                         <  4:212900/149‑1 (MQ=255)
         acCCCAACGTGCGTGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                        <  4:100470/149‑1 (MQ=255)
              aaCGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCcgc                                                                                                                                   <  2:238378/149‑1 (MQ=255)
                  tGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAg                                                                                                                               >  5:152284/1‑149 (MQ=255)
                    cgcgGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAGTGGCCagcttgtc                                                                                                                              <  8:71702/148‑6 (MQ=255)
                             tgtgGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGAt                                                                                                                    <  5:205901/149‑1 (MQ=255)
                                      ggcgAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGcc                                                                                                           <  3:262910/149‑1 (MQ=255)
                                      ggcgAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGcc                                                                                                           >  4:304165/1‑149 (MQ=255)
                                               ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                  >  2:209370/1‑149 (MQ=255)
                                               ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                  >  2:77322/1‑149 (MQ=255)
                                               ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                  >  3:57769/1‑149 (MQ=255)
                                                  aaCGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAg                                                                                               >  5:298639/1‑149 (MQ=255)
                                                    cGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                        <  7:189271/106‑1 (MQ=255)
                                                    cGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                        >  8:189271/1‑106 (MQ=255)
                                                      tttGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTg                                                                                           <  3:304165/149‑1 (MQ=255)
                                                          aGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGagca                                                                                       <  6:152284/149‑1 (MQ=255)
                                                                  tgtgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTg                                                                               >  2:316734/1‑149 (MQ=255)
                                                                   gtgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGag                                                                                         >  4:42905/1‑138 (MQ=255)
                                                                   gtgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGag                                                                                         <  3:42905/138‑1 (MQ=255)
                                                                      ggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTgcgc                                                                            <  1:165238/148‑1 (MQ=255)
                                                                      ggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTgcgc                                                                            <  1:77322/148‑1 (MQ=255)
                                                                      ggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                           >  8:79886/1‑149 (MQ=255)
                                                                       gTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTg                                                                          >  7:178881/1‑149 (MQ=255)
                                                                            ccTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGc                                                                                              <  6:237188/124‑1 (MQ=255)
                                                                            ccTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGc                                                                                              >  5:237188/1‑124 (MQ=255)
                                                                               tGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTg                                                                                                             >  8:281855/1‑106 (MQ=255)
                                                                               tGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTg                                                                                                             <  7:281855/106‑1 (MQ=255)
                                                                                    ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCgg                                                             >  7:15324/1‑149 (MQ=255)
                                                                                     tACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGc                                                            >  3:4655/1‑149 (MQ=255)
                                                                                       cGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGa                                                                                                <  1:176246/111‑1 (MQ=255)
                                                                                       cGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGa                                                                                                >  2:176246/1‑111 (MQ=255)
                                                                                        gCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCt                                                         >  1:95069/1‑149 (MQ=255)
                                                                                          aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCg                                                                                    <  3:222279/120‑1 (MQ=255)
                                                                                          aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCg                                                                                    >  4:222279/1‑120 (MQ=255)
                                                                                           ttGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCaa                                                      <  1:209370/149‑1 (MQ=255)
                                                                                             gTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAac                                                     >  2:157623/1‑148 (MQ=255)
                                                                                                    cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCTGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTa                                             >  6:39547/1‑149 (MQ=255)
                                                                                                    cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTa                                             >  2:164855/1‑149 (MQ=255)
                                                                                                    cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTa                                             >  5:3859/1‑149 (MQ=255)
                                                                                                     tCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                            >  4:315963/1‑149 (MQ=255)
                                                                                                      ccAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacgacg                                                 >  2:317386/1‑143 (MQ=255)
                                                                                                      ccAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacgacg                                                 <  1:317386/143‑1 (MQ=255)
                                                                                                       cAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAcc                                           <  6:3859/148‑1 (MQ=255)
                                                                                                            aaGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGc                                     >  8:98271/1‑149 (MQ=255)
                                                                                                                  ccGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgcg                                 >  4:307415/1‑147 (MQ=255)
                                                                                                                           tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGcgccgc                      >  6:239894/1‑149 (MQ=255)
                                                                                                                               gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcg                  >  6:94086/1‑149 (MQ=255)
                                                                                                                                aCAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTctg                 <  3:307415/149‑3 (MQ=255)
                                                                                                                                aCAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGTctg                 <  4:57769/149‑3 (MQ=255)
                                                                                                                                    gCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGTctggtaa             <  3:315963/149‑7 (MQ=255)
                                                                                                                                      tGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                            <  2:156955/116‑1 (MQ=255)
                                                                                                                                      tGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                            >  1:156955/1‑116 (MQ=255)
                                                                                                                                        gCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTa         >  3:150684/1‑149 (MQ=255)
                                                                                                                                        gCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTa         >  3:158243/1‑149 (MQ=255)
                                                                                                                                         ccagcTCGCGATAGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACg                                      >  3:255993/5‑120 (MQ=255)
                                                                                                                                         ccagcTCGCGATAGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACg                                      <  4:255993/116‑1 (MQ=255)
                                                                                                                                           aCTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                <  1:247884/107‑1 (MQ=255)
                                                                                                                                           aCTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                >  2:247884/1‑107 (MQ=255)
                                                                                                                                               cGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgct   <  4:158243/149‑1 (MQ=255)
                                                                                                                                                gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg  <  8:15324/149‑1 (MQ=255)
                                                                                                                                                gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg  <  1:316734/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGTTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTG  >  NC_002947/980492‑980783

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