Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,429,884 T→C 5.5% intergenic (+480/‑32) 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,8840TC 5.5% 1093.6 / 32.8 351intergenic (+480/‑32)PP_3039/PP_3040hypothetical protein/holin
Reads supporting (aligned to +/- strand):  ref base T (167/156);  new base C (10/9);  total (185/166)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-01

Only 100 of 358 total aligned reads displayed.

TTCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCT  >  NC_002947/3429735‑3430033
                                                                                                                                                     |                                                                                                                                                     
ttCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCt                                                                                                                                                       >  2:4006373/1‑150 (MQ=255)
ttCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCt                                                                                                                                                       >  2:288187/1‑150 (MQ=255)
    aCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCCAACACTGAAAACCCCGGCCCTGgcg                                                                                                                                                   >  1:6036454/1‑150 (MQ=255)
       gACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCg                                                                                                                                                >  2:7437828/1‑150 (MQ=255)
           ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGt                                                                                                                                            <  2:5781775/150‑1 (MQ=255)
           ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGt                                                                                                                                            <  2:890678/150‑1 (MQ=255)
            tCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGtt                                                                                                                                           <  2:3796857/150‑1 (MQ=255)
            tCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGtt                                                                                                                                           <  1:4119971/150‑1 (MQ=255)
            gccccGGACGATTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGtt                                                                                                                                           <  2:5925750/149‑1 (MQ=255)
              cccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGtttt                                                                                                                                         >  1:5381075/1‑150 (MQ=255)
               ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                        <  1:5759556/150‑1 (MQ=255)
               ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                        <  1:3760601/150‑1 (MQ=255)
                cGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttca                                                                                                                                       <  2:800208/150‑3 (MQ=255)
                cGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTAAAAACCCCGGCGCCAGGGCCGGGGttttca                                                                                                                                       <  2:3038232/150‑3 (MQ=255)
                       ttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                <  1:7267568/150‑1 (MQ=255)
                       ttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                <  1:1795353/150‑1 (MQ=255)
                       ttACCACTTATTGCCGGGTGGCAATAAAACTCGCTGGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                <  2:5752107/150‑1 (MQ=255)
                        tACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                               <  2:6315666/150‑1 (MQ=255)
                         aCCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTGTgggg                                                                                                                              >  2:1965752/1‑150 (MQ=255)
                           cACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGa                                                                                                                            <  2:3893353/150‑1 (MQ=255)
                           cACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGa                                                                                                                            >  2:2672572/1‑150 (MQ=255)
                            aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                           >  2:7322580/1‑150 (MQ=255)
                            aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                           >  2:6718188/1‑150 (MQ=255)
                            aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                           <  1:3186721/150‑1 (MQ=255)
                                aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                       <  1:1079797/150‑1 (MQ=255)
                                aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                       <  2:1467229/150‑1 (MQ=255)
                                aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                       >  1:4392507/1‑150 (MQ=255)
                                 ttGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                      <  2:5584610/150‑1 (MQ=255)
                                 ttGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                      <  2:6751228/150‑1 (MQ=255)
                                 ttGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                      <  1:6663456/150‑1 (MQ=255)
                                 ttGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                      >  1:6293012/1‑150 (MQ=255)
                                    ccGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGaa                                                                                                                   <  2:6190608/150‑1 (MQ=255)
                                         tGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttcag                                                                                                                                      <  2:5202462/126‑4 (MQ=39)
                                          ggCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTGTTGCCTTTGGGGGATTCGATGAACAGCGa                                                                                                             >  2:538159/1‑150 (MQ=255)
                                          ggCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTGTTTGCCTTTGGGGGATTCGATGAACAGCGa                                                                                                             >  1:127524/1‑150 (MQ=255)
                                            cATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagc                                                                                                           <  1:1562681/150‑1 (MQ=255)
                                             aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagcc                                                                                                          >  1:697749/1‑149 (MQ=255)
                                             aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                          <  2:3707000/150‑1 (MQ=255)
                                             aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                          >  2:201719/1‑150 (MQ=255)
                                             aTTAAAACTCGCGTGCTGCCACCGGAATCGACCTGTAAATAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                          <  2:6037543/150‑1 (MQ=255)
                                             aTTAAAACACGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACAGCAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCGTTGGGGGATTCGATGAACAGCGagca                                                                                                          >  1:4652887/1‑150 (MQ=255)
                                                 aaaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCaa                                                                                                      >  1:4042188/1‑150 (MQ=255)
                                                   aCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGc                                                                                                    >  1:2966955/1‑150 (MQ=255)
                                                     tCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGt                                                                                                  >  1:4683467/1‑150 (MQ=255)
                                                         ttGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGt                                                                                              <  1:2502496/150‑1 (MQ=255)
                                                           gctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTg                                                                                            <  2:5291646/150‑1 (MQ=255)
                                                            ctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTga                                                                                           <  1:2738519/150‑1 (MQ=255)
                                                            ctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTga                                                                                           <  1:5439403/150‑1 (MQ=255)
                                                            ctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTga                                                                                           >  2:136935/1‑150 (MQ=255)
                                                              gcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                         <  1:6459815/150‑1 (MQ=255)
                                                              gcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                         >  2:2627280/1‑150 (MQ=255)
                                                                cACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATg                                                                                       >  2:2341005/1‑150 (MQ=255)
                                                                 aCCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGc                                                                                      >  1:2206721/1‑150 (MQ=255)
                                                                 aCCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGATTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGc                                                                                      <  2:3140987/150‑1 (MQ=255)
                                                                   cGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGGTGAGATGCCa                                                                                    >  1:3424275/1‑150 (MQ=255)
                                                                         cGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtgg                                                                              <  1:3565078/150‑1 (MQ=255)
                                                                          gACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggt                                                                             >  1:712818/1‑150 (MQ=255)
                                                                          gACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggt                                                                             <  1:7111749/150‑1 (MQ=255)
                                                                           aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                            >  1:1144676/1‑150 (MQ=255)
                                                                           aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                            <  1:1526294/150‑1 (MQ=255)
                                                                              agTAAATCGTACCAAGCTTGGAGACGTGCGGGCTCTCTTAGCGGCCCGCCAAGCACTGAAAACCCCGGGCCTGGCGCCGGGGTTTTTGGCTTTGGTGGATTAGATGTACAGCTAGCAGCAAGAGTTAGTTGAGATGCCAATCtggttggt                                                                         <  1:4834032/149‑1 (MQ=255)
                                                                               gTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTg                                                                        >  1:3256914/1‑150 (MQ=255)
                                                                                tAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGa                                                                       >  2:2416894/1‑150 (MQ=255)
                                                                                  aaaaGTACCCATCTTCTAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCGTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATc                                                                     >  1:6818412/1‑150 (MQ=255)
                                                                                       tACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCt                                                                <  2:1088091/150‑1 (MQ=255)
                                                                                            aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTccc                                                           >  2:4509978/1‑150 (MQ=255)
                                                                                             tCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCt                                                          >  2:6965483/1‑150 (MQ=255)
                                                                                             tCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCt                                                          >  2:4120878/1‑150 (MQ=255)
                                                                                              cTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTg                                                         <  2:3202962/150‑1 (MQ=255)
                                                                                               ttCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTgg                                                        >  1:1031275/1‑150 (MQ=255)
                                                                                                       gTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                                <  1:1879970/150‑1 (MQ=255)
                                                                                                       gTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                                >  1:7424262/1‑150 (MQ=255)
                                                                                                       gTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                                >  1:6131251/1‑150 (MQ=255)
                                                                                                        tGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCgggg                                               <  1:266350/150‑1 (MQ=255)
                                                                                                         gCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGgt                                              <  2:64839/150‑1 (MQ=255)
                                                                                                                cacaAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCGTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCACTCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGc                                       >  2:5010063/1‑150 (MQ=255)
                                                                                                                   aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttgg                                                                          <  2:1453713/112‑1 (MQ=255)
                                                                                                                   aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttgg                                                                          >  1:1453713/1‑112 (MQ=255)
                                                                                                                       cGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGAtgt                                <  1:6372266/150‑1 (MQ=255)
                                                                                                                        ggCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGAtgtg                               >  1:4410902/1‑150 (MQ=255)
                                                                                                                         gCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTgg                              <  1:6602572/150‑1 (MQ=255)
                                                                                                                              ccAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                                  <  1:3385348/109‑1 (MQ=255)
                                                                                                                                 aaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                      >  1:1971920/1‑54 (MQ=255)
                                                                                                                                      tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGTAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg                 >  2:3460129/1‑150 (MQ=255)
                                                                                                                                      tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg                 >  1:2097566/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)
                                                                                                                                        aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCCGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                                <  2:5085074/117‑1 (MQ=255)
                                                                                                                                        aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                                >  1:5085074/1‑117 (MQ=255)
                                                                                                                                           cccccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc            <  2:6857814/149‑1 (MQ=255)
                                                                                                                                           aCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc            <  2:3386586/150‑1 (MQ=255)
                                                                                                                                             cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt          <  2:1609669/150‑1 (MQ=255)
                                                                                                                                             cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt          >  1:3840713/1‑150 (MQ=255)
                                                                                                                                             cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt          >  1:5616580/1‑150 (MQ=255)
                                                                                                                                             cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt          >  1:2265880/1‑150 (MQ=255)
                                                                                                                                             cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGt          >  2:4729889/1‑150 (MQ=255)
                                                                                                                                               cGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTa        <  2:4854844/150‑1 (MQ=255)
                                                                                                                                                 gCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGc      <  2:5551471/150‑1 (MQ=255)
                                                                                                                                                    cTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGc   <  1:5313094/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCT  >  NC_002947/3429735‑3430033

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: