Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 898,295 T→G 5.1% F60C (TTC→TGC)  PP_0781 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947898,2950TG 5.1% 291.9 / 3.8 92F60C (TTC→TGC) PP_0781hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base G (3/2);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGC  >  NC_002947/898149‑898443
                                                                                                                                                  |                                                                                                                                                    
ccTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTc                                                                                                                                                     <  8:153695/148‑1 (MQ=255)
 cTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTa                                                                                                                                                   >  2:258012/1‑149 (MQ=255)
    gTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTAc                                                                                                                                                  <  7:91256/147‑1 (MQ=255)
    gTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCAAGGTCCGCACCCGCCTCACTGAAGCATTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCa                                                                                                                                                <  7:328453/149‑1 (MQ=14)
         ctcctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGCCGCTGTGCTCTCAGCCTGCTACCAGGCt                                                                                                                                            >  6:89587/1‑148 (MQ=255)
            ctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAAt                                                                                                                                        <  6:134829/149‑1 (MQ=255)
            ctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAAt                                                                                                                                        <  5:368712/149‑1 (MQ=255)
            ctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAAt                                                                                                                                        <  5:165665/149‑1 (MQ=255)
            ctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTGGTAGAAGGCTGa                                                                                                                                          <  6:46617/147‑1 (MQ=255)
            ctGGCACGAGATCGAGTCACGTGTCATCGACATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGaa                                                                                                                                         <  4:275420/148‑1 (MQ=255)
            agggCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAAt                                                                                                                                        <  8:111574/147‑1 (MQ=255)
               gCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCTATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGc                                                                                                                                     <  6:310806/149‑1 (MQ=255)
                   gagaTCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGAc                                                                                                                                 >  4:212287/1‑149 (MQ=255)
                     gaTCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACgg                                                                                                                               >  7:11582/1‑149 (MQ=255)
                     gaTCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACgg                                                                                                                               >  3:253677/1‑149 (MQ=255)
                          aGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGcac                                                                                                                          <  1:125897/149‑1 (MQ=255)
                              aCGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCtt                                                                                                                                                      <  4:35412/117‑1 (MQ=255)
                              aCGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCtt                                                                                                                                                      >  3:35412/1‑117 (MQ=255)
                                   tCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGc                                                                                                                 <  2:241847/149‑1 (MQ=255)
                                    cATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGc                                                                                                                 >  1:212245/1‑148 (MQ=255)
                                     aTCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGc                                                                                                                 >  2:170630/1‑147 (MQ=255)
                                          cATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCt                                                                                                          >  7:275118/1‑149 (MQ=255)
                                            tggtggGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCGGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTtc                                                                                                        <  6:100185/149‑1 (MQ=255)
                                              gtggGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTtcgt                                                                                                      >  3:262247/1‑149 (MQ=255)
                                              gtggGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTtcgt                                                                                                      <  4:127929/149‑1 (MQ=255)
                                                gggCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTtcgtc                                                                                                     >  2:252186/1‑148 (MQ=255)
                                                 ggCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGa                                                                                                   >  8:175904/1‑149 (MQ=255)
                                                 ggCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGa                                                                                                   <  1:258012/149‑1 (MQ=255)
                                                   cTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTGGTAGAAGGCTGagc                                                                                                                                        >  4:286360/1‑108 (MQ=17)
                                                    tACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCt                                                                                                >  6:204455/1‑149 (MQ=255)
                                                       ccAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGa                                                                                             >  1:371848/1‑149 (MQ=255)
                                                        cAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGAc                                                                                            <  3:133557/149‑1 (MQ=255)
                                                          aGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCCGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACAt                                                                                          >  5:1292/1‑149 (MQ=255)
                                                           gCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATc                                                                                         <  4:253677/149‑1 (MQ=255)
                                                               tGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACa                                                                                     <  6:346867/149‑1 (MQ=255)
                                                                 gCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACAAGGCTGAATGGCAGACGGCGAACATGCTACGCCAGCGCTTCGACGACCTCGACATCAACAg                                                                                    >  4:312171/1‑148 (MQ=255)
                                                                  ccAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCa                                                                                  >  1:162924/1‑149 (MQ=255)
                                                                   cAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCAt                                                                                 >  4:297594/1‑149 (MQ=255)
                                                                    aGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCAtc                                                                                <  2:303324/149‑1 (MQ=255)
                                                                         cGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCAtcgtcg                                                                            <  2:212245/148‑1 (MQ=255)
                                                                         cGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTcga                                                                           <  4:262247/149‑1 (MQ=255)
                                                                          gCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTcgac                                                                          <  6:53555/149‑1 (MQ=255)
                                                                          gCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTcgac                                                                          <  2:162924/149‑1 (MQ=255)
                                                                            aCCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTcgacg                                                                         >  6:12394/1‑148 (MQ=255)
                                                                                gCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGCCCTCGACATCAACAGCATCGTCGACGAActg                                                                    >  2:66523/1‑149 (MQ=255)
                                                                                 ccTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAAc                                                                      >  7:68615/1‑146 (MQ=255)
                                                                                 ccTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAAc                                                                      <  8:68615/146‑1 (MQ=255)
                                                                                   tCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAActgct                                                                  >  7:163107/1‑148 (MQ=255)
                                                                                    cAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGa                                                                <  7:175904/149‑1 (MQ=255)
                                                                                     aaTGAAGAAGTTGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAGCTGCTGAg                                                               >  2:206399/1‑149 (MQ=255)
                                                                                     aaTGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGa                                                                <  5:25017/148‑1 (MQ=255)
                                                                                     aaTGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAg                                                               >  8:108682/1‑149 (MQ=255)
                                                                                     aaTGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAg                                                               >  4:118643/1‑149 (MQ=255)
                                                                                      aTGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAgg                                                              >  7:143273/1‑149 (MQ=255)
                                                                                       tGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAggg                                                             <  5:204455/149‑1 (MQ=255)
                                                                                       tGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAggg                                                             >  8:24028/1‑149 (MQ=255)
                                                                                         aagaagTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGt                                                            <  1:252186/148‑1 (MQ=255)
                                                                                         aagaagTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTc                                                           >  5:2182/1‑149 (MQ=255)
                                                                                           gaagTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGt                                                         >  4:25715/1‑149 (MQ=255)
                                                                                            aagTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTg                                                        <  3:212287/149‑1 (MQ=255)
                                                                                               tCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTgcag                                                     <  8:108216/149‑1 (MQ=255)
                                                                                               tCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTgcag                                                     >  7:131015/1‑149 (MQ=255)
                                                                                                 gCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTgcag                                                     >  4:347087/1‑147 (MQ=255)
                                                                                                   cTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGa                                                 <  3:271864/149‑1 (MQ=255)
                                                                                                      gCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATgg                                              <  6:2182/149‑1 (MQ=255)
                                                                                                        tGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGc                                             >  1:160255/1‑148 (MQ=255)
                                                                                                          ggCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCga                                           <  2:103777/148‑1 (MQ=255)
                                                                                                            ccTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCgatga                                        >  7:325927/1‑149 (MQ=255)
                                                                                                             cTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCgatgat                                       >  5:261000/1‑149 (MQ=255)
                                                                                                                aaCCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgt                                    >  5:104314/1‑149 (MQ=255)
                                                                                                                    tGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgtcgtc                                >  7:96875/1‑149 (MQ=255)
                                                                                                                     gcgcCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgtcgtcg                               >  4:37323/1‑149 (MQ=255)
                                                                                                                      cgcCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgtcgtcg                               <  7:24028/148‑1 (MQ=255)
                                                                                                                      cgcCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgtcgtc                                <  4:183822/147‑1 (MQ=255)
                                                                                                                       gcCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCgatgat                                       >  7:115677/1‑139 (MQ=255)
                                                                                                                       gcCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCgatgat                                       <  8:115677/139‑1 (MQ=255)
                                                                                                                       gcCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCgg                              >  1:60407/1‑148 (MQ=255)
                                                                                                                        ccGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGgc                                                                                                                              >  8:332008/1‑51 (MQ=255)
                                                                                                                        ccGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGgc                                                                                                                              <  7:332008/51‑1 (MQ=255)
                                                                                                                                aCGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCg                                            <  6:384202/125‑1 (MQ=255)
                                                                                                                                aCGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCg                                            >  5:384202/1‑125 (MQ=255)
                                                                                                                                aCGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGt                     >  7:175399/1‑148 (MQ=255)
                                                                                                                                   cTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTAcc                 <  6:261000/149‑1 (MQ=255)
                                                                                                                                    tgtgCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgtc                                   <  5:3066/130‑1 (MQ=255)
                                                                                                                                    tgtgCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGAtcgtc                                   >  6:3066/1‑130 (MQ=255)
                                                                                                                                      tgCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGt              <  2:160255/149‑1 (MQ=255)
                                                                                                                                      tgCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGt              >  2:141813/1‑149 (MQ=255)
                                                                                                                                      tgCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGt              <  8:175399/149‑1 (MQ=255)
                                                                                                                                        cgcTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTgg            >  8:332785/1‑149 (MQ=255)
                                                                                                                                          ttCAGCCTGGTAGAAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCg                <  7:322879/142‑1 (MQ=255)
                                                                                                                                          ttCAGCCTGGTAGAAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCg                >  8:322879/2‑143 (MQ=255)
                                                                                                                                             aGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGc       <  1:170630/149‑1 (MQ=255)
                                                                                                                                               ccTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCtt     <  8:275118/149‑1 (MQ=255)
                                                                                                                                                cTTCTACCAGGCTGAATTGAAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCTTCGACGAACTGCTGAGGGTCGTGCAGCATATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTTGCCTGCTTg    <  6:353958/149‑1 (MQ=14)
                                                                                                                                                 ttCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTgg   >  5:111482/1‑149 (MQ=255)
                                                                                                                                                  tCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGgc  >  7:146003/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGACATCAACAGCATCGTCGACGAACTGCTGAGGGTCGTGCAGCAGATGGCGATGATCGTCGTCGGCAGCACCGTTACCGGTGGCCTGCTTGGC  >  NC_002947/898149‑898443

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