Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,483,550 A→T 5.6% intergenic (‑88/‑120) PP_3975 ← / → PP_3976 hypothetical protein/isochorismatase family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,483,5500AT 5.6% 279.6 / 4.5 89intergenic (‑88/‑120)PP_3975/PP_3976hypothetical protein/isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (45/39);  new base T (3/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

CCTCTTCCCATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTCCTGA  >  NC_002947/4483402‑4483697
                                                                                                                                                    |                                                                                                                                                   
ccTCTTCCCATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGTTTGGACGCTACCAATGGCGCATTaaa                                                                                                                                                     <  7:220019/149‑1 (MQ=255)
 ctctTCCCATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGTTTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAg                                                                                                                                                    >  8:279474/1‑149 (MQ=255)
 ctctTCCCATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAg                                                                                                                                                    >  8:59648/1‑149 (MQ=255)
 ctctTCCCATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCATAGCACCGGATTGGACGCTACCAATGGCGCATTAAAg                                                                                                                                                    >  8:55153/1‑149 (MQ=255)
       ccATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTc                                                                                                                                              >  4:198833/1‑149 (MQ=255)
        cATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTc                                                                                                                                              >  5:12672/1‑148 (MQ=255)
         aTGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCTTTAATGCGccattg                                                                                                                                            >  5:190268/1‑145 (MQ=255)
         aTGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTcc                                                                                                                                             >  7:154088/1‑148 (MQ=255)
           gAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACTTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCaa                                                                                                                                           <  6:224969/148‑1 (MQ=255)
           gAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGAGCATTAAAGCGCCTcc                                                                                                                                             >  5:112011/1‑146 (MQ=255)
            aaaCGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCTTTAATGCGccattggt                                                                                                                                          <  6:190268/148‑7 (MQ=255)
              aCGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAgc                                                                                                                                       >  4:65804/1‑149 (MQ=255)
              aCGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAgc                                                                                                                                       >  3:115208/1‑149 (MQ=255)
                 ccccAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAgc                                                                                                                                                   >  4:245641/1‑134 (MQ=255)
                 ccccAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAgc                                                                                                                                                   <  3:245641/134‑1 (MQ=255)
                  cccAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTCGATTCTCCTTGGCGAAAAATGGATCCACGTGAACCGACTCGCCGAAATCCTAGCACCGGATTGGACGCAACCAATGGCGCATTAAAGCGCCACCAAGAGCGCa                                                                                                                                    >  3:62783/1‑148 (MQ=255)
                   ccAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCaa                                                                                                                                   >  6:330650/1‑148 (MQ=255)
                       tAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTg                                                                                                                              >  1:226228/1‑149 (MQ=255)
                         aGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTgcg                                                                                                                            <  7:262971/149‑1 (MQ=255)
                         aGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTgc                                                                                                                             <  7:113492/148‑1 (MQ=255)
                         aGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTgc                                                                                                                             <  7:279474/148‑1 (MQ=255)
                                   tcgtTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTc                                                                                                                  >  1:301326/1‑149 (MQ=255)
                                      ttGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCt                                                                                                               >  3:324382/1‑149 (MQ=255)
                                         aaaGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGa                                                                                                             <  1:337641/148‑1 (MQ=255)
                                          aaGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAgt                                                                                                           <  3:276712/149‑1 (MQ=255)
                                             gTCGTGATGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAgc                                                                                                                                       <  6:56814/118‑1 (MQ=255)
                                              tCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGttt                                                                                                       <  3:289311/149‑1 (MQ=255)
                                                gTGAAGGGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCTCCGGTTTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCGACCGATTCGCTGAGTGTTTc                                                                                                      <  4:62783/148‑1 (MQ=255)
                                                gTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTc                                                                                                      <  3:198833/148‑1 (MQ=255)
                                                gTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCg                                                                                                     >  2:107665/1‑149 (MQ=255)
                                                 tGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGc                                                                                                    <  5:330650/149‑1 (MQ=255)
                                                     gCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTgc                                                                                                                             <  2:3107/120‑1 (MQ=255)
                                                     gCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTgc                                                                                                                             >  1:3107/1‑120 (MQ=255)
                                                     gCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATc                                                                                                <  8:114288/149‑1 (MQ=255)
                                                              ggATTCTCCTTCGCGAAAAATGGACCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTg                                                                                        <  6:16942/148‑1 (MQ=255)
                                                               gATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTgg                                                                                       <  1:67133/148‑1 (MQ=255)
                                                               gATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGc                                                                                      <  7:59648/149‑1 (MQ=255)
                                                                aTTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGc                                                                                      <  7:55153/148‑1 (MQ=255)
                                                                       ttGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGc                                                                               >  7:258561/1‑148 (MQ=255)
                                                                       ttGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGc                                                                               >  4:118408/1‑148 (MQ=255)
                                                                            gAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGa                                                                                                             >  4:36339/1‑113 (MQ=255)
                                                                            gAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGa                                                                                                             <  3:36339/113‑1 (MQ=255)
                                                                            gAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCgacgac                                                                         <  6:112011/149‑1 (MQ=255)
                                                                            gAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCgacgac                                                                         >  4:271236/1‑149 (MQ=255)
                                                                            gAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCgacgac                                                                         >  2:114316/1‑149 (MQ=255)
                                                                            gAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGCCTGACCCTCGGcgccatc                                                                         <  1:107665/149‑6 (MQ=255)
                                                                               aaaTGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACccg                                                                      >  5:338600/1‑149 (MQ=255)
                                                                                aaTGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACccgc                                                                     <  3:271236/149‑1 (MQ=255)
                                                                                             gATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCt                                                                                                               <  2:289902/94‑1 (MQ=255)
                                                                                             gATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCt                                                                                                               >  1:289902/1‑94 (MQ=255)
                                                                                              aTCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACccgc                                                                     >  5:283169/1‑135 (MQ=255)
                                                                                              aTCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACccgc                                                                     <  6:283169/135‑1 (MQ=255)
                                                                                              aTCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGc                                                       >  5:160057/1‑149 (MQ=255)
                                                                                                    cTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGt                                                 >  2:156458/1‑149 (MQ=255)
                                                                                                      ggCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGtt                                                >  3:199655/1‑148 (MQ=255)
                                                                                                       gCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGtt                                                >  1:136800/1‑147 (MQ=255)
                                                                                                       gCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCa                                              <  3:65804/149‑1 (MQ=255)
                                                                                                         cGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGc                                            <  8:258561/149‑1 (MQ=255)
                                                                                                          gAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGc                                            >  8:148818/1‑148 (MQ=255)
                                                                                                              tCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAg                                       >  3:71029/1‑149 (MQ=255)
                                                                                                               ccTAGCACCGGATTGGACGCAACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAgg                                      >  6:187886/1‑149 (MQ=255)
                                                                                                                cTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTgg                                                                                            >  2:1932/1‑94 (MQ=255)
                                                                                                                cTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTgg                                                                                            <  1:1932/94‑1 (MQ=255)
                                                                                                                 tAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAggg                                     >  1:13149/1‑148 (MQ=255)
                                                                                                                  aGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAggg                                     >  7:17676/1‑147 (MQ=255)
                                                                                                                         gATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAAt                            >  4:56775/1‑149 (MQ=255)
                                                                                                                           ttGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGa                          <  8:17676/149‑1 (MQ=255)
                                                                                                                              gACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGc                        >  7:279773/1‑148 (MQ=255)
                                                                                                                               aCGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGa                          >  8:219134/1‑145 (MQ=255)
                                                                                                                               aCGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGa                          <  7:219134/145‑1 (MQ=255)
                                                                                                                                 gCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCgacgac                                                                         <  2:224666/96‑1 (MQ=255)
                                                                                                                                 gCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCgacgac                                                                         >  1:224666/1‑96 (MQ=255)
                                                                                                                                 gCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGACACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATc                    <  1:134161/149‑1 (MQ=255)
                                                                                                                                   tACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGACACTCGGCGACGACACGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATcc                   >  4:46821/1‑148 (MQ=255)
                                                                                                                                   tACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGATGACACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCca                  >  4:51116/1‑149 (MQ=255)
                                                                                                                                    aCCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCca                  >  5:119498/1‑148 (MQ=255)
                                                                                                                                     ccAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCcaca                <  2:13149/149‑1 (MQ=255)
                                                                                                                                     ccAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCcac                 <  3:118408/148‑1 (MQ=255)
                                                                                                                                      cAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCa                                              >  1:157696/1‑118 (MQ=255)
                                                                                                                                      cAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCa                                              <  2:157696/118‑1 (MQ=255)
                                                                                                                                       ggaggCGCTTTAATGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCcacaca              >  1:294205/4‑149 (MQ=255)
                                                                                                                                       ggaggCGCTTTAATGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCcacaca              >  4:173251/4‑149 (MQ=255)
                                                                                                                                       aaTGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCcacaca              >  7:152745/1‑149 (MQ=255)
                                                                                                                                            cgcATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGaa         <  3:56775/149‑1 (MQ=255)
                                                                                                                                             gctttAATGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAg        <  3:173251/146‑1 (MQ=255)
                                                                                                                                             gcATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAg        <  4:115208/149‑1 (MQ=255)
                                                                                                                                               aTTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTc      >  4:122533/1‑149 (MQ=255)
                                                                                                                                               aTTAAAGCGCCTCCAAGAGCGAAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGAAGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTc      <  5:187886/149‑1 (MQ=255)
                                                                                                                                                ttAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTcc     <  4:71029/149‑1 (MQ=255)
                                                                                                                                                ttAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTcc     <  2:301326/149‑1 (MQ=255)
                                                                                                                                                  aaaGCTCCTCCAAGATCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGAAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTCCt    <  8:154088/148‑1 (MQ=255)
                                                                                                                                                  aaaGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTCCt    <  4:324382/148‑1 (MQ=255)
                                                                                                                                                   aaGCGCCTCCAAGAGCGCAAAATTGCGCCCCCGATTCGCTGACTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTCCTGa  >  8:192207/1‑149 (MQ=255)
                                                                                                                                                    aGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAggg                                     >  8:13721/1‑113 (MQ=255)
                                                                                                                                                    aGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAggg                                     <  7:13721/113‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTCTTCCCATGAAACGCCCCAGTAAGCAGTGTCGTCGTTGAAAGGTCGTGAAGCGGCCATTGGATTCTCCTTGGCGAAAAATGGATCGACGTGATCCGACTGGCCGAAATCCTAGCACCGGATTGGACGCTACCAATGGCGCATTAAAGCGCCTCCAAGAGCGCAAAATTGCGCCACCGATTCGCTGAGTGTTTCGCTATCTTGGGCTGGCACTCGGCGACGACCCGCCGAAAAACTCTCGCCACCGTTCAGCAAAAGGGAAGAACAATGAGCAATCCACACACTGAAGTCCTGA  >  NC_002947/4483402‑4483697

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