Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,229,387 T→A 7.1% intergenic (‑21/‑426) PP_1964 ← / → holB dGMP/dTMP/dhydroxyCMP kinase/DNA polymerase III subunit delta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,229,3870TA 7.1% 582.3 / 24.3 182intergenic (‑21/‑426)PP_1964/holBdGMP/dTMP/dhydroxyCMP kinase/DNA polymerase III subunit delta'
Reads supporting (aligned to +/- strand):  ref base T (77/92);  new base A (7/6);  total (84/98)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 184 total aligned reads displayed.

TGCTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCGTATC  >  NC_002947/2229238‑2229535
                                                                                                                                                     |                                                                                                                                                    
tGCTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAt                                                                                                                                                      >  1:4086610/1‑150 (MQ=255)
   ttGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTc                                                                                                                                                   <  1:2664626/150‑1 (MQ=255)
     ggCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCgg                                                                                                                                                 >  2:5249648/1‑150 (MQ=255)
      gCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                                <  1:474592/150‑1 (MQ=255)
      gCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGGCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                                <  1:5005767/150‑1 (MQ=255)
      gCAACCTGCTTTGAGCGGATCAGCTAGGGATTATGTCGCGACCTCAGGGTGCTCAATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                                <  1:517604/150‑1 (MQ=255)
       cAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCgggg                                                                                                                                               <  1:3186433/150‑1 (MQ=255)
         aCCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGt                                                                                                                                             <  1:1839643/150‑1 (MQ=255)
          ccTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGtt                                                                                                                                            <  2:4982820/150‑1 (MQ=255)
          ccTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                            <  1:1057077/150‑4 (MQ=255)
          ccTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                            <  2:2712890/150‑4 (MQ=255)
             gCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTaga                                                                                                                                         <  1:3805627/150‑1 (MQ=255)
              cTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTagag                                                                                                                                        <  1:3952364/150‑1 (MQ=255)
                 tGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGa                                                                                                                                     <  2:1770703/150‑1 (MQ=255)
                     cGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTa                                                                                                                                 <  2:1555820/150‑1 (MQ=255)
                        aTCAGCTAGGGCATATGTCGCGACCTCCGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGaa                                                                                                                              <  2:5397632/150‑1 (MQ=255)
                        aTCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGaa                                                                                                                              <  1:4796895/150‑1 (MQ=255)
                        aTCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGaa                                                                                                                              <  1:1808927/150‑1 (MQ=255)
                           aGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTc                                                                                                                           <  1:2118819/150‑1 (MQ=255)
                           aGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTc                                                                                                                           <  1:1459260/150‑1 (MQ=255)
                            gCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCa                                                                                                                          <  1:5421332/150‑1 (MQ=255)
                             cTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCAt                                                                                                                         <  1:1082726/150‑1 (MQ=255)
                             cTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcctgga                                                                                                                                         <  2:1835255/134‑7 (MQ=255)
                              tAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCAtt                                                                                                                        >  2:723303/1‑149 (MQ=255)
                                gggCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGACTTTCGGGGGTTAGAGGGATGTAGAACTCATGtt                                                                                                                      >  2:3859919/1‑150 (MQ=255)
                                gggCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGtt                                                                                                                      <  1:3645688/150‑1 (MQ=255)
                                        gTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCAt                                                                                                                         >  1:3933293/1‑139 (MQ=255)
                                        gTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCAt                                                                                                                         <  2:3933293/139‑1 (MQ=255)
                                          cgcgACCTCAGGGTGCTCCATGAGCATTGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                            <  1:3287093/150‑1 (MQ=255)
                                          cgcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                            <  2:1636815/150‑1 (MQ=255)
                                          cgcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                            >  2:55073/1‑150 (MQ=255)
                                          cgcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                            <  2:3154372/150‑1 (MQ=255)
                                           gcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAg                                                                                                           >  1:2457293/1‑150 (MQ=255)
                                             gACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTGTTAGTa                                                                                                         >  1:3718153/1‑150 (MQ=255)
                                              aCCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTaa                                                                                                        <  1:2081352/150‑1 (MQ=255)
                                               ccTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCAt                                                                                                                         <  1:3098101/132‑1 (MQ=255)
                                                cTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                            >  2:5155922/1‑109 (MQ=255)
                                                cTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                            <  1:5155922/112‑4 (MQ=255)
                                                 tCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGa                                                                                                     <  2:3981462/150‑1 (MQ=255)
                                                 tCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGa                                                                                                     >  1:702193/1‑150 (MQ=255)
                                                     ggTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGt                                                                                                 <  2:2457293/150‑1 (MQ=255)
                                                     ggTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGt                                                                                                 >  1:86953/1‑150 (MQ=255)
                                                      gTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGta                                                                                                <  1:3457339/150‑1 (MQ=255)
                                                          tCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTataat                                                                                            >  2:2050514/1‑150 (MQ=255)
                                                           ccATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACCCATGTTTACATTCTTAGTAAGGAACGTataata                                                                                           <  2:4515758/150‑1 (MQ=255)
                                                            cATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAg                                                                                          >  2:3713712/1‑150 (MQ=255)
                                                            cATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATGTTTGGCTCCAGGCCCCGAATTTCGGGGGTGAGAGGGATGGAGAACTCATGTTTACATTCTTAGTAAGGAACGTATACTAg                                                                                          >  1:128696/1‑150 (MQ=255)
                                                              tGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGc                                                                                        <  2:1177701/150‑1 (MQ=255)
                                                                   aTGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATg                                                                                   >  1:552690/1‑150 (MQ=255)
                                                                      gAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGcc                                                                                       <  2:843559/143‑1 (MQ=255)
                                                                      gAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGcc                                                                                       >  1:843559/1‑143 (MQ=255)
                                                                      gAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGttt                                                                                >  1:4404519/1‑150 (MQ=255)
                                                                      gAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGttt                                                                                <  1:3713712/150‑1 (MQ=255)
                                                                       aaGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTa                                                                               >  1:4764244/1‑150 (MQ=255)
                                                                         gCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACa                                                                             >  2:2253572/1‑150 (MQ=255)
                                                                         gCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACa                                                                             <  2:1444886/150‑1 (MQ=255)
                                                                               gtgtCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACt                                                                       <  2:3394838/150‑1 (MQ=255)
                                                                                tgtCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTc                                                                      <  2:3470611/150‑1 (MQ=255)
                                                                                tgtCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTc                                                                      <  1:5233258/150‑1 (MQ=255)
                                                                                 gtCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGACGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTAt                                                                                    <  1:582935/135‑1 (MQ=255)
                                                                                   cTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAg                                                                   >  2:136487/1‑150 (MQ=255)
                                                                                   cTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAg                                                                   >  2:136493/1‑150 (MQ=255)
                                                                                   cTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAg                                                                   <  2:4195595/150‑1 (MQ=255)
                                                                                    tttACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGa                                                                  >  1:2697139/1‑150 (MQ=255)
                                                                                    tttACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGa                                                                  >  1:2483906/1‑150 (MQ=255)
                                                                                     ttACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTATGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGaa                                                                 <  1:5095544/150‑1 (MQ=255)
                                                                                      tACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAAc                                                                >  2:2335274/1‑150 (MQ=255)
                                                                                        ccGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGa                                                              >  1:2359132/1‑150 (MQ=255)
                                                                                         cGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGaa                                                             <  1:4271720/150‑1 (MQ=255)
                                                                                           gATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATa                                                           <  1:2335274/150‑1 (MQ=255)
                                                                                             tCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACa                                                         <  1:2952833/150‑1 (MQ=255)
                                                                                              cGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACaa                                                        <  1:4690792/150‑1 (MQ=255)
                                                                                                aGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGa                                                      <  2:4703606/150‑1 (MQ=255)
                                                                                                 tcTACTGCGGCGAGGACGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTATAAATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGAc                                                     >  2:2950603/2‑150 (MQ=255)
                                                                                                 gCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGATCTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGAc                                                     <  2:3598151/150‑1 (MQ=255)
                                                                                                   tACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGc                                                   >  1:4860194/1‑150 (MQ=255)
                                                                                                        cggcgAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAGTTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCt                                              >  2:3321447/1‑150 (MQ=255)
                                                                                                            gAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGc                                          <  1:764281/150‑1 (MQ=255)
                                                                                                            gAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGc                                          <  1:3498684/150‑1 (MQ=255)
                                                                                                            gAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGc                                          <  2:1120150/150‑1 (MQ=255)
                                                                                                                  gatgatGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGc                                                   >  1:4261352/1‑135 (MQ=255)
                                                                                                                  gatgatGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCa                                    >  1:442708/1‑150 (MQ=255)
                                                                                                                   atgatgGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                            <  2:3848508/77‑1 (MQ=255)
                                                                                                                      atgGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGggc                                >  2:4064471/1‑150 (MQ=255)
                                                                                                                             gCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATACTAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTc                         >  2:2349933/1‑150 (MQ=255)
                                                                                                                                tttttGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTg                      <  1:4593409/150‑1 (MQ=255)
                                                                                                                                tttttGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTg                      >  2:4813079/1‑150 (MQ=255)
                                                                                                                                       cTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtc               >  1:3522006/1‑150 (MQ=255)
                                                                                                                                       cTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtc               <  1:5249648/150‑1 (MQ=255)
                                                                                                                                       cTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtc               >  2:1873492/1‑150 (MQ=255)
                                                                                                                                        tCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcg              >  2:439860/1‑150 (MQ=255)
                                                                                                                                         ctaacccccGAAATTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATCGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcgg             >  2:3214613/6‑150 (MQ=255)
                                                                                                                                          cAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcggt            <  1:3571070/150‑1 (MQ=255)
                                                                                                                                            acccccGAAATTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcggtcg          >  2:2123072/3‑150 (MQ=255)
                                                                                                                                             gCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGa         >  2:2194557/1‑150 (MQ=255)
                                                                                                                                              ccccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAg        <  1:1641431/150‑1 (MQ=255)
                                                                                                                                               cccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGc       >  2:4288353/1‑150 (MQ=255)
                                                                                                                                               cccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGc       >  2:816724/1‑150 (MQ=255)
                                                                                                                                                ccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCg      >  2:3744541/1‑150 (MQ=255)
                                                                                                                                                 cGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCGt     >  1:2435107/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCGTATC  >  NC_002947/2229238‑2229535

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