Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,426,086 G→C 7.4% G208R (GGC→CGC)  PP_4765 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,426,0860GC 7.4% 244.9 / 9.2 81G208R (GGC→CGC) PP_4765hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base C (3/3);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

AGGTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAC  >  NC_002947/5425940‑5426233
                                                                                                                                                  |                                                                                                                                                   
aGGTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCggg                                                                                                                                                    <  1:155124/148‑1 (MQ=255)
   tGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTg                                                                                                                                                 >  5:297079/1‑148 (MQ=255)
    gAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTg                                                                                                                                                 >  5:108702/1‑147 (MQ=255)
     aaaGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGa                                                                                                                                               <  7:285765/148‑1 (MQ=255)
           tCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  2:281838/1‑148 (MQ=255)
           tCGGTCGTGGTGCTGGCGCTGGGTGGCGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  5:63608/1‑148 (MQ=255)
             ggTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTa                                                                                                                                       <  4:80093/148‑1 (MQ=255)
              gTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  4:56460/149‑1 (MQ=255)
              gTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCGGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  5:127033/149‑1 (MQ=255)
               tCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTc                                                                                                                                    <  2:108069/149‑1 (MQ=255)
                  tggtgCGGGCGCTGGGTGGTGCCAGCGGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATAGCTCGCCTTTGCAGCCGAGCACCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAg                                                                                                                                 <  1:198037/149‑1 (MQ=255)
                     tgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGAc                                                                                                                              <  5:254951/149‑1 (MQ=255)
                        tGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGAc                                                                                                                              >  3:8417/1‑146 (MQ=255)
                        tGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGAc                                                                                                                              >  3:12777/1‑146 (MQ=255)
                          gggcTGGGTGGTGGCAGCTGGGGGCGGCTGGGTTCTGACGGCAGCTGGGATGCGGTGCTGGGTGATCGGGCTGTGGTTATCTCGCCTTTGCTGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGActa                                                                                                                            <  4:8417/144‑3 (MQ=255)
                          gcgcTGGGTTGTGCCAGCTTGGCGCGGATGGGTTCTGACGGAAGCTGGCTACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCTGTTTTTAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACaa                                                                                                                            <  4:12777/146‑1 (MQ=255)
                            gcTGGGTGGTGCCAGCTGGGGGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcg                                                                                                                       <  7:271485/149‑1 (MQ=255)
                            gcTGGGTGGTGCCAGCTGGGGGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcg                                                                                                                       <  7:275796/149‑1 (MQ=255)
                             cTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     >  3:214045/1‑134 (MQ=255)
                             cTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  4:214045/134‑1 (MQ=255)
                                       ccAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCg                                                                                                            <  5:263862/149‑1 (MQ=255)
                                           cTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTg                                                                                                         >  3:202108/1‑148 (MQ=255)
                                                cgcgcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGCAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTgaagaa                                                                                                    >  1:303749/1‑148 (MQ=255)
                                                   gcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACtt                                                                                                                                    >  4:79537/1‑112 (MQ=255)
                                                   gcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACtt                                                                                                                                    <  3:79537/113‑2 (MQ=255)
                                                     cTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTAc                                                                                                                                      <  5:182346/109‑1 (MQ=255)
                                                     cTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTAc                                                                                                                                      >  6:182346/1‑109 (MQ=255)
                                                      tGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATTCCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         <  1:84759/105‑1 (MQ=255)
                                                      tGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  2:84759/1‑105 (MQ=255)
                                                      tGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGACTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCC                                                                                             >  6:198132‑M1/1‑149 (MQ=255)
                                                       gggTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCC                                                                                            >  4:116158‑M1/1‑149 (MQ=255)
                                                        ggTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCC                                                                                             <  4:202108‑M1/147‑1 (MQ=255)
                                                         gTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCT                                                                                           >  5:30684‑M1/1‑148 (MQ=255)
                                                             tGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAg                                                                                                                           <  1:278166/112‑1 (MQ=255)
                                                             tGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAg                                                                                                                           >  2:278166/1‑112 (MQ=255)
                                                              gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcg                                                                                                                                          <  1:291866/96‑1 (MQ=255)
                                                              gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcg                                                                                                                                          >  2:291866/1‑96 (MQ=255)
                                                              gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGt                                                                                     >  2:97726/1‑149 (MQ=255)
                                                                 ggCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTccc                                                                                  >  7:93192/1‑149 (MQ=255)
                                                                  gcctcTGGCAACCGCTGCTGTCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCg                                                                                 <  8:148515/145‑1 (MQ=255)
                                                                       tGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGc                                                                                       >  3:224977/1‑138 (MQ=255)
                                                                       tGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGc                                                                                       <  4:224977/138‑1 (MQ=255)
                                                                         gCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAgcg                                                                          >  5:320688/1‑149 (MQ=255)
                                                                                         gAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAagggcga                                                          >  1:175897‑M2/1‑142 (MQ=255)
                                                                                         gAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGAGTACCGGGCAGCTAAACGCGCAagg                                                              <  7:67382‑M2/145‑4 (MQ=255)
                                                                                            cGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGtt                                                       <  1:281838/149‑1 (MQ=255)
                                                                                               gCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGCGCGCACCGCTGAAGAACATCGCCCTtgcgcggtgcgctgcccgggacgctgagggcgatgttct                                                    <  6:108702‑M1/149‑40 (MQ=255)
                                                                                               gCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGCAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAt                                                    <  2:303749/149‑1 (MQ=255)
                                                                                                cTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   <  8:93192/149‑1 (MQ=255)
                                                                                                    ggATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCa                                               >  6:331302/1‑149 (MQ=255)
                                                                                                       tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAag                                                               <  7:144250‑M2/130‑3 (MQ=255)
                                                                                                       tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAag                                                               >  8:144250‑M2/1‑128 (MQ=255)
                                                                                                       tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCg                                            >  3:257907/1‑149 (MQ=255)
                                                                                                         tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   <  2:281624/140‑1 (MQ=255)
                                                                                                         tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   >  1:281624/1‑140 (MQ=255)
                                                                                                          ctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCa                                                                                         <  5:60465/101‑1 (MQ=255)
                                                                                                          ctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCa                                                                                         >  6:60465/1‑101 (MQ=255)
                                                                                                           tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTtgcgcggtgcgctgcccgggacgctgagggcgatgttc                                                     <  7:163076‑M1/136‑39 (MQ=255)
                                                                                                           tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTtgcgcggtgcgctgcccgggacgctgagggcgatgttc                                                     >  8:163076‑M1/1‑98 (MQ=255)
                                                                                                           tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAagg                                                              >  3:32754‑M2/1‑124 (MQ=255)
                                                                                                           tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAagg                                                              <  4:32754‑M2/127‑4 (MQ=255)
                                                                                                             gCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAg                                      <  3:116158/149‑1 (MQ=255)
                                                                                                              ccTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAgg                                     >  2:211360/1‑149 (MQ=255)
                                                                                                              ccTTTGCAGCCGAGCAACTGCGGGTTTGAAGTTGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCTCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGATCATTCTCACCGCTCAAgg                                     <  7:269119/149‑1 (MQ=255)
                                                                                                                  tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGTCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAggatt                                  <  6:320688/148‑4 (MQ=255)
                                                                                                                  tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGt                                  <  5:331302/148‑1 (MQ=255)
                                                                                                                   gCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTg                                 <  2:220337/148‑1 (MQ=255)
                                                                                                                     aGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTgg                                >  5:103131/1‑147 (MQ=255)
                                                                                                                     aGCCGAGCAACTGCGGTTTTGAAGTAGCGGACTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGaa                              <  5:198132/149‑1 (MQ=255)
                                                                                                                       ccGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                            <  6:103131/149‑1 (MQ=255)
                                                                                                                       ccGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCAACGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                            >  8:130748/1‑149 (MQ=255)
                                                                                                                        cGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                            >  1:242804/1‑148 (MQ=255)
                                                                                                                             aaCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAACTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCt                      >  7:212005/1‑149 (MQ=255)
                                                                                                                              aCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTg                     <  6:297079/149‑1 (MQ=255)
                                                                                                                              aCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTg                     >  1:45250/1‑149 (MQ=255)
                                                                                                                                tGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                   >  1:271293/1‑149 (MQ=255)
                                                                                                                                  cGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                   <  1:211360/147‑1 (MQ=255)
                                                                                                                                    gTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTAt               <  6:30684/149‑1 (MQ=255)
                                                                                                                                     ttttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATg              >  3:251589/1‑149 (MQ=255)
                                                                                                                                        ttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCt           >  7:224235/4‑149 (MQ=255)
                                                                                                                                          aaGTAGAGTTCTGGAGCGATCTACTCAAGGACAAGTTCGCTGGCGAACCGATGAAGAACATCGCCTTCATCGTCCCGGGCAGCGAACCGCGAAAGGGCTAGTTAATTATCAACGATAAATGCGTGGATGGCAGCCTTGTGTATGCCTg          <  8:67382/148‑1 (MQ=255)
                                                                                                                                                cGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCa   >  1:204673/1‑149 (MQ=255)
                                                                                                                                                 gggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAc  <  2:204673/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGGTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAC  >  NC_002947/5425940‑5426233

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: