Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,955,480 A→C L175R (CTG→CGG)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,955,4800AC97.8% 326.9 / ‑1.2 90L175R (CTG→CGG) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (42/46);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAA  >  NC_000913/2955334‑2955625
                                                                                                                                                  |                                                                                                                                                 
gcgcGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCcgcg                                                                                                                                                <  1:1078467/150‑1 (MQ=255)
gcgcGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCcgcg                                                                                                                                                <  1:1015283/150‑1 (MQ=255)
   cgTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTa                                                                                                                                            <  2:1007304/151‑1 (MQ=255)
      gACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCa                                                                                                                                         <  1:535528/151‑1 (MQ=255)
         gggAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTg                                                                                                                                      <  2:1109681/151‑1 (MQ=255)
         gggAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTg                                                                                                                                      >  2:306620/1‑151 (MQ=255)
            aaGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGc                                                                                                                                    <  2:335810/150‑1 (MQ=255)
            aaGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCg                                                                                                                                   >  1:448470/1‑151 (MQ=255)
              gCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGAc                                                                                                                                 <  2:915883/151‑1 (MQ=255)
              gCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGAc                                                                                                                                 <  2:792812/151‑1 (MQ=255)
              gCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGAc                                                                                                                                 >  1:209069/1‑151 (MQ=255)
               ccAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACg                                                                                                                                >  1:249326/1‑151 (MQ=255)
                  gCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGcc                                                                                                                              >  1:1063896/1‑150 (MQ=255)
                  gCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCa                                                                                                                             <  2:171362/151‑1 (MQ=255)
                   cGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAg                                                                                                                            <  2:262594/151‑1 (MQ=255)
                    gTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGcc                                                                                                                              >  1:1146818/1‑148 (MQ=255)
                      ttCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGaa                                                                                                                          >  1:319985/1‑150 (MQ=255)
                         atcatcGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCg                                                                                                                      <  2:959798/151‑1 (MQ=255)
                          tcatcGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGCTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCgg                                                                                                                     >  2:509805/1‑151 (MQ=255)
                            atcGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGgcg                                                                                                                   <  1:819753/151‑1 (MQ=255)
                                ggcggcGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAgg                                                                                                               <  1:373501/151‑1 (MQ=255)
                                  cggcggTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGcc                                                                                                             >  1:893990/1‑151 (MQ=255)
                                        tGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgt                                                                                                         >  2:241676/1‑149 (MQ=255)
                                         gCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtag                                                                                                       <  1:765805/150‑1 (MQ=255)
                                               aTAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCa                                                                                                <  1:264278/151‑1 (MQ=255)
                                                 aaCCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtagcgtagc                                                                                                 >  1:1292797/1‑148 (MQ=255)
                                                    cGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGAc                                                                                           >  1:622703/1‑151 (MQ=255)
                                                    cGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGAc                                                                                           <  2:771233/151‑1 (MQ=255)
                                                      gcgcTTCGCCTTACAGATAACGATTAATATCGCGCAGCAACGCCTGCGACCCTTTCCAGGTTTCAAAGACAACCTGAGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACt                                                                                          <  1:735312/150‑1 (MQ=255)
                                                       cgcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCa                                                                                        <  1:1213146/151‑1 (MQ=255)
                                                        gcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCAt                                                                                       >  2:972928/1‑151 (MQ=255)
                                                         cTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATc                                                                                      >  2:869206/1‑151 (MQ=255)
                                                          ttCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCt                                                                                     >  1:821535/1‑151 (MQ=255)
                                                          ttCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCt                                                                                     >  1:328515/1‑151 (MQ=255)
                                                          ttCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCt                                                                                     >  2:605872/1‑151 (MQ=255)
                                                           tCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCtt                                                                                    >  2:1028831/1‑151 (MQ=255)
                                                            cGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTc                                                                                   >  1:24668/1‑151 (MQ=255)
                                                            cGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTc                                                                                   >  2:1292021/1‑151 (MQ=255)
                                                                  gCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAg                                                                             <  1:746605/151‑1 (MQ=255)
                                                                    aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgct                                                                           >  2:628883/1‑151 (MQ=255)
                                                                    aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTGTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgc                                                                            <  1:1172018/150‑1 (MQ=255)
                                                                      aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgctgc                                                                         <  2:644393/151‑1 (MQ=255)
                                                                       tAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgctgct                                                                        >  2:824333/1‑151 (MQ=255)
                                                                        aaCGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgctgctg                                                                       >  2:208053/1‑151 (MQ=255)
                                                                        aaCGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGAAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgctgct                                                                        >  2:669115/1‑150 (MQ=255)
                                                                          cGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAgctgctgct                                                                     <  2:755749/151‑1 (MQ=255)
                                                                           gATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACCCATCTTCAATCAGCTGCTGCTc                                                                    >  2:532600/1‑151 (MQ=255)
                                                                             ttAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCa                                                                   >  1:195518/1‑150 (MQ=255)
                                                                             ttAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCa                                                                   >  1:38194/1‑150 (MQ=255)
                                                                              tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCaa                                                                  >  2:732847/1‑150 (MQ=255)
                                                                                atatCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAAc                                                                <  1:509791/150‑1 (MQ=255)
                                                                                atatCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACa                                                               <  1:638479/151‑1 (MQ=255)
                                                                                  atCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAg                                                              <  2:558429/150‑1 (MQ=255)
                                                                                    cgcgCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCAt                                                           >  2:1232222/1‑151 (MQ=255)
                                                                                      cgcAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGc                                                         >  2:160485/1‑151 (MQ=255)
                                                                                      cgcAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGc                                                         >  1:469673/1‑151 (MQ=255)
                                                                                      cgcAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGc                                                         <  2:657417/151‑1 (MQ=255)
                                                                                      cgcAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGc                                                         >  1:55991/1‑151 (MQ=255)
                                                                                       gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGcc                                                        >  1:557458/1‑151 (MQ=255)
                                                                                           caACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGAt                                                    <  2:617464/151‑1 (MQ=255)
                                                                                            aaCGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGAtt                                                   <  2:68615/151‑1 (MQ=255)
                                                                                            aaCGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGAtt                                                   <  1:899296/151‑1 (MQ=255)
                                                                                               gCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCaa                                                <  2:749611/151‑1 (MQ=255)
                                                                                                  tGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAgg                                             >  1:436272/1‑151 (MQ=255)
                                                                                                   gCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAgg                                             <  1:298600/150‑1 (MQ=255)
                                                                                                   gCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAg                                              <  2:1107783/149‑1 (MQ=255)
                                                                                                    cGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCa                                           <  2:489589/151‑1 (MQ=255)
                                                                                                     ggCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCAt                                          >  2:55441/1‑151 (MQ=255)
                                                                                                        ccTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCa                                       <  1:646120/151‑1 (MQ=255)
                                                                                                          tttCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAgg                                     >  2:1281390/1‑151 (MQ=255)
                                                                                                           ttacAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGt                                    >  1:249814/4‑151 (MQ=255)
                                                                                                           ttCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAg                                      <  2:683864/149‑1 (MQ=255)
                                                                                                           ttCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAg                                      <  2:215881/149‑1 (MQ=255)
                                                                                                             ccAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTg                                  >  1:1183571/1‑151 (MQ=255)
                                                                                                              cAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGtt                                   <  1:1225/149‑1 (MQ=255)
                                                                                                                      aaaGACGACCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATgcg                          <  2:353985/150‑1 (MQ=255)
                                                                                                                         gacgacCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTgg                      >  1:4974/1‑151 (MQ=255)
                                                                                                                          acgacCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGc                     <  2:448481/151‑1 (MQ=255)
                                                                                                                           cgacCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCa                    <  1:1198801/151‑1 (MQ=255)
                                                                                                                            gacCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCaa                   <  2:783217/151‑1 (MQ=255)
                                                                                                                            gacCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCaa                   <  1:1281311/151‑1 (MQ=255)
                                                                                                                            gacCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCaa                   <  2:698149/151‑1 (MQ=255)
                                                                                                                             acCTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCaaa                  <  2:783582/151‑1 (MQ=255)
                                                                                                                               cTGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAACATGCGCTGGCaaa                  >  1:1288206/1‑149 (MQ=255)
                                                                                                                                tGGGCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGc               >  1:790001/1‑151 (MQ=255)
                                                                                                                                 gggCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGcc              >  2:956061/1‑151 (MQ=255)
                                                                                                                                 gggCTATTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGcc              <  1:756960/151‑1 (MQ=255)
                                                                                                                                      aTTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTg           <  1:234537/149‑1 (MQ=255)
                                                                                                                                      aTTTCACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGa          >  2:800806/1‑150 (MQ=255)
                                                                                                                                          cACGCGGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGt     <  2:404659/151‑1 (MQ=255)
                                                                                                                                            cgcgGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTaa   <  2:980060/151‑1 (MQ=255)
                                                                                                                                             gcgGCCGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTaaa  <  1:1130002/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAA  >  NC_000913/2955334‑2955625

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