Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,761 C→A 6.9% P68Q (CCG→CAG)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,7610CA 6.9% 271.5 / 7.0 87P68Q (CCG→CAG) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (35/46);  new base A (4/2);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

GATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGC  >  NC_002947/980613‑980909
                                                                                                                                                    |                                                                                                                                                     
gATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTCcgcc                                                                                                                                                       >  6:114161/1‑149 (MQ=255)
  tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGcgccgc                                                                                                                                                     >  2:4695/1‑149 (MQ=255)
    tcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGt                                                                                                                                                   >  5:137694/1‑149 (MQ=255)
      gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGcgccg                                                                                                                                                      <  8:117732/144‑1 (MQ=255)
      gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTc                                                                                                                                                  <  8:165710/148‑1 (MQ=255)
      gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGTc                                                                                                                                                  <  2:91936/148‑1 (MQ=255)
      gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACCCTGCGCAGTGCGCCGCGTcg                                                                                                                                                 >  6:102097/1‑149 (MQ=255)
         aaGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgc                                                                                                                                              <  8:97776/149‑1 (MQ=255)
           gCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                                                                                                                                            >  8:67491/1‑149 (MQ=255)
           gCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                                                                                                                                            >  2:214800/1‑149 (MQ=255)
           gCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCAACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                                                                                                                                            >  1:176114/1‑149 (MQ=255)
            cTGGCCACTCGCTATCTCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgcg                                                                                                                                           <  2:299917/149‑1 (MQ=255)
            cTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgcg                                                                                                                                           <  4:299086/149‑1 (MQ=255)
                 cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGTGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                                                                                                                                            >  7:286490/1‑143 (MQ=255)
                 cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGTGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                                                                                                                                            <  8:286490/143‑1 (MQ=255)
                      gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                                 <  4:88799/149‑1 (MQ=255)
                      gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                                 <  4:30461/149‑1 (MQ=255)
                      gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                                 <  8:332810/149‑1 (MQ=255)
                      gCTATCGCGAACTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgc                                                                                                                                   <  3:107515/147‑1 (MQ=255)
                         ttCGCGGGCCAGCTCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGc                                                                                                                               >  6:185419/2‑148 (MQ=255)
                            gcgAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt                                                                                                                            <  1:303765/148‑1 (MQ=255)
                             cgAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCa                                                                                                                          <  3:214674/149‑1 (MQ=255)
                                gCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCg                                                                                                                       <  2:37278/149‑1 (MQ=255)
                                gCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCg                                                                                                                       <  1:214800/149‑1 (MQ=255)
                                      ccGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCAc                                                                                                                  >  8:165750/1‑148 (MQ=255)
                                      ccGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCAc                                                                                                                  >  3:244161/1‑148 (MQ=255)
                                        gcgAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGaa                                                                                                               <  7:165750/149‑1 (MQ=255)
                                         cgAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAg                                                                                                              >  7:88096/1‑149 (MQ=255)
                                            aGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCTCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTTCCTGCTGGCCATCAGCGGCCACGAAGTg                                                                                                            <  1:23109/148‑1 (MQ=255)
                                            aGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGc                                                                                                           <  2:350026/149‑1 (MQ=255)
                                              tCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgc                                                                                                                                   >  8:311438/1‑123 (MQ=255)
                                              tCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgc                                                                                                                                   <  7:311438/123‑1 (MQ=255)
                                                 aTCCTGATCGCTTTGCTGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGt                                                                                                                                                   >  7:129743/1‑104 (MQ=255)
                                                 aTCCTGATCGCTTTGCTGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGt                                                                                                                                                   <  8:129743/104‑1 (MQ=255)
                                                 aTCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCag                                                                                                      >  2:162369/1‑149 (MQ=255)
                                                    cTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCaggag                                                                                                   <  8:88096/149‑1 (MQ=255)
                                                      gATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCaggagg                                                                                                  >  8:252751/1‑148 (MQ=255)
                                                           cTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAAGTCAggg                                                                                            >  3:42706/1‑149 (MQ=255)
                                                            tttGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGt                                                                                           >  6:383264/1‑149 (MQ=255)
                                                              tGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACg                                                                                             <  5:116706/145‑1 (MQ=255)
                                                              tGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACg                                                                                             >  6:116706/1‑145 (MQ=255)
                                                              tGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTc                                                                                          <  3:251708/148‑1 (MQ=255)
                                                                ccGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                                       <  1:227475/149‑1 (MQ=255)
                                                                ccGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                                       >  6:116281/1‑149 (MQ=255)
                                                                 cGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                                       >  4:228249/1‑148 (MQ=255)
                                                                 cGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                                       >  3:355726/1‑148 (MQ=255)
                                                                   aCGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑c                                                                                   <  2:144156/149‑1 (MQ=255)
                                                                       cTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCt                                                                                >  4:276305/1‑148 (MQ=255)
                                                                       cTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCt                                                                                >  2:122152/1‑148 (MQ=255)
                                                                        tCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGa                                                                              >  7:243324/1‑149 (MQ=255)
                                                                        tCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGa                                                                              >  4:119344/1‑149 (MQ=255)
                                                                          cGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACCCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGa                                                                              <  5:102097/147‑1 (MQ=255)
                                                                           gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑ccc                                                                                 >  3:236495/1‑143 (MQ=255)
                                                                           gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑ccc                                                                                 <  4:236495/143‑1 (MQ=255)
                                                                             gCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTcc                                                                         <  6:2823/149‑1 (MQ=255)
                                                                                tGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAAGTCAGGGTCGTGGA‑CCGTGACCTCctgg                                                                      <  4:42706/149‑2 (MQ=255)
                                                                                  agcagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAAGTCAGGGTCGTGGA‑CCGTGACCTCctggg                                                                     <  5:359905/148‑3 (MQ=38)
                                                                                              cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCCTGACTTCCTGCTGCAGCAGATGCa                                                         <  1:4695/149‑1 (MQ=255)
                                                                                              cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAg                                                        <  3:57819/149‑1 (MQ=255)
                                                                                                 tGGAATAATCCTCGGACCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTc                                                                                          <  7:130060/113‑1 (MQ=255)
                                                                                                   gAAAAATCCGGGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTgg                                                     <  3:184993/147‑1 (MQ=255)
                                                                                                   gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                                    <  5:114161/148‑1 (MQ=255)
                                                                                                   gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                                    <  3:276305/148‑1 (MQ=255)
                                                                                                       aaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCTTCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                               <  1:122152/149‑1 (MQ=255)
                                                                                                       aaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                               <  4:244161/149‑1 (MQ=255)
                                                                                                         tCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                               >  6:46354/1‑147 (MQ=255)
                                                                                                         tCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                               >  3:57166/1‑147 (MQ=255)
                                                                                                          ccGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                            >  5:324424/1‑149 (MQ=255)
                                                                                                               gCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCg                                        >  5:207271/1‑148 (MQ=255)
                                                                                                                cccTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaa                                      <  5:46354/149‑1 (MQ=255)
                                                                                                                 ccTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCg                                        >  2:253259/1‑146 (MQ=255)
                                                                                                                   tCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                                   >  7:80942/1‑149 (MQ=255)
                                                                                                                   tCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                                   >  4:305267/1‑149 (MQ=255)
                                                                                                                     aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                                   <  7:252751/147‑1 (MQ=255)
                                                                                                                     aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                                   <  2:176114/147‑1 (MQ=255)
                                                                                                                      acgacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                                <  1:253259/149‑1 (MQ=255)
                                                                                                                      acgacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                                >  4:57419/1‑149 (MQ=255)
                                                                                                                         acgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAg                             <  3:228249/149‑1 (MQ=255)
                                                                                                                           gCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGa                           <  5:185419/149‑1 (MQ=255)
                                                                                                                             ttACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGAt                          >  5:232341/1‑148 (MQ=255)
                                                                                                                               gcgcgACGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATg                         >  6:129260/5‑147 (MQ=255)
                                                                                                                                ccAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg                      <  7:67491/149‑1 (MQ=255)
                                                                                                                                 cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGc                     >  3:149916/1‑149 (MQ=255)
                                                                                                                                    aCGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa                  <  6:137694/149‑1 (MQ=255)
                                                                                                                                    aCGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa                  >  1:379586/1‑149 (MQ=255)
                                                                                                                                        tGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCtt              >  4:199825/1‑149 (MQ=255)
                                                                                                                                         gcgcACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTg             >  5:215368/1‑149 (MQ=255)
                                                                                                                                          cgcAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTg             <  6:207271/148‑1 (MQ=255)
                                                                                                                                           gcAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGAc           <  3:119344/149‑1 (MQ=255)
                                                                                                                                            cAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACg          <  3:305267/149‑1 (MQ=255)
                                                                                                                                                 cgccgcGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                                >  8:101274/1‑122 (MQ=255)
                                                                                                                                                 cgccgcGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                                <  7:101274/122‑1 (MQ=255)
                                                                                                                                                    cgcgTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGc  >  3:201305/1‑149 (MQ=255)
                                                                                                                                                    cgcgTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGc  <  6:232341/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGA‑CCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGC  >  NC_002947/980613‑980909

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