Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,943,665 C→A 5.2% R407R (CGC→CGA rbbA → ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,943,6650CA 5.2% 1080.5 / 28.9 317R407R (CGC→CGArbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base C (156/144);  new base A (11/6);  total (167/150)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

ACCCGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCA  >  NC_002947/5943516‑5943813
                                                                                                                                                     |                                                                                                                                                    
 cccGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGcc                                                                                                                                                     >  1:6270386/1‑150 (MQ=255)
  ccGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCt                                                                                                                                                    >  2:315273/1‑150 (MQ=255)
     gCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTc                                                                                                                                                 >  2:3576374/1‑150 (MQ=255)
        gAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCt                                                                                                                                              <  1:3975168/150‑1 (MQ=255)
        gAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCt                                                                                                                                              <  1:1641534/150‑1 (MQ=255)
         aaCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTg                                                                                                                                             <  2:3953685/150‑1 (MQ=255)
         aaCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTg                                                                                                                                             <  2:5019476/150‑1 (MQ=255)
          aCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTgg                                                                                                                                            <  1:442070/150‑1 (MQ=255)
           ccTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGAGCGCCGCGCATCGATGAGCTGATCGAGCGCTTCGCCCGCGCCGCTATCGCCGACCAACCGTCTGGCTCCCTGCCCCTTGGTTTGCTCGAGCGACAGTCGCTGGc                                                                                                                                           <  2:4016038/150‑1 (MQ=255)
            cTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCa                                                                                                                                          <  2:1785159/150‑1 (MQ=255)
                 ccTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGc                                                                                                                                     >  2:5891692/1‑150 (MQ=255)
                  cTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCg                                                                                                                                    <  2:6291832/150‑1 (MQ=255)
                   tGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCgg                                                                                                                                   <  2:4663548/150‑1 (MQ=255)
                   ggCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCgg                                                                                                                                   <  2:3414733/149‑1 (MQ=255)
                    gCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGt                                                                                                                                  <  2:5404559/150‑1 (MQ=255)
                     cATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGtt                                                                                                                                 >  2:1730813/1‑150 (MQ=255)
                      aTGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTc                                                                                                                                >  1:223316/1‑150 (MQ=255)
                      aTGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTc                                                                                                                                <  1:789690/150‑1 (MQ=255)
                      aTGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTc                                                                                                                                <  1:5096081/150‑1 (MQ=255)
                      aTGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTc                                                                                                                                <  2:6129602/150‑1 (MQ=255)
                        gCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTg                                                                                                                              <  1:5928135/150‑1 (MQ=255)
                          ccGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCa                                                                                                                            <  1:2092568/150‑1 (MQ=255)
                              cTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGc                                                                                                                        >  2:316995/1‑150 (MQ=255)
                                    gACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGaa                                                                                                                  <  1:6126420/150‑1 (MQ=255)
                                           ccAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGa                                                                                                           >  1:6428550/1‑150 (MQ=255)
                                            cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGAt                                                                                                          >  2:529929/1‑150 (MQ=255)
                                            cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGAt                                                                                                          >  1:1595589/1‑150 (MQ=255)
                                            cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGAt                                                                                                          >  2:4557766/1‑150 (MQ=255)
                                             aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATc                                                                                                         <  1:3724954/150‑1 (MQ=255)
                                             aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATc                                                                                                         >  1:6719305/1‑150 (MQ=255)
                                             aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATc                                                                                                         >  2:4339719/1‑150 (MQ=255)
                                             aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATc                                                                                                         >  1:5706057/1‑150 (MQ=255)
                                              aaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATcc                                                                                                        <  2:2073190/150‑1 (MQ=255)
                                              aaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATcc                                                                                                        >  1:2010426/1‑150 (MQ=255)
                                                  cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  2:2478324/150‑1 (MQ=255)
                                                  cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  1:5877987/150‑1 (MQ=255)
                                                  cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  2:473762/150‑1 (MQ=255)
                                                  cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  1:2504477/150‑1 (MQ=255)
                                                  cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  1:1739756/150‑1 (MQ=255)
                                                   gagaGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcgac                                                                                                   <  1:5378466/150‑1 (MQ=255)
                                                       gcgcCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGACGTCTTGATCCTCGACGAAc                                                                                               <  1:1569095/150‑1 (MQ=255)
                                                        cgcCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAAcc                                                                                              <  1:4999239/150‑1 (MQ=255)
                                                         gcCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCg                                                                                             >  2:3811998/1‑150 (MQ=255)
                                                            cAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGAcc                                                                                          >  1:5284426/1‑150 (MQ=255)
                                                            cAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGAcc                                                                                          >  1:6171323/1‑150 (MQ=255)
                                                               cgcATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTcc                                                                                       <  2:2561799/150‑1 (MQ=255)
                                                                     gATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTc                                                                                 >  1:6247468/1‑150 (MQ=255)
                                                                            tGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccg                                                                          <  1:1443436/150‑1 (MQ=255)
                                                                             gATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgc                                                                         >  1:5514041/1‑150 (MQ=255)
                                                                               tCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgccg                                                                       >  1:4642632/1‑150 (MQ=255)
                                                                               tCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgccg                                                                       >  2:1367294/1‑150 (MQ=255)
                                                                               tCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgccg                                                                       >  2:1403778/1‑150 (MQ=255)
                                                                                 gAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCgccgcc                                                                     >  2:4983055/1‑150 (MQ=255)
                                                                                 gAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCgccgcc                                                                     >  2:824762/1‑150 (MQ=255)
                                                                                  aGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAAACGACCTCCGGCGTCGACCCCGccgccc                                                                    <  1:5103762/150‑1 (MQ=255)
                                                                                   gcgcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGccgcccg                                                                   >  1:1363396/1‑150 (MQ=255)
                                                                                      cTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcga                                                                <  2:4854448/150‑1 (MQ=255)
                                                                                       ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgac                                                               >  2:1584880/1‑150 (MQ=255)
                                                                                       ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgac                                                               >  2:2204263/1‑150 (MQ=255)
                                                                                       ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgac                                                               <  2:1901084/150‑1 (MQ=255)
                                                                                       ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgac                                                               >  2:6194121/1‑150 (MQ=255)
                                                                                           aCCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACt                                                           >  1:559371/1‑150 (MQ=255)
                                                                                           aCCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACt                                                           >  1:5836712/1‑150 (MQ=255)
                                                                                             cTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTCCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTc                                                         >  2:6582682/1‑150 (MQ=255)
                                                                                              tCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCt                                                        >  1:4724804/1‑150 (MQ=255)
                                                                                              tCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCt                                                        <  1:3245405/150‑1 (MQ=255)
                                                                                              tCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCt                                                        <  1:3876878/150‑1 (MQ=255)
                                                                                                gccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTgg                                                      >  1:3408944/1‑150 (MQ=255)
                                                                                                gccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTgg                                                      >  1:2813066/1‑150 (MQ=255)
                                                                                                gccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTgg                                                      >  2:1605290/1‑150 (MQ=255)
                                                                                                gccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCGTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGCTCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACGTCTgg                                                      >  1:2290897/1‑150 (MQ=255)
                                                                                                 ccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGgc                                                     <  2:5843384/150‑1 (MQ=255)
                                                                                                    cTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGcc                                                  <  1:5303554/150‑1 (MQ=255)
                                                                                                          ccGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTgg                                            <  2:5345263/150‑1 (MQ=255)
                                                                                                            gACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTa                                          >  1:427482/1‑150 (MQ=255)
                                                                                                            gACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTa                                          <  1:3507869/150‑1 (MQ=255)
                                                                                                               cAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGaa                                       <  1:2708923/150‑1 (MQ=255)
                                                                                                               cAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGaa                                       <  1:6864582/150‑1 (MQ=255)
                                                                                                                aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                      >  1:3022346/1‑150 (MQ=255)
                                                                                                                 aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACt                                     >  2:1184526/1‑150 (MQ=255)
                                                                                                                   cGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGt                                   >  1:5984900/1‑150 (MQ=255)
                                                                                                                    gTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTc                                  <  1:1154155/150‑1 (MQ=255)
                                                                                                                     gtcggCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTcc                                 >  2:2704539/4‑150 (MQ=255)
                                                                                                                      cTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTccc                                <  2:721645/150‑1 (MQ=255)
                                                                                                                        ggCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgc                              >  2:4637834/1‑150 (MQ=255)
                                                                                                                        ggCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgc                              <  1:3595503/150‑1 (MQ=255)
                                                                                                                         gcgcCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgcg                             <  1:103631/150‑1 (MQ=255)
                                                                                                                             ccTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCa                         >  2:3644741/1‑150 (MQ=255)
                                                                                                                               tGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAgg                       >  1:2603295/1‑150 (MQ=255)
                                                                                                                               tGCCCCTTGGTGTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAgg                       <  2:4129278/150‑1 (MQ=255)
                                                                                                                                   ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGt                   <  2:27810/150‑1 (MQ=255)
                                                                                                                                    cTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTg                  >  1:4525934/1‑150 (MQ=255)
                                                                                                                                     ttGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                 <  1:4081246/150‑1 (MQ=255)
                                                                                                                                        gccactgcCAGCGACAGGCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCa              >  2:2908098/7‑150 (MQ=255)
                                                                                                                                             cgcCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTc         <  2:4260499/150‑1 (MQ=255)
                                                                                                                                             cgcCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTc         <  2:5747341/150‑1 (MQ=255)
                                                                                                                                                cAGCGACAGGCGCTGGCGGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCctc      >  2:894294/1‑150 (MQ=255)
                                                                                                                                                    gcCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTcca  <  2:6243139/150‑1 (MQ=255)
                                                                                                                                                    gcCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTcca  >  1:24792/1‑150 (MQ=255)
                                                                                                                                                    gcCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTcca  <  2:1262711/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCCGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCA  >  NC_002947/5943516‑5943813

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: