Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,699,452 C→T 7.4% R36R (CGC→CGT hslU → protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,699,4520CT 7.4% 1194.6 / 56.2 393R36R (CGC→CGThslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (183/180);  new base T (21/8);  total (205/188)
Fisher's exact test for biased strand distribution p-value = 3.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 399 total aligned reads displayed.

GCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTC  >  NC_002947/5699303‑5699599
                                                                                                                                                     |                                                                                                                                                   
  ttttCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCAt                                                                                                                                                   >  2:8841537/1‑150 (MQ=255)
  ttttCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCAt                                                                                                                                                   >  2:1954350/1‑150 (MQ=255)
    ttCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGc                                                                                                                                                 >  1:3710088/1‑150 (MQ=255)
    ttCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCGGGCGGCGCATGc                                                                                                                                                 <  2:3820630/150‑1 (MQ=255)
     tCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCa                                                                                                                                                >  2:4198384/1‑150 (MQ=255)
            tGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCt                                                                                                                                         <  2:2513340/150‑1 (MQ=255)
             gACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTg                                                                                                                                        >  2:8264744/1‑150 (MQ=255)
              aCTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGc                                                                                                                                       >  2:8570751/1‑150 (MQ=255)
               cTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGcc                                                                                                                                      >  2:3434312/1‑150 (MQ=255)
                tCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCg                                                                                                                                     >  2:4384965/1‑150 (MQ=255)
                 cACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGa                                                                                                                                    <  1:6030947/150‑1 (MQ=255)
                 cACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGa                                                                                                                                    <  1:3161213/150‑1 (MQ=255)
                  aCCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                   >  2:4833593/1‑150 (MQ=255)
                  aCCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAg                                                                                                                                   >  2:2476576/1‑150 (MQ=255)
                   ccTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGc                                                                                                                                  >  2:318353/1‑150 (MQ=255)
                      gCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCTAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGc                                                                                                                               <  2:1439746/150‑1 (MQ=255)
                      gCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGc                                                                                                                               >  1:311364/1‑150 (MQ=255)
                       cccAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCg                                                                                                                              >  2:3387151/1‑150 (MQ=255)
                          aaGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGc                                                                                                                           <  1:1342190/150‑1 (MQ=255)
                               ccGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGaa                                                                                                                       <  1:8309804/149‑1 (MQ=255)
                                  tATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCa                                                                                                                   >  1:6457823/1‑150 (MQ=255)
                                   aTTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAc                                                                                                                  >  1:8748902/1‑150 (MQ=255)
                                     tGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCAccc                                                                                                                <  1:4244697/150‑1 (MQ=255)
                                        tCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCg                                                                                                                         <  2:2743260/138‑1 (MQ=255)
                                        tCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCa                                                                                                             >  1:3538484/1‑150 (MQ=255)
                                        tCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCa                                                                                                             >  2:5335123/1‑150 (MQ=255)
                                            gTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCaaaaa                                                                                                         >  2:6561034/1‑150 (MQ=255)
                                              ccATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACa                                                                                                       >  2:726696/1‑150 (MQ=255)
                                              ccATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACa                                                                                                       >  2:1087964/1‑150 (MQ=255)
                                                aTGACCCCCCGCGAGATCGTCCACGAACTCCACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATc                                                                                                     >  1:6778663/1‑150 (MQ=255)
                                                aTGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATc                                                                                                     <  1:8105448/150‑1 (MQ=255)
                                                aTGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATc                                                                                                     >  1:769967/1‑150 (MQ=255)
                                                 tGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATcc                                                                                                    <  1:8144838/150‑1 (MQ=255)
                                                 tGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATcc                                                                                                    <  1:4630802/150‑1 (MQ=255)
                                                  gACCCCCCGCTAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTCCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGAGCCGAAGTCACCCCCAAAAACATCCt                                                                                                   <  2:7932672/150‑1 (MQ=255)
                                                  gACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCt                                                                                                   <  1:570093/150‑1 (MQ=255)
                                                     cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                >  1:3852614/1‑150 (MQ=255)
                                                     cccccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCtgat                                                                                                >  1:2716904/1‑150 (MQ=255)
                                                        ccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTcc                                                                                                                                           >  2:6844792/1‑104 (MQ=255)
                                                        ccGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTgatgat                                                                                             >  1:696289/1‑150 (MQ=255)
                                                         cgtgaaaTCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATc                                                                                            >  1:754789/7‑150 (MQ=255)
                                                         cgcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATc                                                                                            <  1:725366/150‑1 (MQ=255)
                                                         cgcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATc                                                                                            <  2:854792/150‑1 (MQ=255)
                                                          gcgAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCg                                                                                           >  1:1583075/1‑150 (MQ=255)
                                                            gagaTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGc                                                                                         <  1:5084559/150‑1 (MQ=255)
                                                            gaaaTCGTCTCCGAGCTTGACCGCTTCATCATCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCgga                                                                                         <  1:654722/147‑2 (MQ=255)
                                                             agaTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGcc                                                                                        >  1:2089757/1‑150 (MQ=255)
                                                                 cGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTAc                                                                                    <  2:5178757/150‑1 (MQ=255)
                                                                  gTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTAcc                                                                                   >  2:4709637/1‑150 (MQ=255)
                                                                     cACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGc                                                                                >  2:8426533/1‑150 (MQ=255)
                                                                        gatcTCAACCGCCACATCATCGGGCAGGTCGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTc                                                                             <  2:4920022/147‑1 (MQ=255)
                                                                        gaccTCACCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTc                                                                             <  2:7592754/147‑1 (MQ=255)
                                                                            tCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                         >  1:4528957/1‑150 (MQ=255)
                                                                            tCAACCGCCACATCATCGGCCAGGACGACGCCAAGAGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCa                                                                         >  2:2116110/1‑150 (MQ=255)
                                                                              aaCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAg                                                                       >  2:5902117/1‑150 (MQ=255)
                                                                              aaCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAg                                                                       >  2:6092151/1‑150 (MQ=255)
                                                                               aCCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGa                                                                      >  1:2272033/1‑150 (MQ=255)
                                                                                    ttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGaa                                                                 >  1:7420090/3‑150 (MQ=255)
                                                                                    ttcatcatCGGCCAGAACGACGCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGaa                                                                 >  2:6971940/3‑150 (MQ=255)
                                                                                    cacaTCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCg                                                                            <  1:5277116/139‑1 (MQ=255)
                                                                                       atcatcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATc                                                              <  2:3601183/150‑1 (MQ=255)
                                                                                         catcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgc                                                            <  1:5604379/150‑1 (MQ=255)
                                                                                          atcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgcc                                                           >  1:1067708/1‑150 (MQ=255)
                                                                                           tcGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATcgccc                                                          >  2:1705838/1‑150 (MQ=255)
                                                                                                cAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCgg                                                     <  1:4600188/150‑1 (MQ=255)
                                                                                                 aGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGc                                                    >  2:4325772/1‑150 (MQ=255)
                                                                                                   gacgacGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAACTCGCCCGCCGGCTg                                                  >  2:5624167/1‑150 (MQ=255)
                                                                                                    acgacgCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTgg                                                 <  1:8986869/150‑1 (MQ=255)
                                                                                                    acgacgCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTgg                                                 <  1:4729933/150‑1 (MQ=255)
                                                                                                      gacgCCAAGCGCGCCGTGGCTATTGCTCTGCGCAATCGCTGGCGTCGTCAGCAGCTGGAAGGCCAGATGCGCGAAGAGGTCATGCCAAAAAACATTCTGATGATCGGACCGACCGGCGTTGGCAAGACGGAAATCTCCCGTCGTCTGGcc                                               >  2:283751/1‑150 (MQ=255)
                                                                                                            aaGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTg                                         <  1:101573/150‑1 (MQ=255)
                                                                                                             cgcgcgcCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTgg                                        <  1:996703/149‑1 (MQ=255)
                                                                                                              gcgcgcCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGc                                       <  2:1280230/150‑1 (MQ=255)
                                                                                                              gcgcgcCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGc                                       <  1:1819661/150‑1 (MQ=255)
                                                                                                                gcgcCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGACACCCCCAAAAACATCCAGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCa                                     >  1:2168260/1‑150 (MQ=255)
                                                                                                                    cGTGGCCATTGCCCTGCGCATCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAAcgc                                 <  2:5425524/150‑1 (MQ=255)
                                                                                                                     gTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAAcgcg                                <  1:1705838/150‑1 (MQ=255)
                                                                                                                      tGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACgcgc                               <  1:2495624/150‑1 (MQ=255)
                                                                                                                       ggCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGcc                              >  1:851445/1‑150 (MQ=255)
                                                                                                                        gCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCCACGCGCCg                             >  2:6675145/1‑150 (MQ=255)
                                                                                                                        gCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCg                             <  1:3387152/150‑1 (MQ=255)
                                                                                                                          cATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGtt                           >  2:51413/1‑150 (MQ=255)
                                                                                                                            ttGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTcc                         <  1:2361662/150‑1 (MQ=255)
                                                                                                                             tGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCt                        <  2:8718734/150‑1 (MQ=255)
                                                                                                                               cccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCgg                                                                           >  1:5371769/1‑97 (MQ=255)
                                                                                                                                ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                     >  2:5833162/1‑150 (MQ=255)
                                                                                                                                ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                     <  2:8942719/150‑1 (MQ=255)
                                                                                                                                ccTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCaa                     <  1:589691/150‑1 (MQ=255)
                                                                                                                                 cTGCGCAACCGATGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAg                    <  1:4198077/150‑1 (MQ=255)
                                                                                                                                    cgcAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAACAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCGAAACTGGCCTACGCGCCGTTCCTCAAGGTc                 <  2:7751023/150‑1 (MQ=255)
                                                                                                                                    cgcAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTc                 >  1:9098170/1‑150 (MQ=255)
                                                                                                                                     gcAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCg                <  1:8186724/150‑1 (MQ=255)
                                                                                                                                      cAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGa               <  1:8570894/150‑1 (MQ=255)
                                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAACCTGGCCAACGCGCCGTTCCTCAAGGTCgcg              >  2:6146261/1‑148 (MQ=255)
                                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg              <  2:809655/150‑1 (MQ=255)
                                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg              <  2:1540101/150‑1 (MQ=255)
                                                                                                                                       aaCCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAg              <  1:4946667/150‑1 (MQ=255)
                                                                                                                                            cTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGccac         <  1:7961824/150‑1 (MQ=255)
                                                                                                                                               gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa      >  1:5791413/1‑150 (MQ=255)
                                                                                                                                               gcggcgCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCaa      <  2:8988190/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGCTGGCGGCGCATGCAGCTCCCTGCCGAGCTGCGTGCCGAAGTCACCCCCAAAAACATCCTGATGATCGGCCCTACCGGCGTCGGCAAGACCGAAATCGCCCGCCGGCTGGCCAAACTGGCCAACGCGCCGTTCCTCAAGGTCGAGGCCACCAAGTTC  >  NC_002947/5699303‑5699599

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