Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,488,389 G→A 7.6% R88H (CGC→CAC)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,488,3890GA 7.6% 1045.9 / 48.7 342R88H (CGC→CAC) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base G (146/170);  new base A (5/21);  total (151/191)
Fisher's exact test for biased strand distribution p-value = 7.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 348 total aligned reads displayed.

CAAGCGCTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCA  >  NC_002947/2488240‑2488538
                                                                                                                                                     |                                                                                                                                                     
   gcgcTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCt                                                                                                                                                    <  2:2640487/150‑1 (MQ=255)
    cgcTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTg                                                                                                                                                   <  2:4075050/150‑1 (MQ=255)
       tACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTg                                                                                                                                                >  2:3404310/1‑150 (MQ=255)
        aCCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGc                                                                                                                                               <  1:5321022/150‑1 (MQ=255)
         ccccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTAGTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGcc                                                                                                                                              <  2:6750130/150‑1 (MQ=255)
           ccctGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGACCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            <  1:5681886/147‑1 (MQ=255)
           ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGACCTGCGTGCTGCTGGGTCCCCAGGGCGGCCTTTACCCGCTTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            <  2:4525991/150‑1 (MQ=255)
           ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            <  2:5424568/150‑1 (MQ=255)
           ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            <  1:2266307/150‑1 (MQ=255)
           ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            <  1:2522640/150‑1 (MQ=255)
           ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            <  1:2811759/150‑1 (MQ=255)
           ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGACGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                            >  1:3845994/1‑150 (MQ=255)
            cATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGCCggg                                                                                                                                           <  2:1564947/150‑1 (MQ=255)
              tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGAGCCGGATGCTGCCCGGCACCAGGCGCCgggca                                                                                                                                         <  2:5770002/150‑3 (MQ=255)
              tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCACCAGGCGCCgggca                                                                                                                                         <  2:5311648/150‑3 (MQ=255)
              tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCACCAGGCGCCgggca                                                                                                                                         <  2:1932520/150‑3 (MQ=255)
              tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCGCCTGGTGCCGGGTa                                                                                                                                         <  2:5236555/150‑1 (MQ=255)
                gACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGCCgggcagc                                                                                                                                       <  1:4927840/150‑5 (MQ=255)
                     gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                  <  1:4414666/150‑1 (MQ=255)
                     gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                  >  2:3632780/1‑150 (MQ=255)
                      cTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAg                                                                                                                                 >  1:2717775/1‑150 (MQ=255)
                       tGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                                >  1:203097/1‑150 (MQ=255)
                       tGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                                >  2:307042/1‑150 (MQ=255)
                         gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                              <  1:4879844/150‑1 (MQ=255)
                         gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                              <  2:4554398/150‑1 (MQ=255)
                         gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                              >  1:5951702/1‑150 (MQ=255)
                         gAAAAAGTAGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                              >  2:3872043/1‑150 (MQ=255)
                          aaaaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCgg                                                                                                                             <  2:3724187/150‑1 (MQ=255)
                           aaaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGt                                                                                                                            <  1:3202298/150‑1 (MQ=255)
                           aaaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGt                                                                                                                            >  2:2816809/1‑150 (MQ=255)
                            aaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTg                                                                                                                           >  1:3739599/1‑150 (MQ=255)
                             aaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTgc                                                                                                                          <  2:1703675/150‑1 (MQ=255)
                              aGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTgcg                                                                                                                         >  1:5241503/1‑150 (MQ=255)
                                tCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCt                                                                                                                       >  1:3527402/1‑150 (MQ=255)
                                  gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTgg                                                                                                                     >  2:5518983/1‑150 (MQ=255)
                                  gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTgg                                                                                                                     >  1:6430812/1‑150 (MQ=255)
                                  gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTgg                                                                                                                     >  1:755211/1‑150 (MQ=255)
                                   cAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTggg                                                                                                                    <  1:1200267/150‑1 (MQ=255)
                                           gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                            <  2:4456872/150‑1 (MQ=255)
                                           gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                            <  2:6553128/150‑1 (MQ=255)
                                           gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                            <  1:6216251/150‑1 (MQ=255)
                                             gcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGt                                                                                                          <  2:2352807/150‑1 (MQ=255)
                                             gcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGt                                                                                                          >  2:101676/1‑150 (MQ=255)
                                             gcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTg                                                                                                         >  1:6742065/1‑150 (MQ=255)
                                               ctctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGa                                                                                                        >  2:5812823/1‑150 (MQ=255)
                                                tctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGAt                                                                                                       >  1:5730015/1‑150 (MQ=255)
                                                 ctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATc                                                                                                      <  1:3224814/150‑1 (MQ=255)
                                                 ctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATc                                                                                                      <  2:642727/150‑1 (MQ=255)
                                                           cGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGc                                                                                            >  2:3530107/1‑150 (MQ=255)
                                                              aGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGa                                                                                         <  2:4113186/150‑1 (MQ=255)
                                                                   tGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACg                                                                                    >  2:2745550/1‑150 (MQ=255)
                                                                     cGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTc                                                                                  <  1:1460297/150‑1 (MQ=255)
                                                                       tgctgcTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCg                                                                                >  1:100063/1‑150 (MQ=255)
                                                                        gctgctGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCTGGGAACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGa                                                                               >  2:6647131/1‑150 (MQ=255)
                                                                             tgGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCa                                                                          >  2:796489/1‑150 (MQ=255)
                                                                               ggTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGACGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACt                                                                        <  1:1661112/150‑1 (MQ=255)
                                                                                gTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTg                                                                       >  1:4067093/1‑150 (MQ=255)
                                                                                   cccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCCCTGCaa                                                                    >  2:8444/1‑150 (MQ=255)
                                                                                    ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                   >  1:5931381/1‑150 (MQ=255)
                                                                                      aGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGc                                                                 >  2:5004019/1‑150 (MQ=255)
                                                                                         gcggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGa                                                              <  2:433256/150‑1 (MQ=255)
                                                                                         gcggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGa                                                              <  1:707432/150‑1 (MQ=255)
                                                                                           ggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATc                                                            <  1:1232514/150‑1 (MQ=255)
                                                                                           ggcCTTTACCCGATGGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATc                                                            <  2:3459893/150‑1 (MQ=255)
                                                                                             ccTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGa                                                          >  2:1259303/1‑150 (MQ=255)
                                                                                                ttACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAAtt                                                       <  2:4399350/150‑1 (MQ=255)
                                                                                                 tACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTc                                                      <  1:4718787/150‑1 (MQ=255)
                                                                                                  aCCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCg                                                     <  1:3119350/150‑1 (MQ=255)
                                                                                                   cccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGa                                                    <  2:6314009/150‑1 (MQ=255)
                                                                                                   cccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGa                                                    <  2:6083529/150‑1 (MQ=255)
                                                                                                      gATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGcc                                                 >  1:4556209/1‑150 (MQ=255)
                                                                                                        ttGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTc                                               <  2:1463042/150‑1 (MQ=255)
                                                                                                           gTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTg                                            >  1:3633970/1‑150 (MQ=255)
                                                                                                               cTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGt                                        <  1:5518983/150‑1 (MQ=255)
                                                                                                                tATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTg                                       >  1:241504/1‑150 (MQ=255)
                                                                                                                tATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTg                                       >  1:4186439/1‑150 (MQ=255)
                                                                                                                       tCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgcgc                                <  1:1557433/150‑1 (MQ=255)
                                                                                                                         gCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCgt                              <  2:3178425/150‑1 (MQ=255)
                                                                                                                          ccGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCgtg                             <  2:2117036/150‑1 (MQ=255)
                                                                                                                             aCGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttc                          <  2:4864578/150‑1 (MQ=255)
                                                                                                                             aCGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttc                          >  2:5091768/1‑150 (MQ=255)
                                                                                                                              cGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttct                         <  2:2418328/150‑1 (MQ=255)
                                                                                                                               ggCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttctt                        <  1:3003785/150‑1 (MQ=255)
                                                                                                                                gCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttcttc                       <  1:2470534/150‑1 (MQ=255)
                                                                                                                                 cGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCa                      <  1:4441452/150‑1 (MQ=255)
                                                                                                                                    cccGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCg                   >  2:6787626/1‑150 (MQ=255)
                                                                                                                                          gcgtacccGGCACCAGGCGCCGGGCACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGc             >  1:5022533/6‑150 (MQ=255)
                                                                                                                                           gctgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATc                                                                                       >  1:5206738/1‑75 (MQ=255)
                                                                                                                                           gctgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGcc            <  2:2910794/150‑1 (MQ=255)
                                                                                                                                           gctgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGcc            >  2:5994598/1‑150 (MQ=255)
                                                                                                                                             tgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCaa          >  1:677717/1‑150 (MQ=255)
                                                                                                                                             tacccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCaa          >  1:895448/3‑150 (MQ=255)
                                                                                                                                              gcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTACAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAg         >  1:5077951/1‑150 (MQ=255)
                                                                                                                                              acccGGCACCAGGCGCCGGGCACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAg         >  1:1402299/2‑150 (MQ=255)
                                                                                                                                               cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGc        <  2:272100/150‑1 (MQ=255)
                                                                                                                                                  ggCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTg     <  1:4284902/150‑1 (MQ=255)
                                                                                                                                                  ggCGCCTGGTGCCGGGGACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTg     <  1:1848035/150‑1 (MQ=255)
                                                                                                                                                   gcgcCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTgg    <  2:1315575/150‑1 (MQ=255)
                                                                                                                                                   gcgcCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTgg    <  2:4761214/150‑1 (MQ=255)
                                                                                                                                                     gcCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCa  >  2:2497086/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGCGCTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCA  >  NC_002947/2488240‑2488538

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: