Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,039 C→T 5.2% R403C (CGC→TGC)  imuB → DNA repair DNA linked protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,0390CT 5.2% 291.1 / 4.2 96R403C (CGC→TGC) imuBDNA repair DNA linked protein
Reads supporting (aligned to +/- strand):  ref base C (38/53);  new base T (3/2);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGGCC  >  NC_002947/3527893‑3528183
                                                                                                                                                  |                                                                                                                                                
gggCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTAAAGGGGCATGGCGCAGCCGCCGATCACCGCGCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCTCAAGGCAACATGCCGGTTGCACCGGGCAgccacc                                                                                                                                               >  2:985437/1‑150 (MQ=255)
    gcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                           <  1:2714485/150‑1 (MQ=255)
    gcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                           >  2:1345270/1‑150 (MQ=255)
    gcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                           <  1:1789555/150‑1 (MQ=255)
    gcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                           <  1:2811067/150‑1 (MQ=255)
      aaTTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGt                                                                                                                                         >  1:513600/1‑150 (MQ=255)
        ttGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTgg                                                                                                                                       <  2:1416191/150‑1 (MQ=255)
          gcgcGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCt                                                                                                                                     <  1:204031/150‑1 (MQ=255)
          gcgcGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCt                                                                                                                                     <  1:1201146/150‑1 (MQ=255)
             cgAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgtt                                                                                                                                  <  1:448108/150‑1 (MQ=255)
               aGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgttgt                                                                                                                                <  1:1345270/150‑1 (MQ=255)
               aGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgttgt                                                                                                                                <  2:740286/150‑1 (MQ=255)
                gCGGTTGCCCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGAGCCGGGCATTTGGCGGAGCAGGGTGCCCANGGCAACATGNCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTc                                                                                                                               <  1:985437/150‑1 (MQ=255)
                 cGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGAACACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCg                                                                                                                              <  2:371308/150‑1 (MQ=255)
                     tGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCAc                                                                                                                          <  2:712755/150‑1 (MQ=255)
                     tGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCAc                                                                                                                          <  2:651718/150‑1 (MQ=255)
                        gcgcCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCACCCCGGGCGGCtgcccg                                                                                                                                >  2:407717/1‑137 (MQ=255)
                        gcgcCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCACCCCGGGCGGCtgcccg                                                                                                                                <  1:407717/141‑5 (MQ=255)
                         cgcccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGAt                                                                                                                      <  1:906879/150‑1 (MQ=255)
                          gcccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATc                                                                                                                     <  2:513600/150‑1 (MQ=255)
                          gcccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATc                                                                                                                     <  2:2602067/150‑1 (MQ=255)
                            ccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATc                                                                                                                     >  1:2388911/1‑148 (MQ=255)
                            ccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATc                                                                                                                     <  2:2388911/148‑1 (MQ=255)
                              gccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCGGCTGCCCggt                                                                                                                                           <  1:1328640/124‑2 (MQ=255)
                              gccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCGGCTGCCCggt                                                                                                                                           >  2:1328640/1‑123 (MQ=255)
                               ccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCt                                                                                                                <  2:2838120/150‑1 (MQ=255)
                                 tGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTgg                                                                                                              <  2:2807026/150‑1 (MQ=255)
                                   ggggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGTCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGc                                                                                                                                      >  2:1324583/1‑124 (MQ=255)
                                   ggggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGTCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGc                                                                                                                                      <  1:1324583/124‑1 (MQ=255)
                                     ggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                           >  2:861413/1‑117 (MQ=255)
                                     ggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCggg                                                                                                                                           <  1:861413/117‑1 (MQ=255)
                                     ggATGAGGCCGTCAAGGGGCTTGGCGCAGACGCCGATCACCGCCCCGAATGCGGCTGGCATTTGGCGGATCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATgt                                                                                                          >  2:59770/1‑149 (MQ=255)
                                      gATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAc                                                                                                         >  2:2059437/1‑150 (MQ=255)
                                        tGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCt                                                                                                                >  2:527857/1‑141 (MQ=255)
                                        tGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCt                                                                                                                <  1:527857/141‑1 (MQ=255)
                                          aGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCTCAAGTCAACATGTCGGTTGCCCCGGGCGGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAcgccg                                                                                                     >  1:690468/1‑150 (MQ=255)
                                           ggCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAcgccgc                                                                                                    <  2:1906550/150‑1 (MQ=255)
                                               gTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                <  2:1901290/150‑1 (MQ=255)
                                               gTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                <  1:2059437/150‑1 (MQ=255)
                                               gTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                >  2:65007/1‑150 (MQ=255)
                                                tCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAcg                                                                                                        <  2:226821/141‑1 (MQ=255)
                                                tCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGAcg                                                                                                        >  1:226821/1‑141 (MQ=255)
                                                  aaGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAg                                                                                             >  2:1386310/1‑150 (MQ=255)
                                                   aGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAgg                                                                                            <  2:1350417/150‑1 (MQ=255)
                                                     gggCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTg                                                                                          >  1:193642/1‑150 (MQ=255)
                                                     gggCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTg                                                                                          >  1:193265/1‑150 (MQ=255)
                                                        cTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAg                                                                                       <  2:2949249/150‑1 (MQ=255)
                                                         ttGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAgg                                                                                      <  1:2581593/150‑1 (MQ=255)
                                                          tGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGGACCGGGCAGTCGTCAGTGGTTGCTTTTGTCGCCACCGATCGCTCTGGATGACGCCGCCTACCAGGTGCAggg                                                                                     <  1:59770/150‑1 (MQ=255)
                                                                    gccgATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGaa                                                                           >  2:1956838/1‑150 (MQ=255)
                                                                      cgATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACg                                                                         >  2:1297349/1‑150 (MQ=255)
                                                                        aTCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCa                                                                       >  2:956318/1‑150 (MQ=255)
                                                                        aTCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCa                                                                       >  2:281423/1‑150 (MQ=255)
                                                                         tCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCAt                                                                      <  1:1297349/150‑1 (MQ=255)
                                                                            ccGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGa                                                                   >  2:669420/1‑150 (MQ=255)
                                                                             cGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgtt                                                                                                                                  >  2:1319592/1‑86 (MQ=255)
                                                                             cGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgtt                                                                                                                                  <  1:1319592/86‑1 (MQ=255)
                                                                              gCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                >  1:2874212/1‑119 (MQ=255)
                                                                              gCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                <  2:2874212/119‑1 (MQ=255)
                                                                                cccGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCa                                                               >  1:943347/1‑150 (MQ=255)
                                                                                     aTGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtg                                                          >  1:769362/1‑150 (MQ=255)
                                                                                         gcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtgg                                                      <  1:956318/150‑1 (MQ=255)
                                                                                         gcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtgg                                                      <  1:1386310/150‑1 (MQ=255)
                                                                                          ccTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTggg                                                     >  1:1637587/1‑150 (MQ=255)
                                                                                           cTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGa                                                    <  2:1079634/150‑1 (MQ=255)
                                                                                            tGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGAt                                                   >  1:1792068/1‑150 (MQ=255)
                                                                                            tGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGAt                                                   <  2:1242328/150‑1 (MQ=255)
                                                                                              gCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATgg                                                 >  2:842153/1‑150 (MQ=255)
                                                                                                atgtGGCGGAGCAGGGTTCCCAAGGCACGGTACCCGCGGCAGGGAACAGCCGCCCGCAGTGGCTGTTCCCGGCCCCCGTCGCTCTGGATGACGCCCTCTACCAGGTGTAGCGGTAGGCCGAACGCATCGAATCAGGTTGGTGGGATggcg                                               <  2:246179/147‑1 (MQ=255)
                                                                                                aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAggg                                                                                     >  1:374913/1‑112 (MQ=255)
                                                                                                aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAggg                                                                                     <  2:374913/112‑1 (MQ=255)
                                                                                                aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcg                                               >  1:2153418/1‑150 (MQ=255)
                                                                                                     gCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAt                                          <  1:1956838/150‑1 (MQ=255)
                                                                                                     gCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAt                                          <  1:784825/150‑1 (MQ=255)
                                                                                                       ggatcAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtgt                                        <  1:266791/146‑1 (MQ=255)
                                                                                                         ggCAGGGTGTCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGc                                      <  2:324935/149‑1 (MQ=255)
                                                                                                             gggTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGACGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTcg                                  <  2:769362/150‑1 (MQ=255)
                                                                                                              ggTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTcgc                                 <  2:1795727/150‑1 (MQ=255)
                                                                                                              ggTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTcgc                                 <  2:1637587/150‑1 (MQ=255)
                                                                                                                   ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGacta                            <  2:2036651/150‑1 (MQ=255)
                                                                                                                   ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGacta                            >  1:390177/1‑150 (MQ=255)
                                                                                                                    cAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGTAGCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGactac                           >  1:1753781/1‑150 (MQ=255)
                                                                                                                       gggaaCATGCTGGTGGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACtactac                        <  2:1753781/147‑1 (MQ=255)
                                                                                                                             aTGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATg                                                                                >  1:2706765/1‑88 (MQ=255)
                                                                                                                             aTGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATg                                                                                <  2:2706765/88‑1 (MQ=255)
                                                                                                                             aTGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAt                                          >  1:813755/1‑126 (MQ=255)
                                                                                                                             aTGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAt                                          <  2:813755/126‑1 (MQ=255)
                                                                                                                                      gCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGt         <  1:541359/150‑1 (MQ=255)
                                                                                                                                        cccGGGCGGCTGCCCGGTGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAAc                                                                          >  2:1049737/2‑83 (MQ=255)
                                                                                                                                        cccGGGCGGCTGCCCGGTGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAAc                                                                          <  1:1049737/82‑1 (MQ=255)
                                                                                                                                        cccGGGCGGCTGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa       >  1:172585/2‑150 (MQ=255)
                                                                                                                                        aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa       >  2:347041/1‑150 (MQ=255)
                                                                                                                                        aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACAGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGACGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGtca       >  1:565233/1‑148 (MQ=255)
                                                                                                                                         ccGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGCTCGCGCGGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGAATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGAc      >  2:1446479/1‑150 (MQ=255)
                                                                                                                                          cGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACg     >  1:1342464/1‑150 (MQ=255)
                                                                                                                                          cGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACg     >  1:1058623/1‑150 (MQ=255)
                                                                                                                                             gCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGGcc  <  2:402351/150‑1 (MQ=255)
                                                                                                                                                 ccgccCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGac                              <  2:2444258/118‑1 (MQ=255)
                                                                                                                                                 ccgccCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGac                              >  1:2444258/1‑118 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGGCC  >  NC_002947/3527893‑3528183

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