Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,397,657 A→T 8.4% L223F (TTA→TTT hsdS → type I restriction modification system specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,397,6570AT 8.4% 226.7 / 10.5 84L223F (TTA→TTThsdStype I restriction modification system specificity protein
Reads supporting (aligned to +/- strand):  ref base A (42/35);  new base T (4/3);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

GGCCGGCGCCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTG  >  NC_002947/5397513‑5397804
                                                                                                                                                |                                                                                                                                                   
ggCCGGCGCCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAAt                                                                                                                                                  <  4:160542/148‑1 (MQ=255)
ggCCGGCGCCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATg                                                                                                                                                 >  5:111389/1‑149 (MQ=255)
 gCCGGCGCCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGt                                                                                                                                                >  3:241809/1‑149 (MQ=255)
    ggCGCCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTgg                                                                                                                                             >  4:211021/1‑149 (MQ=255)
     gcgcCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGc                                                                                                                                            >  1:248264/1‑149 (MQ=255)
       gcCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCt                                                                                                                                          <  1:137290/149‑1 (MQ=255)
        ccAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTTAATGTTGGctgg                                                                                                                                         <  6:179848/149‑4 (MQ=255)
               tGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGAtt                                                                                                                                      >  5:304562/1‑145 (MQ=255)
               tGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGAtt                                                                                                                                      <  6:304562/145‑1 (MQ=255)
                 ccATAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGTCAACATTAAATGTTGGCCTGATTCGTAg                                                                                                                                 <  1:201614/148‑1 (MQ=255)
                 ccAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGt                                                                                                                                >  3:359153/1‑149 (MQ=255)
                 ccAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGt                                                                                                                                <  6:111389/149‑1 (MQ=255)
                               cGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGATAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTg                                                                                                                                              >  5:214334/1‑121 (MQ=255)
                               cGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTg                                                                                                                                              <  6:214334/121‑1 (MQ=255)
                               cGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACctct                                                                                                                  <  4:241809/149‑1 (MQ=255)
                                   tAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGc                                                                                                                <  4:32729/147‑1 (MQ=255)
                                   tAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGc                                                                                                                >  3:32729/1‑147 (MQ=255)
                                     gTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTc                                                                                                             <  2:248264/148‑1 (MQ=255)
                                      tCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGcc                                                                                                                                           >  4:86005/1‑117 (MQ=255)
                                      tCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGcc                                                                                                                                           <  3:86005/117‑1 (MQ=255)
                                       ccccAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGc                                                                                                                        >  1:216204/1‑135 (MQ=255)
                                       ccccAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGc                                                                                                                        <  2:216204/135‑1 (MQ=255)
                                        cccAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTTAATGTTGGctgg                                                                                                                                         >  4:253275/1‑114 (MQ=255)
                                        cccAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTTAATGTTGGctgg                                                                                                                                         <  3:253275/117‑4 (MQ=255)
                                               tgCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGccca                                                                                                   >  1:128004/1‑146 (MQ=255)
                                                gCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGaa                                                                                                  >  7:159568/1‑148 (MQ=255)
                                                 cATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGaac                                                                                                 >  5:111248/1‑148 (MQ=255)
                                                    gTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAacaaca                                                                                             >  2:47127/1‑149 (MQ=255)
                                                      ctcGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGGGAAAGCAGTTCTTGGGAGATAATCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCACCTAGCCGAACAACAcc                                                                                           >  5:266124/1‑149 (MQ=255)
                                                       tcGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCg                                                                                          >  5:305227/1‑149 (MQ=255)
                                                         gAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCa                                                                                        >  3:86062/1‑149 (MQ=255)
                                                         gAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCACCGAACAACACCGCa                                                                                        >  5:308021/1‑149 (MQ=255)
                                                             gACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGt                                                                                    >  6:372353/1‑149 (MQ=255)
                                                              aCTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGtt                                                                                   <  3:211021/149‑1 (MQ=255)
                                                               tctcATAATTTGGTTGCTACAAATTGACCTTATGAGAAAGCAGTTATTGGTAGATACTCGAAACCTCGCCCAGACAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTTCCTCCCCTCGCCGAACAACACCGCATCGtt                                                                                   <  4:249062/145‑1 (MQ=255)
                                                               cTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAAGATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTg                                                                                  <  6:81218/149‑1 (MQ=255)
                                                               cTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTg                                                                                  <  4:221785/149‑1 (MQ=255)
                                                                  cATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTcccc                                                                                                          <  4:130045/122‑1 (MQ=255)
                                                                  cATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTcccc                                                                                                          >  3:130045/1‑122 (MQ=255)
                                                                   aTAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACATCTGCCTCCCCTCGCCGAAAAACACCGCATCGtt                                                                                   >  3:249062/1‑144 (MQ=255)
                                                                            ttGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCgaag                                                                      >  7:333511/1‑148 (MQ=255)
                                                                              gCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCgaagaa                                                                    >  8:221354/1‑148 (MQ=255)
                                                                                tACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAg                                                                                                                                 <  5:320604/85‑1 (MQ=255)
                                                                                tACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAg                                                                                                                                 >  6:320604/1‑85 (MQ=255)
                                                                                tACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCcgc                                                                                                                           >  7:337957/1‑91 (MQ=255)
                                                                                tACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCcgc                                                                                                                           <  8:337957/91‑1 (MQ=255)
                                                                                tACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTc                                                                                                        <  2:245054/110‑1 (MQ=255)
                                                                                tACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTc                                                                                                        >  1:245054/1‑110 (MQ=255)
                                                                                  cAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAAt                                                               >  4:279624/1‑149 (MQ=255)
                                                                                    aaGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATgg                                                             <  6:305227/149‑1 (MQ=255)
                                                                                     aGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGc                                                            >  4:51248/1‑149 (MQ=255)
                                                                                      gTGACCTTATGCGAAAGCAGTTCTTTGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGcc                                                           <  8:159568/149‑1 (MQ=255)
                                                                                      gTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCACCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGcc                                                           <  6:308021/149‑1 (MQ=255)
                                                                                           cTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATg                                                      <  2:128004/149‑1 (MQ=255)
                                                                                           cTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATg                                                      >  3:173849/1‑149 (MQ=255)
                                                                                           cTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATg                                                      <  3:279624/149‑1 (MQ=255)
                                                                                           cTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATg                                                      >  8:283965/1‑149 (MQ=255)
                                                                                            ttATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGc                                                                                                                                            >  4:197437/1‑62 (MQ=255)
                                                                                                gCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGAcc                                                 >  3:216621/1‑149 (MQ=255)
                                                                                                 cGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTc                                                                                                             <  1:36723/88‑1 (MQ=255)
                                                                                                 cGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTc                                                                                                             >  2:36723/1‑88 (MQ=255)
                                                                                                      gCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTgg                                           <  6:111248/149‑1 (MQ=255)
                                                                                                        aGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGaa                                         <  7:221354/149‑1 (MQ=255)
                                                                                                          ttCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGc                                       >  8:345029/1‑149 (MQ=255)
                                                                                                                  agaTACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAgg                               <  7:345029/149‑1 (MQ=255)
                                                                                                                   gaTACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAgg                               <  4:173849/148‑1 (MQ=255)
                                                                                                                        tCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTc                                                                          >  7:77814/1‑100 (MQ=255)
                                                                                                                        tCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTc                                                                          <  8:77814/100‑1 (MQ=255)
                                                                                                                          gAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGcc                       >  1:250029/1‑149 (MQ=255)
                                                                                                                          gAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGcc                       <  7:283965/149‑1 (MQ=255)
                                                                                                                            aCCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGa                     <  4:359153/149‑1 (MQ=255)
                                                                                                                             cccTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGcc                       >  3:195687/1‑146 (MQ=255)
                                                                                                                             cccTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGcc                       <  4:195687/146‑1 (MQ=255)
                                                                                                                              ccTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGt                   >  4:18690/1‑149 (MQ=255)
                                                                                                                                  gCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAAAAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATCCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGc               >  3:210898/1‑149 (MQ=255)
                                                                                                                                  atcaggCCAACATTTAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGc               >  2:19403/6‑149 (MQ=255)
                                                                                                                                   cccAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGcc              >  6:88759/1‑149 (MQ=255)
                                                                                                                                    ccAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGAc                          <  7:269741/136‑1 (MQ=255)
                                                                                                                                    ccAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGAc                          >  8:269741/1‑136 (MQ=255)
                                                                                                                                     aggCCAACATTTAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGaa                                         <  3:276007/118‑1 (MQ=255)
                                                                                                                                     aggCCAACATTTAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGaa                                         >  4:276007/3‑120 (MQ=255)
                                                                                                                                        ccAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCACTCGCCGAACAACACCGCATCGTTGCCAAAGTAGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCTCCCACGt         <  6:266124/149‑1 (MQ=255)
                                                                                                                                           aCATTTAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACAAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTcc       >  8:320499/1‑148 (MQ=255)
                                                                                                                                           aCATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTcc       >  8:69234/1‑148 (MQ=255)
                                                                                                                                                aaaTGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTg  <  2:250029/148‑1 (MQ=255)
                                                                                                                                                aaaTGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTg  <  8:333511/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCCGGCGCCAACATTGCCAGAGCTATTTGCCGTGTAGTCCCCAGTGTGCATGTCTCGAAGGACTACATAATTTGGTTGCTACAAAGTGACCTTATGCGAAAGCAGTTCTTGGGAGATACTCGAACCCTCGCCCAGCCAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTG  >  NC_002947/5397513‑5397804

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: