Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,709,374 A→G 7.0% E1131G (GAA→GGA)  PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,709,3740AG 7.0% 248.6 / 6.9 86E1131G (GAA→GGA) PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base G (3/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

GTACGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGAC  >  NC_002947/5709228‑5709521
                                                                                                                                                  |                                                                                                                                                   
gTACGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAAc                                                                                                                                                   <  1:277892/149‑1 (MQ=255)
   cGGTGCAGCCATCTAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGATACTGCCGAGCAGATACAGCAATACCTCGATTATCGCTAACgg                                                                                                                                                 <  8:14567/148‑1 (MQ=255)
         aGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                           >  8:134926/1‑148 (MQ=255)
           ccATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGAt                                                                                                                                        >  3:224778/1‑149 (MQ=255)
           ccATCGAAAACGAAATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCACGCGTTCAAGAGACTTCCGAGAAGATACAGCAATACCTCGATTATCGCGAACTGGTCCGCGAt                                                                                                                                        >  4:231242/1‑149 (MQ=255)
            cATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGAtt                                                                                                                                       <  4:55135/149‑1 (MQ=255)
              tCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTcc                                                                                                                                     <  4:113531/149‑1 (MQ=255)
               cGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTc                                                                                                                                      >  7:246660/1‑147 (MQ=255)
                 aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGGACCCGTTCGCGAtaatcga                                                                                                                                  <  3:253575/149‑7 (MQ=255)
                 aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                  >  4:35163/1‑149 (MQ=255)
                   aaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagcca                                                                                                                                <  2:48825/149‑1 (MQ=255)
                   aaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                  >  6:287004/1‑147 (MQ=255)
                    aCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagccag                                                                                                                               >  3:6948/1‑149 (MQ=255)
                      gACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAgccaggc                                                                                                                             >  5:65142/1‑149 (MQ=255)
                        cATGCTCCTGTGCCGCCGCTTGCTCGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGCGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATAGCGAACCGGTCCGCGATTCCAgccaggcca                                                                                                                           >  2:92166/1‑149 (MQ=255)
                           gctgctGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACgg                                                                                                                                                 <  7:261166/124‑1 (MQ=255)
                           gctgctGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACgg                                                                                                                                                 >  8:261166/1‑124 (MQ=255)
                            ctgctgTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccg                                                                                                                       >  4:227902/1‑149 (MQ=255)
                              gctgTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGAGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgcc                                                                                                                     >  1:68185/1‑149 (MQ=255)
                                    ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGc                                                                                                               >  4:303250/1‑149 (MQ=255)
                                    ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGc                                                                                                               >  1:125427/1‑149 (MQ=255)
                                     cgccgcTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCt                                                                                                              >  8:258370/1‑149 (MQ=255)
                                       ccgcTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTg                                                                                                            <  2:68185/149‑1 (MQ=255)
                                         gcTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtc                                                                                                          <  6:64298/149‑1 (MQ=255)
                                          cTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCt                                                                                                                 >  7:30068/1‑141 (MQ=255)
                                          cTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCt                                                                                                                 <  8:30068/141‑1 (MQ=255)
                                          cTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcg                                                                                                         <  7:31137/149‑1 (MQ=255)
                                           ttGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        <  4:6948/149‑1 (MQ=255)
                                           ttGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGAGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        >  8:135920/1‑149 (MQ=255)
                                             gCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCTCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        <  8:246660/147‑1 (MQ=255)
                                             gCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        >  6:36209/1‑147 (MQ=255)
                                             gCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        >  6:1019/1‑147 (MQ=255)
                                              cTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtca                                                                                                      <  1:21486/148‑1 (MQ=255)
                                               tGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtcagt                                                                                                    >  5:206802/1‑149 (MQ=255)
                                                     aGGTCGACCACCAGGCGCTGAGCGCACAGTACATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCgg                                                                                              <  3:35163/149‑1 (MQ=255)
                                                       gTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        <  1:18547/137‑1 (MQ=255)
                                                       gTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                        >  2:18547/1‑137 (MQ=255)
                                                              accaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcgg                                                                                     <  7:134926/149‑1 (MQ=255)
                                                               ccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTTCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcgga                                                                                    <  2:125427/149‑1 (MQ=255)
                                                               ccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtca                                                                                                      <  6:130789/131‑1 (MQ=255)
                                                               ccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtca                                                                                                      >  5:130789/1‑131 (MQ=255)
                                                                           cgcACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGaaa                                                                        <  3:26851/149‑1 (MQ=255)
                                                                              acaGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAacac                                                                     <  3:303250/149‑1 (MQ=255)
                                                                              acaGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAacac                                                                     >  2:38849/1‑149 (MQ=255)
                                                                              acaGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAacac                                                                     >  2:34130/1‑149 (MQ=255)
                                                                                aGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTg                                                                   <  5:1019/149‑1 (MQ=255)
                                                                                     aTCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGa                                                              >  2:221061/1‑149 (MQ=255)
                                                                                                   cgTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGAtt                                                <  1:254761/149‑1 (MQ=255)
                                                                                                       ccAGAGACTGCCGAGCAGATACAGCAATCCCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTc                                             <  5:287004/148‑1 (MQ=255)
                                                                                                       ccAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGTCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAg                                                                <  1:260060/129‑1 (MQ=255)
                                                                                                       ccAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAg                                                                >  2:260060/1‑129 (MQ=255)
                                                                                                        cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                            >  1:313749/1‑148 (MQ=255)
                                                                                                        cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                            <  6:206802/148‑1 (MQ=255)
                                                                                                        cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                            <  5:36209/148‑1 (MQ=255)
                                                                                                        cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTc                                             <  5:302357/147‑1 (MQ=255)
                                                                                                        cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTc                                             >  6:302357/1‑147 (MQ=255)
                                                                                                            gaCTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATa                                       <  7:258370/149‑1 (MQ=255)
                                                                                                              atGCCGAGCATATACATCAATACCTCGATTATCGCGAACGTGTCCGCGATTCCAGCCAGGACATCCGCCGCCTGATTGTCGTCAGTCAGCGGCTGGATCTTATCCTGGCTGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATaaa                                     <  3:231242/148‑1 (MQ=255)
                                                                                                                 ccGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGcc                                  <  1:92166/149‑1 (MQ=255)
                                                                                                                 ccGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGc                                   >  3:46743/1‑148 (MQ=255)
                                                                                                                     gCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACg                               >  8:84247/1‑148 (MQ=255)
                                                                                                                        gATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAAc                           >  5:81079/1‑149 (MQ=255)
                                                                                                                         aTACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGaaa                            <  6:5950/147‑1 (MQ=255)
                                                                                                                         aTACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGaaa                            >  5:5950/1‑147 (MQ=255)
                                                                                                                          tACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTa                         <  2:264526/149‑1 (MQ=255)
                                                                                                                           aCAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        >  8:56055/1‑149 (MQ=255)
                                                                                                                           aCAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        >  8:51677/1‑149 (MQ=255)
                                                                                                                            cAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACGCTGCAGGCCCCCAAGATCGATTTTCCCAATAACGCCACGAAGCTAGc                       >  6:306713/1‑149 (MQ=255)
                                                                                                                            cAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGc                       <  1:221061/149‑1 (MQ=255)
                                                                                                                             aGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgc                                                                                                                      <  8:243519/53‑1 (MQ=255)
                                                                                                                                 aTACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGa                  >  3:22772/1‑149 (MQ=255)
                                                                                                                                    ccTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGGTCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCa               >  4:275070/1‑149 (MQ=255)
                                                                                                                                     cacGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGACAGCCCCCGCATGCTTCTCGTAAGTCACCGGCTGGATCGGATCCTGGCGGAAACACTGAAGGACCCCAAGAACGATATTCCCAATAAATCCAACAAAATATCCAAGAACAt              >  4:38491/3‑149 (MQ=255)
                                                                                                                                     cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAt              <  4:22772/149‑1 (MQ=255)
                                                                                                                                     cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAt              >  4:209074/1‑149 (MQ=255)
                                                                                                                                      tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAg                   <  8:287937/143‑1 (MQ=255)
                                                                                                                                      tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAg                   >  7:287937/1‑143 (MQ=255)
                                                                                                                                        gATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             <  6:65142/147‑1 (MQ=255)
                                                                                                                                        gATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             <  7:84247/147‑1 (MQ=255)
                                                                                                                                         ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                    >  7:6098/4‑139 (MQ=255)
                                                                                                                                         ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                    >  5:149139/4‑139 (MQ=255)
                                                                                                                                         ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                    <  6:149139/136‑1 (MQ=255)
                                                                                                                                         ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                    <  8:6098/136‑1 (MQ=255)
                                                                                                                                             tCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAAATAGCCAAGACCATTGAACg        >  7:108056/1‑147 (MQ=255)
                                                                                                                                             tCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACg        >  8:205846/1‑147 (MQ=255)
                                                                                                                                                 gAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTATCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGAc  >  1:61718/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTACGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGAC  >  NC_002947/5709228‑5709521

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: ‑