Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,709,381 C→T 7.2% V1133V (GTC→GTT)  PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,709,3810CT 7.2% 244.6 / 5.9 83V1133V (GTC→GTT) PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/35);  new base T (3/3);  total (45/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 = 9.12e-01

AGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACG  >  NC_002947/5709237‑5709522
                                                                                                                                                |                                                                                                                                             
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)
                   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                                                                                                         >  2:18547/1‑137 (MQ=255)
                                              gTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtcgt                                                                                                         <  1:18547/137‑1 (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                                                                      >  2:38849/1‑149 (MQ=255)
                                                                     acaGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAacac                                                                      <  3:303250/149‑1 (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                                             <  5:36209/148‑1 (MQ=255)
                                                                                               cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                             >  1:313749/1‑148 (MQ=255)
                                                                                               cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                             <  6:206802/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              <  7:84247/147‑1 (MQ=255)
                                                                                                                               gATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt              <  6:65142/147‑1 (MQ=255)
                                                                                                                                ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                     >  5:149139/4‑139 (MQ=255)
                                                                                                                                ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                     <  8:6098/136‑1 (MQ=255)
                                                                                                                                ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                     >  7:6098/4‑139 (MQ=255)
                                                                                                                                ggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                     <  6:149139/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)
                                                                                                                                          aCGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGTTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACg  >  3:173464/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                             
AGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACG  >  NC_002947/5709237‑5709522

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