Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,928 C→G 5.3% G501R (GGC→CGC)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,9280CG 5.3% 323.2 / 7.6 95G501R (GGC→CGC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base C (43/47);  new base G (2/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

TAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960780‑2961076
                                                                                                                                                    |                                                                                                                                                    
tAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATAGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGcc                                                                                                                                                      >  1:777980/1‑149 (MQ=255)
    gTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAt                                                                                                                                                  >  1:893173/1‑149 (MQ=255)
    gTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGGTGTCGTCATCTATGCCCCAACGAATGCCGGAt                                                                                                                                                  >  1:160869/1‑149 (MQ=255)
      aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTc                                                                                                                                                <  2:50984/149‑1 (MQ=255)
        aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGt                                                                                                                                              >  2:785632/1‑149 (MQ=255)
         acagATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                             <  2:1267063/146‑1 (MQ=255)
         aCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                             >  2:1320860/1‑149 (MQ=255)
           cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGa                                                                                                                                           >  2:1047067/1‑149 (MQ=255)
           cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGGGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGa                                                                                                                                           <  1:1100434/149‑1 (MQ=255)
            gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACACCATACGCGTCAGGCCAACTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggg                                                                                                                                          <  1:325081/149‑3 (MQ=255)
            gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggg                                                                                                                                          <  1:1178958/149‑3 (MQ=255)
                gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggggcat                                                                                                                                      <  2:1415996/149‑7 (MQ=255)
                 ccACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                     >  1:369882/1‑149 (MQ=255)
                 ccACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCCTCAGGCCAAATCGCCAGGTGTGTTGTACGGTGGCGGGGAGATCCAGCTCGCGAACGTTGTCGGCATCTTTGCCCCAACGCTTGCCGGATTCGTTGACCCACt                                                                                                                                     >  2:1522220/1‑149 (MQ=255)
                   aCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                                   >  2:1019480/1‑149 (MQ=255)
                    cTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGc                                                                                                                                  >  2:580713/1‑149 (MQ=255)
                       gCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTa                                                                                                                               >  1:87943/1‑149 (MQ=255)
                         ccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                             <  1:623714/149‑1 (MQ=255)
                          cTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACCTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCCCTGGCGttatt                                                                                                                            >  1:40415/1‑145 (MQ=255)
                                  tctcCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGCCCCACTGGCGTAAATAAcgccgc                                                                                                                    >  2:66641/1‑149 (MQ=255)
                                    tcCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTACGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGccc                                                                                                                  >  2:1025468/1‑149 (MQ=255)
                                      cATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCctt                                                                                                                >  1:54936/1‑149 (MQ=255)
                                           cgTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCg                                                                                                           <  1:1047067/149‑1 (MQ=255)
                                            gTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGGGTTGTCCGGGGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTTAATAACGCAGCCCTTCTTCgg                                                                                                          <  2:1016559/149‑1 (MQ=255)
                                             tAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGt                                                                                                         <  2:1098062/149‑1 (MQ=255)
                                              aGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGGGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTg                                                                                                        <  2:777980/149‑1 (MQ=255)
                                               gCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                       <  1:1267716/149‑1 (MQ=255)
                                               gCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                       >  1:468123/1‑149 (MQ=255)
                                                cccATCAACATACGCGGCAGGCCTAATCGGCATGTGTGTTGGCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGGCCCAACGAATGCCGGATTCGTTGACCCACTTGCGTAAATAACGCAGCCGTTCTTAGGTGAt                                                                                                      <  2:819886/149‑1 (MQ=255)
                                                cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGAt                                                                                                      <  2:893173/149‑1 (MQ=255)
                                                  caacaaCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGt                                                                                                    <  1:910296/149‑1 (MQ=255)
                                                  caacaaCATACGCGTCAGGCCAAATCGCCAGGTGTGTAGAACTCTGTCGTGGAGTTACATCTAGAGAACATTGTCGTCATCTATGCCCCAACGAATGACGGATTCTTTGACCCACTGGCGTAAATAACGAAGCCCGTCTTCGGTGCTGt                                                                                                    >  2:1437227/1‑149 (MQ=255)
                                                   aacaacATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                   >  1:933175/1‑149 (MQ=255)
                                                   aacaacATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                   <  1:614428/149‑1 (MQ=255)
                                                      atcaTACGCGTCAGGCCAATTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGGTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCTGTGATGTCaaa                                                                                                <  2:498603/147‑1 (MQ=255)
                                                      aacaTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTg                                                                                                                                    <  2:183093/113‑1 (MQ=255)
                                                      aacaTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTg                                                                                                                                    >  1:183093/1‑113 (MQ=255)
                                                      aacaTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCaaa                                                                                                >  1:1387744/1‑149 (MQ=255)
                                                      aacaTACGCGTCAGGCCAAATCGCAAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGGAGCCGGTCTTCGGTGATGTCaaa                                                                                                >  1:1048830/1‑149 (MQ=255)
                                                          tACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCTCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgc                                                                                            <  2:377400/149‑1 (MQ=255)
                                                           aCGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgcg                                                                                           >  1:1430619/1‑149 (MQ=255)
                                                               gTCTGGCCATATCGCCAGGTGTGTTGTCCGGTGGCGTGGAGTTCCAGCTCGAGAACGTTGTCGTCATCTATGACCAAACGAATGCCTGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccgc                                                                                       <  1:988115/149‑1 (MQ=255)
                                                                 cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCa                                                                                     <  1:580713/149‑1 (MQ=255)
                                                                    gCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTTTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCAACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCgccgcccggc                                                                                  >  2:781781/1‑145 (MQ=255)
                                                                    gCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACg                                                                                                                        <  2:1061880/111‑1 (MQ=255)
                                                                    gCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACg                                                                                                                        >  1:1061880/1‑111 (MQ=255)
                                                                    gCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCAAGATATAGAACGTTGTCGTCATCAATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCATCTTCGGTGATGTAAAACCGCGCCGCCAGcc                                                                                  >  2:459097/1‑148 (MQ=255)
                                                                     ccAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAc                                                                                 <  1:1478256/149‑1 (MQ=255)
                                                                       aaaTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCg                                                                               >  2:1362317/1‑149 (MQ=255)
                                                                          tCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGca                                                                            >  1:15017/1‑149 (MQ=255)
                                                                              cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGTCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATc                                                                        <  2:767169/149‑1 (MQ=255)
                                                                              cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATc                                                                        >  2:317140/1‑149 (MQ=255)
                                                                               aGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  1:1019480/149‑1 (MQ=255)
                                                                                     gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAAc                                                                 <  2:1020201/149‑1 (MQ=255)
                                                                                      ttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                                >  1:600956/1‑149 (MQ=255)
                                                                                      ttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                                >  1:60928/1‑149 (MQ=255)
                                                                                            ggTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGca                                                          <  1:781781/149‑1 (MQ=255)
                                                                                             gTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCTAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGcac                                                         <  1:785632/149‑1 (MQ=255)
                                                                                              tGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGcaca                                                        <  1:813060/149‑1 (MQ=255)
                                                                                              tGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGcaca                                                        >  2:1251049/1‑149 (MQ=255)
                                                                                                 cGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATcc                                                     >  1:1119979/1‑149 (MQ=255)
                                                                                                   gggAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCTGCACATCCAGCAACGCCAGCACATCctc                                                   <  1:1025468/149‑1 (MQ=255)
                                                                                                       gTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                               <  2:369882/149‑1 (MQ=255)
                                                                                                         tCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGacac                                             <  2:160869/149‑1 (MQ=255)
                                                                                                           cAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGCCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGCTGTCAAACCGCGCCGCCAGCCCCGGCACCTCCCGCAACGCCAGCACATCCTCTGACACaa                                           >  2:102362/1‑149 (MQ=255)
                                                                                                             gCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACaaaa                                         <  1:499161/149‑1 (MQ=255)
                                                                                                              cTCGCGAACGTTGTCTTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAAc                                        <  1:1431383/149‑1 (MQ=255)
                                                                                                              cTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAAc                                        <  1:1565599/149‑1 (MQ=255)
                                                                                                                  cgATCGTTTTCGTCATCTATGCCCCAACGAATGCCGTATTCGTTGACCCACTGGCGTAAATATCGCTGACCTTCTTGGGTGATGTCATACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACt                                    <  2:78848/149‑1 (MQ=255)
                                                                                                                   gAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCCGCACATCCTCTGACACAAAACGCCTg                                   >  1:546456/1‑149 (MQ=255)
                                                                                                                    aaCGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGt                                  <  2:87943/149‑1 (MQ=255)
                                                                                                                     aCGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTTACCCACTTGCGTAAATAACGCAGCCCTTCTTCTGTTATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTc                                 <  2:1134186/149‑1 (MQ=255)
                                                                                                                       gTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGGGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAg                               >  2:1383400/1‑149 (MQ=255)
                                                                                                                          gtcgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCa                            >  1:500758/1‑149 (MQ=255)
                                                                                                                          gtcgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCa                            >  1:557670/1‑149 (MQ=255)
                                                                                                                            cgtaATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCTTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGt                          <  1:1211409/149‑1 (MQ=255)
                                                                                                                               cATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt                       <  1:1010652/149‑1 (MQ=255)
                                                                                                                               aaTCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt                       <  1:553286/148‑1 (MQ=255)
                                                                                                                               aaTCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt                       <  1:1170261/148‑1 (MQ=255)
                                                                                                                                aTCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACt                                    <  2:737923/135‑1 (MQ=255)
                                                                                                                                aTCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACt                                    >  1:737923/1‑135 (MQ=255)
                                                                                                                                   tATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACa                   <  1:459097/149‑1 (MQ=255)
                                                                                                                                    aTGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAg                  >  2:1325052/1‑149 (MQ=255)
                                                                                                                                    aTGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAg                  >  2:1071452/1‑149 (MQ=255)
                                                                                                                                     tGCCCCAACGAATGCCGGTTTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGACATCTCATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg                 <  2:1048830/149‑1 (MQ=255)
                                                                                                                                     tGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg                 <  2:1387744/149‑1 (MQ=255)
                                                                                                                                      gCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGc                <  1:1251049/149‑1 (MQ=255)
                                                                                                                                       ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt               <  2:468123/149‑1 (MQ=255)
                                                                                                                                       ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt               >  2:1237619/1‑149 (MQ=255)
                                                                                                                                         tcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCCGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGATGa             >  1:576965/2‑149 (MQ=255)
                                                                                                                                         tcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa             >  2:912791/2‑149 (MQ=255)
                                                                                                                                          cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt            <  2:600956/149‑1 (MQ=255)
                                                                                                                                          cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt            >  1:508633/1‑149 (MQ=255)
                                                                                                                                            aCGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATaa          >  1:1045463/1‑149 (MQ=255)
                                                                                                                                              gAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAAc        <  2:933175/149‑1 (MQ=255)
                                                                                                                                                 tGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCt     >  2:7858/1‑149 (MQ=255)
                                                                                                                                                    cGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  >  1:1215804/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960780‑2961076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑