Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,755 G→C 6.5% S66T (AGT→ACT)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,7550GC 6.5% 296.2 / 8.9 92S66T (AGT→ACT) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base G (39/47);  new base C (4/2);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

GGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACG  >  NC_002947/980612‑980901
                                                                                                                                               |                                                                                                                                                   
ggATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTgcgc                                                                                                                                                >  6:359905/1‑149 (MQ=255)
ggATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTgcgc                                                                                                                                                >  3:299086/1‑149 (MQ=255)
 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                                                                                                                                    <  8:286490/143‑1 (MQ=255)
                  cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGTGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                                                                                                                                    >  7:286490/1‑143 (MQ=255)
                       gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                         <  4:30461/149‑1 (MQ=255)
                       gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                         <  4:88799/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                                                                                                               <  1:214800/149‑1 (MQ=255)
                                 gCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCg                                                                                                               <  2:37278/149‑1 (MQ=255)
                                       ccGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCAc                                                                                                          >  3:244161/1‑148 (MQ=255)
                                       ccGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCAc                                                                                                          >  8:165750/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                                                                                                                           <  7:311438/123‑1 (MQ=255)
                                               tCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgc                                                                                                                           >  8:311438/1‑123 (MQ=255)
                                                  aTCCTGATCGCTTTGCTGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGt                                                                                                                                           <  8:129743/104‑1 (MQ=255)
                                                  aTCCTGATCGCTTTGCTGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGt                                                                                                                                           >  7:129743/1‑104 (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                                                                                     >  6:116706/1‑145 (MQ=255)
                                                               tGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACg                                                                                     <  5:116706/145‑1 (MQ=255)
                                                               tGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTc                                                                                  <  3:251708/148‑1 (MQ=255)
                                                                 ccGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                               >  6:116281/1‑149 (MQ=255)
                                                                 ccGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                               <  1:227475/149‑1 (MQ=255)
                                                                  cGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                               >  3:355726/1‑148 (MQ=255)
                                                                  cGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                               >  4:228249/1‑148 (MQ=255)
                                                                    aCGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGAc                                                                           <  2:144156/149‑1 (MQ=255)
                                                                        cTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCt                                                                        >  2:122152/1‑148 (MQ=255)
                                                                        cTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCt                                                                        >  4:276305/1‑148 (MQ=255)
                                                                         tCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGa                                                                      >  7:243324/1‑149 (MQ=255)
                                                                         tCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGa                                                                      >  4:119344/1‑149 (MQ=255)
                                                                           cGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACCCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGa                                                                      <  5:102097/147‑1 (MQ=255)
                                                                            gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGAccc                                                                         >  3:236495/1‑143 (MQ=255)
                                                                            gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGAccc                                                                         <  4:236495/143‑1 (MQ=255)
                                                                              gCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTcc                                                                 <  6:2823/149‑1 (MQ=255)
                                                                                 tGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAAGTCAGGGTCGTGGACCGTGACCTCctgg                                                              <  4:42706/149‑2 (MQ=255)
                                                                                   agcagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAAGTCAGGGTCGTGGACCGTGACCTCctggg                                                             <  5:359905/148‑3 (MQ=38)
                                                                                               cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCCTGACTTCCTGCTGCAGCAGATGCa                                                 <  1:4695/149‑1 (MQ=255)
                                                                                               cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAg                                                <  3:57819/149‑1 (MQ=255)
                                                                                                  tGGAATAATCCTCGGACCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTc                                                                                  <  7:130060/113‑1 (MQ=255)
                                                                                                    gAAAAATCCGGGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTgg                                             <  3:184993/147‑1 (MQ=255)
                                                                                                    gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                            <  5:114161/148‑1 (MQ=255)
                                                                                                    gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                            <  3:276305/148‑1 (MQ=255)
                                                                                                        aaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCTTCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                       <  1:122152/149‑1 (MQ=255)
                                                                                                        aaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                       <  4:244161/149‑1 (MQ=255)
                                                                                                          tCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                       >  3:57166/1‑147 (MQ=255)
                                                                                                          tCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                       >  6:46354/1‑147 (MQ=255)
                                                                                                           ccGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                    >  5:324424/1‑149 (MQ=255)
                                                                                                                gCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCg                                >  5:207271/1‑148 (MQ=255)
                                                                                                                 cccTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaa                              <  5:46354/149‑1 (MQ=255)
                                                                                                                  ccTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCg                                >  2:253259/1‑146 (MQ=255)
                                                                                                                    tCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                           >  7:80942/1‑149 (MQ=255)
                                                                                                                    tCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                           >  4:305267/1‑149 (MQ=255)
                                                                                                                      aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                           <  2:176114/147‑1 (MQ=255)
                                                                                                                      aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaagaa                           <  7:252751/147‑1 (MQ=255)
                                                                                                                       acgacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                        <  1:253259/149‑1 (MQ=255)
                                                                                                                       acgacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                        >  4:57419/1‑149 (MQ=255)
                                                                                                                          acgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAg                     <  3:228249/149‑1 (MQ=255)
                                                                                                                            gCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGa                   <  5:185419/149‑1 (MQ=255)
                                                                                                                              ttACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGAt                  >  5:232341/1‑148 (MQ=255)
                                                                                                                                gcgcgACGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATg                 >  6:129260/5‑147 (MQ=255)
                                                                                                                                 ccAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg              <  7:67491/149‑1 (MQ=255)
                                                                                                                                  cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGc             >  3:149916/1‑149 (MQ=255)
                                                                                                                                     aCGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa          <  6:137694/149‑1 (MQ=255)
                                                                                                                                     aCGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa          >  1:379586/1‑149 (MQ=255)
                                                                                                                                         tGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCtt      >  4:199825/1‑149 (MQ=255)
                                                                                                                                          gcgcACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTg     >  5:215368/1‑149 (MQ=255)
                                                                                                                                           cgcAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTg     <  6:207271/148‑1 (MQ=255)
                                                                                                                                            gcAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGAc   <  3:119344/149‑1 (MQ=255)
                                                                                                                                             cAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACg  <  3:305267/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACG  >  NC_002947/980612‑980901

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: