Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,429,887 C→G 5.4% intergenic (+483/‑29) PP_3039 → / → PP_3040 hypothetical protein/holin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,429,8870CG 5.4% 1202.0 / 42.6 350intergenic (+483/‑29)PP_3039/PP_3040hypothetical protein/holin
Reads supporting (aligned to +/- strand):  ref base C (168/163);  new base G (10/9);  total (178/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 364 total aligned reads displayed.

AACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTGG  >  NC_002947/3429738‑3430035
                                                                                                                                                     |                                                                                                                                                    
aaCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGgc                                                                                                                                                      >  2:3628450/1‑150 (MQ=255)
aaCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGgc                                                                                                                                                      >  1:5551471/1‑150 (MQ=255)
    gACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCg                                                                                                                                                  >  2:7437828/1‑150 (MQ=255)
    gACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCg                                                                                                                                                  <  1:4053183/150‑1 (MQ=255)
    gACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATATTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAAACCCGGCCCTTGCGCCg                                                                                                                                                  >  1:3241867/1‑150 (MQ=255)
     aCATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAAACCCGGCCCTGGCGCCgg                                                                                                                                                 >  2:308865/1‑150 (MQ=255)
        ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGt                                                                                                                                              <  2:5781775/150‑1 (MQ=255)
        ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGt                                                                                                                                              <  2:890678/150‑1 (MQ=255)
        ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACATGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAACGCGGCCCGCCAAACACTGAAAACCCCGGCACTGGCGCCGGGGt                                                                                                                                              <  1:2959118/150‑1 (MQ=255)
         tCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGACACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGtt                                                                                                                                             <  1:1538634/150‑1 (MQ=255)
            ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                          <  1:7250879/150‑1 (MQ=255)
            ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                          <  1:65654/150‑1 (MQ=255)
             cGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttca                                                                                                                                         <  2:800208/150‑3 (MQ=255)
             cGGACGTTTACCACTTAGTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttca                                                                                                                                         <  1:4596964/150‑3 (MQ=255)
                 cGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCtt                                                                                                                                     <  2:5881571/150‑1 (MQ=255)
                   tttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTg                                                                                                                                   <  2:4700517/150‑1 (MQ=255)
                    ttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                  <  1:6894816/150‑1 (MQ=255)
                    ttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                  <  1:1795353/150‑1 (MQ=255)
                     tACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                                 <  2:6276524/150‑1 (MQ=255)
                     tACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                                 <  2:524888/150‑1 (MQ=255)
                      aCCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTGTgggg                                                                                                                                >  2:1965752/1‑150 (MQ=255)
                       ccACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCAGGCGCCGGGGTTTTTGCCTTTggggg                                                                                                                               >  2:6642746/1‑150 (MQ=255)
                         aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                             >  2:7322580/1‑150 (MQ=255)
                         aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                             >  2:2692817/1‑150 (MQ=255)
                             aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                         <  2:1467229/150‑1 (MQ=255)
                             aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                         <  1:1079797/150‑1 (MQ=255)
                              ttGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                        <  2:6751228/150‑1 (MQ=255)
                                gCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGa                                                                                                                      >  2:4994111/1‑150 (MQ=255)
                                 ccGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGaa                                                                                                                     <  1:3387872/150‑1 (MQ=255)
                                  cGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAAc                                                                                                                    >  1:4662478/1‑150 (MQ=255)
                                      tGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttcag                                                                                                                                        <  2:5202462/126‑4 (MQ=39)
                                      tGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttcag                                                                                                                                        >  1:5202462/1‑123 (MQ=39)
                                       ggCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTGTTTGCCTTTGGGGGATTCGATGAACAGCGa                                                                                                               >  1:127524/1‑150 (MQ=255)
                                         cATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagc                                                                                                             <  1:1562681/150‑1 (MQ=255)
                                         cATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagc                                                                                                             >  2:5508498/1‑150 (MQ=255)
                                          aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                            >  2:201719/1‑150 (MQ=255)
                                          aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                            <  1:3759510/150‑1 (MQ=255)
                                          aTTAAAACTCGCGTGCTGCCACCGGAATCGACCTGTAAATAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                            <  2:6037543/150‑1 (MQ=255)
                                            tAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCAATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAAACCCGGCCCTGGCGCCGGGGTTTTTGACTTTGGGGGATTCGATGAACAGCGagcagc                                                                                                          >  1:397628/1‑150 (MQ=255)
                                                  tCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGt                                                                                                    >  1:5945155/1‑150 (MQ=255)
                                                    gCTTGCTGCCACTGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTa                                                                                                  <  2:6959132/150‑1 (MQ=255)
                                                      ttGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                  >  2:2437170/1‑116 (MQ=255)
                                                         ctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTga                                                                                             <  1:2738519/150‑1 (MQ=255)
                                                           gcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                           >  2:2627280/1‑150 (MQ=255)
                                                           gcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                           <  1:4811030/150‑1 (MQ=255)
                                                             cACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATg                                                                                         >  2:2341005/1‑150 (MQ=255)
                                                                 ggAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCaa                                                                                     <  1:211030/150‑1 (MQ=255)
                                                                   aaTCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGAGGCCCCCCAAACACTTAAAAAACCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAACCCTTAGTTGAGATTCCaaac                                                                                   >  1:2481793/1‑148 (MQ=255)
                                                                    aTCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCt                                                                                  >  2:3911873/1‑150 (MQ=255)
                                                                      cGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtgg                                                                                <  1:3565078/150‑1 (MQ=255)
                                                                        aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                              <  1:2253341/150‑1 (MQ=255)
                                                                        aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                              <  1:5345303/150‑1 (MQ=255)
                                                                        aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                              <  1:6749073/150‑1 (MQ=255)
                                                                        aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                              >  1:1144676/1‑150 (MQ=255)
                                                                          cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttgg                                                                            >  1:6347339/1‑150 (MQ=255)
                                                                                  aGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                                    <  1:1959711/150‑1 (MQ=255)
                                                                                       ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTc                                                               >  2:5658977/1‑150 (MQ=255)
                                                                                        cATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTcc                                                              >  1:5674820/1‑150 (MQ=255)
                                                                                         aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTccc                                                             >  2:4509978/1‑150 (MQ=255)
                                                                                          tCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCt                                                            >  1:4186246/1‑150 (MQ=255)
                                                                                          tCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCt                                                            <  2:325056/150‑1 (MQ=255)
                                                                                          tCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCt                                                            >  2:4120878/1‑150 (MQ=255)
                                                                                            ttCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTgg                                                          >  1:1031275/1‑150 (MQ=255)
                                                                                                agaCGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTcgg                                                      >  1:4864777/1‑150 (MQ=255)
                                                                                                 gaCGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTcggc                                                     <  1:3075286/150‑1 (MQ=255)
                                                                                                   cGTGCGGGCGCACAAAGCGGCCCACCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTcggcgg                                                   <  1:7009868/150‑1 (MQ=255)
                                                                                                     tGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCgggg                                                 <  1:266350/150‑1 (MQ=255)
                                                                                                      gCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGgt                                                <  2:3151989/150‑1 (MQ=255)
                                                                                                      gCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGgt                                                <  2:64839/150‑1 (MQ=255)
                                                                                                         ggCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAAGCTGGGTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTc                                             >  2:5570783/1‑150 (MQ=255)
                                                                                                          gcgcACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCa                                            >  1:608597/1‑150 (MQ=255)
                                                                                                                aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttgg                                                                            >  1:1453713/1‑112 (MQ=255)
                                                                                                                aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCgag                                      >  2:1388207/1‑150 (MQ=255)
                                                                                                                 aaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCgaga                                     <  2:3839412/150‑1 (MQ=255)
                                                                                                                     ggCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGAtgtg                                 <  2:3420099/150‑1 (MQ=255)
                                                                                                                              aaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                        >  1:1971920/1‑54 (MQ=255)
                                                                                                                              aaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGa                        <  1:877162/150‑1 (MQ=255)
                                                                                                                                   tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg                   >  1:6438144/1‑150 (MQ=255)
                                                                                                                                   tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg                   >  2:1172300/1‑150 (MQ=255)
                                                                                                                                   tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCGTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg                   >  1:4393885/1‑150 (MQ=255)
                                                                                                                                    aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTGGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                  >  2:2808557/2‑150 (MQ=255)
                                                                                                                                    aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                  >  1:2297622/2‑150 (MQ=255)
                                                                                                                                    aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                  >  1:3887363/2‑150 (MQ=255)
                                                                                                                                    aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                  >  1:5415987/2‑150 (MQ=255)
                                                                                                                                    aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                  >  1:6627435/2‑150 (MQ=255)
                                                                                                                                     aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCCGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                                  <  2:5085074/117‑1 (MQ=255)
                                                                                                                                     aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGcc                 >  1:7279595/1‑150 (MQ=255)
                                                                                                                                     aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGcc                 >  1:7116855/1‑150 (MQ=255)
                                                                                                                                      aaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCg                >  2:2882997/1‑150 (MQ=255)
                                                                                                                                        cccccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc              <  2:6857814/149‑1 (MQ=255)
                                                                                                                                        aCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc              <  2:5381075/150‑1 (MQ=255)
                                                                                                                                          cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt            >  1:3840713/1‑150 (MQ=255)
                                                                                                                                          cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt            <  2:1609669/150‑1 (MQ=255)
                                                                                                                                          cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt            >  1:6240968/1‑150 (MQ=255)
                                                                                                                                          cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt            >  2:3718044/1‑150 (MQ=255)
                                                                                                                                           ccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGGGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGtt           >  1:3288774/1‑150 (MQ=255)
                                                                                                                                            cGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTa          <  2:4854844/150‑1 (MQ=255)
                                                                                                                                                   ggCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTg   >  2:6179632/1‑150 (MQ=255)
                                                                                                                                                    gcgcCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTgg  <  2:1463739/150‑1 (MQ=255)
                                                                                                                                                    gcgcCGGGGTTTTTGCCTTGGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTgg  <  1:2267155/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTGG  >  NC_002947/3429738‑3430035

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