Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,723,602 G→T 6.3% intergenic (‑39/+25) PP_4178 ← / ← htpG dienelactone hydrolase family protein/chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,723,6020GT 6.3% 626.6 / 22.9 189intergenic (‑39/+25)PP_4178/htpGdienelactone hydrolase family protein/chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base G (91/86);  new base T (5/7);  total (96/93)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 195 total aligned reads displayed.

GCCATCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAA  >  NC_002947/4723454‑4723747
                                                                                                                                                     |                                                                                                                                                 
     cACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCgg                                                                                                                                             <  2:97407/150‑1 (MQ=255)
       ccAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGc                                                                                                                                           >  1:709926/1‑150 (MQ=255)
        cAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCTTCATAAAAAAGCCCGCCGAAGCGGGCt                                                                                                                                          >  2:2675790/1‑150 (MQ=255)
         aGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAGAAAGCCCGCCGAAGCGGGCtt                                                                                                                                         >  2:4736251/1‑150 (MQ=255)
         aGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCtt                                                                                                                                         <  1:5134845/150‑1 (MQ=255)
          gCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCttt                                                                                                                                        <  1:305891/150‑1 (MQ=255)
          gCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCttt                                                                                                                                        <  1:2701262/150‑1 (MQ=255)
           cccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTg                                                                                                                                       <  2:2195468/150‑1 (MQ=255)
             cTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCt                                                                                                                                     <  2:3354456/150‑1 (MQ=255)
              tGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtttttta                                                                                                                                    <  1:3593715/150‑5 (MQ=255)
              tGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCtg                                                                                                                                    >  2:5432742/1‑150 (MQ=255)
              tGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCtg                                                                                                                                    >  2:1464819/1‑150 (MQ=255)
                 gcggcCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGc                                                                                                                                 <  1:195301/150‑1 (MQ=255)
                  cggcCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCa                                                                                                                                <  2:53206/150‑1 (MQ=255)
                      ccAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCaa                                                                                                                            <  1:5252744/150‑1 (MQ=255)
                       cAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAAc                                                                                                                           >  1:1037175/1‑150 (MQ=255)
                             aCCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGc                                                                                                                     <  2:5290363/150‑1 (MQ=255)
                              cccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCa                                                                                                                    <  1:5187387/150‑1 (MQ=255)
                              cccGGCTCATAGACCAAACGGCTTTAGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATTAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGc                                                                                                                     <  2:2308959/150‑1 (MQ=255)
                               ccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAg                                                                                                                   >  1:1575322/1‑150 (MQ=255)
                               ccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAg                                                                                                                   >  2:2086163/1‑150 (MQ=255)
                                  gCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGTATCAACTCAAGCAGACa                                                                                                                <  2:3024538/150‑1 (MQ=255)
                                  gCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACa                                                                                                                <  1:2086163/150‑1 (MQ=255)
                                      aTAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                            <  1:1464701/150‑1 (MQ=255)
                                      aTAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                            <  1:298662/150‑1 (MQ=255)
                                          aCCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTTGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGAcc                                                                                                        >  2:1057418/1‑150 (MQ=255)
                                             aaaCGGCTTTCGAAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGt                                                                                                     >  1:4913092/1‑150 (MQ=255)
                                              caCGGCTTTCGTAGGTCGTGCCAGACAGGTGATACACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCATTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTa                                                                                                    <  2:3062562/149‑1 (MQ=255)
                                              aaCGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTa                                                                                                    >  2:5384103/1‑150 (MQ=255)
                                              aaCGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTa                                                                                                    >  1:2274602/1‑150 (MQ=255)
                                                   cTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTg                                                                                               <  2:2408699/150‑1 (MQ=255)
                                                    tttCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAg                                                                                                                      >  2:51702/1‑126 (MQ=255)
                                                    tttCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAg                                                                                                                      <  1:51702/126‑1 (MQ=255)
                                                    tttCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGt                                                                                              >  1:1839614/1‑150 (MQ=255)
                                                    tttCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGt                                                                                              >  2:1366586/1‑150 (MQ=255)
                                                    gttCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGt                                                                                              >  1:127234/2‑150 (MQ=255)
                                                      tCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCCTCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTc                                                                                            >  1:4397678/1‑150 (MQ=255)
                                                       cGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCa                                                                                           <  1:4052306/150‑1 (MQ=255)
                                                         tAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAgg                                                                                         >  1:2898107/1‑150 (MQ=255)
                                                          aGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGTACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                        <  2:388172/150‑1 (MQ=255)
                                                          aGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                        >  2:255132/1‑150 (MQ=255)
                                                          aGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                        >  2:4482021/1‑150 (MQ=255)
                                                           ggTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtt                                                                                                                                         <  2:1739947/100‑1 (MQ=255)
                                                           ggTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtt                                                                                                                                         >  1:1739947/1‑100 (MQ=255)
                                                           ggTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcg                                                                                       >  2:4200192/1‑150 (MQ=255)
                                                           ggTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcg                                                                                       >  2:75124/1‑150 (MQ=255)
                                                            gTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcga                                                                                      >  1:3670756/1‑150 (MQ=255)
                                                            gTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcga                                                                                      <  1:2025373/150‑1 (MQ=255)
                                                              cTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcgacg                                                                                    >  1:2637584/1‑150 (MQ=255)
                                                               ttGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcgacga                                                                                   >  2:4174371/1‑150 (MQ=255)
                                                                     gACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTaa                                                                             <  1:1390335/150‑1 (MQ=255)
                                                                     gACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTaa                                                                             <  2:2018434/150‑1 (MQ=255)
                                                                            gATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGccgccg                                                                      >  1:4744040/1‑150 (MQ=255)
                                                                             aTAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCgg                                                                     >  1:1629019/1‑150 (MQ=255)
                                                                                  ccAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                            <  1:4338676/106‑1 (MQ=255)
                                                                                  ccAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCt                                                                >  1:2564969/1‑150 (MQ=255)
                                                                                  ccAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCt                                                                >  2:2036390/1‑150 (MQ=255)
                                                                                   cAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCtt                                                               >  2:2946465/1‑150 (MQ=255)
                                                                                   cAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCtt                                                               >  2:1564615/1‑150 (MQ=255)
                                                                                     gCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCa                                                             <  2:5232229/150‑1 (MQ=255)
                                                                                         ttCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCt                                                         <  2:2679475/150‑1 (MQ=255)
                                                                                            gacgaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                            <  1:966209/96‑1 (MQ=255)
                                                                                            gacgaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGt                                                                                                     >  2:664093/1‑103 (MQ=255)
                                                                                            gacgaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGt                                                                                                     <  1:664093/103‑1 (MQ=255)
                                                                                               gaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCaa                                                                                                                       <  1:4195847/82‑1 (MQ=255)
                                                                                               gaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCAcc                                                   <  2:127202/150‑1 (MQ=255)
                                                                                                     tttGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCaa                                                                                                                       >  1:4572941/1‑76 (MQ=255)
                                                                                                     tttGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCaa                                                                                                                       <  2:4572941/76‑1 (MQ=255)
                                                                                                     tttGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGgccgc                                             >  1:2096736/1‑150 (MQ=255)
                                                                                                     tttGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGgccgc                                             >  1:947560/1‑150 (MQ=255)
                                                                                                          tcatcaAGCACTCCTTTGGCACTCGTCATAAAATAGCCCGCCGAAGCGGGCTTTGCTGTGCATCATCTCAAGCAGTCAGCTCGACCAGTAGCTTGTTCAGTCGACGTACGTAAGCCGCCGGGTCCTTCAGGCTGTCAGCGGCCGCCAggg                                        >  2:4986784/1‑150 (MQ=255)
                                                                                                          tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                        >  2:4424330/1‑150 (MQ=255)
                                                                                                          tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                        >  2:4655872/1‑150 (MQ=255)
                                                                                                          tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                        >  1:2635570/1‑150 (MQ=255)
                                                                                                          tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggc                                        <  1:2836450/150‑2 (MQ=255)
                                                                                                          tcatcaAGCACTCCTTTGGCACTCGACATAAAAAAGCCCGCCGAAGCGGGCTTTGATGTGCAACAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCATTCAGGCTGTCACCGGCCGCCAggg                                        >  2:702378/1‑150 (MQ=255)
                                                                                                                aGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCt                                  <  1:744427/150‑1 (MQ=255)
                                                                                                                aGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCt                                  >  2:4301337/1‑150 (MQ=255)
                                                                                                                aGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCt                                  >  2:1508896/1‑150 (MQ=255)
                                                                                                                 gCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTg                                 >  2:4818653/1‑150 (MQ=255)
                                                                                                                    cTCCGTGGGCACTCGTCATAACAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCTGACAGCTCGCCCAGTAGCTTGTTCAGGCGACGACCGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATc                              <  2:3729014/150‑1 (MQ=255)
                                                                                                                       cTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGaa                           >  2:168611/1‑150 (MQ=255)
                                                                                                                         ttGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACa                         <  2:4234644/150‑1 (MQ=255)
                                                                                                                             cACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAt                     <  1:1508896/150‑1 (MQ=255)
                                                                                                                             cACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAt                     >  2:4773091/1‑150 (MQ=255)
                                                                                                                             cACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAt                     <  1:3884304/150‑1 (MQ=255)
                                                                                                                                tCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAtgtg                  <  1:5384103/150‑1 (MQ=255)
                                                                                                                                tCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAtgtg                  <  1:1751256/150‑1 (MQ=255)
                                                                                                                                tCGTCACAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTTTTCAGGCGCCGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAtgtg                  <  1:769426/150‑1 (MQ=255)
                                                                                                                                 cGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGc                 <  2:4909242/150‑1 (MQ=255)
                                                                                                                                 cGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGc                 <  2:2813706/150‑1 (MQ=255)
                                                                                                                                        agcaaaGCCCGCTTCGGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGc           >  2:2102747/4‑150 (MQ=255)
                                                                                                                                        aaaaaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGc           >  2:598586/1‑150 (MQ=255)
                                                                                                                                           aaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCg        <  1:3781435/150‑1 (MQ=255)
                                                                                                                                              gCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCg     >  1:290607/1‑150 (MQ=255)
                                                                                                                                              gCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCg     >  1:1203845/1‑150 (MQ=255)
                                                                                                                                               cccGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCgg                                                                     >  1:1440499/1‑85 (MQ=255)
                                                                                                                                               cccGCCGAAGCGGGCGGTGCTGTGCATCAACTCAAGCAGACCGCTCGACCCGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACGGGATGTGCGACAGCTCGGCGa    >  2:4469531/1‑150 (MQ=255)
                                                                                                                                                 cgccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGaaa  <  1:2675839/150‑1 (MQ=255)
                                                                                                                                                 cgccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGaaa  <  1:2675790/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCCATCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAA  >  NC_002947/4723454‑4723747

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