Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,723,510 T→C 6.5% Q18Q (CAA→CAG)  htpG ← chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,723,5100TC 6.5% 176.2 / 2.8 61Q18Q (CAA→CAG) htpGchaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

TTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCT  >  NC_002947/4723364‑4723649
                                                                                                                                                  |                                                                                                                                           
ttCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGTGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTg                                                                                                                                           >  3:84297/1‑149 (MQ=255)
 tCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTAGAGCGAGTGAATCATGAGGTGCAGCAGTTGc                                                                                                                                          >  8:213771/1‑149 (MQ=255)
    ggTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTc                                                                                                                                       <  1:56795/149‑1 (MQ=255)
     gTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCa                                                                                                                                      >  8:166038/1‑149 (MQ=255)
          cATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTca                                                                                                                                 >  1:219086/1‑148 (MQ=255)
            tCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTcatg                                                                                                                               <  2:219086/149‑4 (MQ=255)
              gCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTc                                                                                                                              <  2:178597/148‑1 (MQ=255)
              gCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCt                                                                                                                             <  6:176484/149‑1 (MQ=255)
                    gAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAg                                                                                                                        <  2:18959/148‑1 (MQ=255)
                    gAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGc                                                                                                                       >  7:28349/1‑149 (MQ=255)
                        aGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAg                                                                                                                   <  8:229954/149‑1 (MQ=255)
                        aGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAg                                                                                                                   >  4:166848/1‑149 (MQ=255)
                             tGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGttt                                                                                                              <  7:166038/149‑1 (MQ=255)
                                    gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCtttt                                                                                                        <  5:298489/148‑1 (MQ=255)
                                    gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtg                                                                                                       <  5:243016/149‑1 (MQ=255)
                                    gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtg                                                                                                       >  6:243016/1‑149 (MQ=255)
                                     ccAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                      <  6:255368/149‑1 (MQ=255)
                                     ccAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                      <  3:166848/149‑1 (MQ=255)
                                       aGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTgtgt                                                                                                    >  5:49938/1‑149 (MQ=255)
                                        gggCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGtt                                                                                                   <  8:165564/149‑1 (MQ=255)
                                        gggCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGtt                                                                                                   >  4:101983/1‑149 (MQ=255)
                                         ggCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGttt                                                                                                  <  2:43514/149‑1 (MQ=255)
                                           cTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTc                                                                                                 >  1:116072/1‑148 (MQ=255)
                                           cTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTc                                                                                                 >  5:284090/1‑148 (MQ=255)
                                              cGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCtttt                                                                                                        <  2:167263/138‑1 (MQ=255)
                                              cGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCtttt                                                                                                        >  1:167263/1‑138 (MQ=255)
                                                  aCGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCa                                                                                         >  2:244332/1‑149 (MQ=255)
                                                    gCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATg                                                                                       <  3:101983/149‑1 (MQ=255)
                                                     cAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCt                                                                                                           >  6:58699/1‑128 (MQ=255)
                                                              ggcggcGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCg                                                                                                                                 <  2:110849/97‑1 (MQ=255)
                                                              ggcggcGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCg                                                                                                                                 >  1:110849/1‑97 (MQ=255)
                                                               gcggcgTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAAc                                                                            >  1:116169/1‑149 (MQ=255)
                                                                  gcgTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTcc                                                                         >  3:263356/1‑149 (MQ=255)
                                                                    gTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCGTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGa                                                                       <  8:28349/149‑1 (MQ=255)
                                                                       ggAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                      <  8:281075/115‑1 (MQ=255)
                                                                       ggAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                      >  7:281075/1‑115 (MQ=255)
                                                                           gCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATg                                                                >  5:69435/1‑149 (MQ=255)
                                                                              ttGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGCTCTTCGAAACTCCGATCTGATGGCg                                                             >  1:66555/1‑149 (MQ=255)
                                                                                 gAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTg                                                          <  4:84297/149‑1 (MQ=255)
                                                                                   aaTCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGgc                                                        <  6:69435/149‑1 (MQ=255)
                                                                                     tCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGgcgc                                                      <  2:116072/149‑1 (MQ=255)
                                                                                       aGCTCCCGCAGGAAGATTTCCTTGTTCGAGTATAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCt                                                    <  7:213771/149‑1 (MQ=255)
                                                                                              gCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTggcggc                                              >  1:159319/1‑148 (MQ=255)
                                                                                               cAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCg                                            <  6:284090/149‑1 (MQ=255)
                                                                                                aGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGCTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCgg                                           <  2:66555/149‑1 (MQ=255)
                                                                                                  gAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGtt                                         <  1:244332/149‑1 (MQ=255)
                                                                                                    aGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCg                                       <  2:116169/149‑1 (MQ=255)
                                                                                                       tttCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                     <  8:157720/148‑1 (MQ=255)
                                                                                                           cTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCg                                >  5:256788/1‑149 (MQ=255)
                                                                                                                tCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTaa                            <  6:256788/148‑1 (MQ=255)
                                                                                                                  gAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTggc                                                 <  8:113921/125‑1 (MQ=255)
                                                                                                                  gAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTggc                                                 >  7:113921/1‑125 (MQ=255)
                                                                                                                     tACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAAcggcg                       >  1:288050/1‑148 (MQ=255)
                                                                                                                        aGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATg                   >  3:15328/1‑149 (MQ=255)
                                                                                                                                  gaccgAGGTGAAGCAACTGCTGCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCg                                            <  8:151596/110‑1 (MQ=39)
                                                                                                                                  gaccgAGGTGAAGCAACTGCTGCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCg                                            >  7:151596/5‑114 (MQ=39)
                                                                                                                                   cATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATg        >  4:148849/1‑149 (MQ=255)
                                                                                                                                    aTGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATgg       <  4:15328/149‑1 (MQ=255)
                                                                                                                                     tGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTg                                                          >  8:108441/1‑97 (MQ=255)
                                                                                                                                     tGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTg                                                          <  7:108441/97‑1 (MQ=255)
                                                                                                                                         gTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCt  >  4:244383/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCT  >  NC_002947/4723364‑4723649

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