Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,413,439 T→A 5.3% N155I (AAC→ATC) ‡ PP_1235 ← AhpC/TSA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,413,4390TA 5.3% 232.5 / 3.9 76N155I (AAC→ATC) ‡PP_1235AhpC/TSA family protein
Reads supporting (aligned to +/- strand):  ref base T (40/32);  new base A (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

TTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGCTTG  >  NC_002947/1413295‑1413588
                                                                                                                                                |                                                                                                                                                     
ttttGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCCATTTCAAGCCTTGTTCAAg                                                                                                                                                  <  2:1513784/150‑1 (MQ=255)
     ggggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTg                                                                                                                                             >  1:1117257/1‑150 (MQ=255)
     ggggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTg                                                                                                                                             >  1:1215510/1‑150 (MQ=255)
       ggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATCGAAGTTGTATGCAGCACTTGTTGTTTGCAAGCCTTGTTCAAGGCTTggg                                                                                                                                           >  1:1374641/1‑150 (MQ=255)
        gATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggc                                                                                                                                          >  1:979713/1‑150 (MQ=255)
        gATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGGGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGCCAGCCCCATTCAATTCCAAGCCTTGAACAAGGCTTgaaa                                                                                                                                          <  2:62895/150‑4 (MQ=255)
         atatTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcg                                                                                                                                         <  1:665867/150‑1 (MQ=255)
         atatTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcg                                                                                                                                         <  1:1625376/150‑1 (MQ=255)
            ttcttcGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgcc                                                                                                                                      >  2:257020/1‑150 (MQ=255)
              cttcGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgccgc                                                                                                                                    <  1:1662602/150‑1 (MQ=255)
               ttcGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAGTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgccgcc                                                                                                                                   >  1:398685/1‑150 (MQ=255)
                    cATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGCGGCCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATAc                                                                                                                              >  2:1698763/1‑150 (MQ=255)
                          gCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGACTAGCGTTTCACCGATGCTTTGAATGTAATTGGTGCTGCGAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                                                        <  2:481539/150‑1 (MQ=255)
                              gCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCCCCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCcaca                                                                                                                    <  2:509498/150‑1 (MQ=255)
                                 aCGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATCCGGCATCCACATGc                                                                                                                 <  1:457185/150‑1 (MQ=255)
                                     tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                                             >  1:638770/1‑150 (MQ=255)
                                      gATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGTAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTCCAAGCCTTGAACAAGGCTTgaaatt                                                                                                                                        >  2:225558/1‑117 (MQ=255)
                                      gATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGTAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTCCAAGCCTTGAACAAGGCTTgaaatt                                                                                                                                        <  1:225558/122‑6 (MQ=255)
                                          aaGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACt                                                                                                        >  1:646901/1‑150 (MQ=255)
                                           aGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACtt                                                                                                       >  1:1328080/1‑150 (MQ=255)
                                            ggCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGACGCTGGATTTAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACTGCATCCACATGCCCGGGCACttt                                                                                                      >  1:798039/1‑150 (MQ=255)
                                              ccGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTca                                                                                                    <  1:841494/150‑1 (MQ=255)
                                               cGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAACGGTTTCGCATGCAATTTCCAGCCATGTTCAAGTCTTGGGCGGCCTCCAATACGGCATCCACATGCCCGGGCCCTTTcac                                                                                                   >  1:425595/1‑150 (MQ=255)
                                                   cTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcacgcca                                                                                               >  2:210488/1‑150 (MQ=255)
                                                    tGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcacgccac                                                                                              <  1:107390/150‑1 (MQ=255)
                                                        tGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAACTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCTGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCAcgccacgcca                                                                                          <  2:979713/150‑1 (MQ=255)
                                                              gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTg                                                                                    >  1:300082/1‑150 (MQ=255)
                                                              gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTg                                                                                    >  2:1080186/1‑150 (MQ=255)
                                                              gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACACCACGCCATTCCTg                                                                                    <  2:1215510/150‑1 (MQ=255)
                                                                ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                                  >  1:211686/1‑150 (MQ=255)
                                                                ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                                  >  2:1867764/1‑150 (MQ=255)
                                                                 cAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgcg                                                                                 >  1:917761/1‑150 (MQ=255)
                                                                        gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCggg                                                                                                            >  2:1723575/1‑116 (MQ=255)
                                                                        gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCggg                                                                                                            <  1:1723575/116‑1 (MQ=255)
                                                                        gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcac                                                                                                   <  1:872316/125‑1 (MQ=255)
                                                                        gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTaac                                                                                                   >  2:872316/1‑125 (MQ=255)
                                                                             tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTc                                                                     <  2:300082/150‑1 (MQ=255)
                                                                                  tttCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                                                                <  1:257020/150‑1 (MQ=255)
                                                                                  tttCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                                                                <  1:1414906/150‑1 (MQ=255)
                                                                                  tttCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                                                                >  2:1285884/1‑150 (MQ=255)
                                                                                      aCCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGac                                                            >  2:1923304/1‑149 (MQ=255)
                                                                                        cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCa                                                          <  2:1374641/150‑1 (MQ=255)
                                                                                          aTGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAgg                                                        >  2:733425/1‑150 (MQ=255)
                                                                                             cTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAg                                                     >  1:1757884/1‑150 (MQ=255)
                                                                                              tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCca                                                                                                                      >  1:808206/1‑84 (MQ=255)
                                                                                              tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCca                                                                                                                      <  2:808206/84‑1 (MQ=255)
                                                                                              tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAgg                                                    >  2:1230852/1‑150 (MQ=255)
                                                                                              tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAgg                                                    >  1:1238918/1‑150 (MQ=255)
                                                                                                 gAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgc                                                 >  2:1540103/1‑150 (MQ=255)
                                                                                                  aaTGTAATTGGGGCTGCAAGGCAGTCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgct                                                >  2:219889/1‑150 (MQ=255)
                                                                                                     gTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCATTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAGCCCGTCCTTGTCGATCAGGAAGGTGCTGCg                                             <  2:1117257/150‑1 (MQ=255)
                                                                                                       aaTTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGt                                           <  2:371133/150‑1 (MQ=255)
                                                                                                                tGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCa                                  >  1:198571/1‑150 (MQ=255)
                                                                                                                tGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCa                                  >  1:1347068/1‑150 (MQ=255)
                                                                                                                tGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTCGTCAGCGAGGAAGGTGCTGCGGTCAACGCCCa                                  >  1:1433961/1‑150 (MQ=255)
                                                                                                                tGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCACGGTCACTTTCACGCCACGCCATTCCTGGCGCATCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCa                                  >  1:1376294/1‑150 (MQ=255)
                                                                                                                     ggCNGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTAc                             >  2:1193485/1‑150 (MQ=255)
                                                                                                                     ggCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGCTGTTGCTGCGGTCAACGCCCATGTAc                             <  1:1698763/150‑1 (MQ=255)
                                                                                                                     ggCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTAc                             <  1:1769726/150‑1 (MQ=255)
                                                                                                                       gaGCCTCATTCAATTTCAAGCCTTGTTCACGGATTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTc                           <  2:398685/149‑1 (MQ=255)
                                                                                                                          ccccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCtt                        <  2:198571/150‑1 (MQ=255)
                                                                                                                            ccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGc                      >  2:1121768/1‑150 (MQ=255)
                                                                                                                            ccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGc                      <  2:1757884/150‑1 (MQ=255)
                                                                                                                            acATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGc                      <  1:435807/149‑1 (MQ=255)
                                                                                                                             cATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGcc                     <  2:211686/150‑1 (MQ=255)
                                                                                                                               ttCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCAt                   <  1:210488/150‑1 (MQ=255)
                                                                                                                                 cAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATAc                 <  2:1238918/150‑1 (MQ=255)
                                                                                                                                    gcccAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGATCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGc              >  1:1842544/4‑150 (MQ=255)
                                                                                                                                       cAAGCCTTGTTCAAGGCTTGGGCGGTCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCCTGTACTCCTTGCCATACAGcttc           <  2:952201/150‑1 (MQ=255)
                                                                                                                                           ccTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCa       >  1:1456256/1‑150 (MQ=255)
                                                                                                                                            cTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAg      >  2:1540957/1‑150 (MQ=255)
                                                                                                                                            cTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAg      >  1:1609850/1‑150 (MQ=255)
                                                                                                                                            cTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGACCATTGCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAg      <  1:1367125/150‑1 (MQ=255)
                                                                                                                                             ttGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGc     <  1:1230852/150‑1 (MQ=255)
                                                                                                                                             ttGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGCCGGCGCGCCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGc     >  1:1562741/1‑150 (MQ=255)
                                                                                                                                              tGTTCAAGGCTTGGGCGGCCGCCAATACGGGATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGGCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGCt    >  2:1119804/1‑150 (MQ=255)
                                                                                                                                                ttCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGCTTg  >  1:190691/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGCTTG  >  NC_002947/1413295‑1413588

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