Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,767,351 A→T 8.1% M1174K (ATG→AAG) ‡ PP_5060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,767,3510AT 8.1% 207.8 / 7.9 74M1174K (ATG→AAG) ‡PP_5060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base T (2/4);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

GCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATG  >  NC_002947/5767207‑5767497
                                                                                                                                                 |                                                                                                                                                  
gCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGc                                                                                                                                                >  1:121547/1‑149 (MQ=255)
     aGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGa                                                                                                                                           <  8:200202/149‑1 (MQ=255)
     aGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggt                                                                                                                                           <  5:73457/149‑2 (MQ=255)
      ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtc                                                                                                                                          <  5:274981/149‑3 (MQ=255)
      ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTgg                                                                                                                                            <  2:86003/147‑1 (MQ=255)
      ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACACTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTgg                                                                                                                                            >  1:86003/1‑147 (MQ=255)
      ggCAACTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtc                                                                                                                                          <  4:245504/149‑3 (MQ=255)
       gCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaa                                                                                                                                          >  2:250602/1‑148 (MQ=255)
            tCACAGAGCATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCAAGGCCTGGAAAAACCg                                                                                                                                    <  5:19773/149‑1 (MQ=255)
             cacaGAACATCAGCAGTTGATGCGTTTGGGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTACGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAAcc                                                                                                                                     <  1:95190/148‑1 (MQ=255)
              acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGa                                                                                                                                  >  8:245444/1‑149 (MQ=255)
               caGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGAt                                                                                                                                 <  2:121547/149‑1 (MQ=255)
                aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGCCGCTGACAGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGCGTATGCGTGCGTTTGTGCCACTGCAGGGCCAGCTCCAACACGCACAATGCCTGCAAAAAGCGCAtc                                                                                                                                >  1:304265/1‑148 (MQ=255)
                aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATg                                                                                                                                >  6:1691/1‑149 (MQ=255)
                aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATg                                                                                                                                >  3:212411/1‑149 (MQ=255)
                   aCATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGtt                                                                                                                              <  4:91920/148‑1 (MQ=255)
                    cATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGa                                                                                                                            <  7:245444/149‑1 (MQ=255)
                             atgatgCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGgtcg                                                                                                                   >  7:86714/1‑149 (MQ=255)
                                    gTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACg                                                                                                            >  7:257441/1‑149 (MQ=255)
                                    gTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGACGTCTTACg                                                                                                            >  3:137192/1‑149 (MQ=255)
                                         tgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTAc                                                                                                             <  1:209451/143‑1 (MQ=255)
                                         tgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTAc                                                                                                             >  2:209451/1‑143 (MQ=255)
                                          gCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTcgc                                                                                                      >  5:268328/1‑149 (MQ=255)
                                           cGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTcgc                                                                                                      <  6:268328/148‑1 (MQ=255)
                                                  ccGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGc                                                                                              >  7:44998/1‑149 (MQ=255)
                                                    gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                            <  8:257441/149‑1 (MQ=255)
                                                    gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                            <  7:65532/149‑1 (MQ=255)
                                                      aGACGCTGACTGGCGAGCGCATGGCGGAAGATCGTGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGTAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                          <  1:243576/149‑1 (MQ=255)
                                                      aGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                          <  6:188002/149‑1 (MQ=255)
                                                      aGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                          <  1:250602/149‑1 (MQ=255)
                                                          gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATa                                                                                                  >  6:292506/1‑137 (MQ=255)
                                                          gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATa                                                                                                  <  5:292506/137‑1 (MQ=255)
                                                          gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACg                                                                                      >  6:152388/1‑149 (MQ=255)
                                                           cTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGACGTCTTACGGCTCGCGATAATGCTCGAGACgg                                                                                     <  4:137192/149‑1 (MQ=255)
                                                                tGGCTAGCGCATTGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCaga                                                                                <  7:292131/149‑1 (MQ=255)
                                                                        gcATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTAtt                                                                        <  4:320578/149‑1 (MQ=255)
                                                                        gcATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTAtt                                                                        >  3:196554/1‑149 (MQ=255)
                                                                         cATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTg                                                                       <  4:212411/149‑1 (MQ=255)
                                                                           tGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGAc                                                                     <  5:152388/149‑1 (MQ=255)
                                                                               ggAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTAAGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGt                                                                   >  4:74/1‑147 (MQ=255)
                                                                                gAAGATCATGTAGTGCTGGGTGGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAg                                                                >  5:330893/1‑149 (MQ=255)
                                                                                gAAGATCATGTAGTGCTGGGTGGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAg                                                                >  5:335345/1‑149 (MQ=255)
                                                                                       aTGTTGTGGTGGGTAGTGGTGCGTTTGTGCCAGTGCAGGGCCAGCTCTAACAGGCCCATGGACTTGAAAAACCGGATGTTGATCGTGGTCGTGTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCAGtt                                                          <  2:304265/148‑1 (MQ=255)
                                                                                        tGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTg                                                         >  1:204323/1‑148 (MQ=255)
                                                                                                tGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACtt                                                          >  4:195296/1‑139 (MQ=255)
                                                                                                tGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACtt                                                          <  3:195296/139‑1 (MQ=255)
                                                                                                 gggTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgagg                                               >  4:38146/1‑147 (MQ=255)
                                                                                                  ggTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgatgt                                              >  5:244663/1‑149 (MQ=255)
                                                                                                  ggTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgat                                                <  2:204323/147‑1 (MQ=255)
                                                                                                          gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACTGTCAGATTTGTATTGACGTAAg                                                                <  2:317446/124‑1 (MQ=255)
                                                                                                          gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTTGTATTGACGTAAg                                                                >  1:317446/1‑124 (MQ=255)
                                                                                                          gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGAt                                       >  2:33304/1‑149 (MQ=255)
                                                                                                           tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCaga                                                                                <  8:159882/107‑1 (MQ=255)
                                                                                                           tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCaga                                                                                >  7:159882/1‑107 (MQ=255)
                                                                                                           tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATg                                      <  4:196554/149‑1 (MQ=255)
                                                                                                             ccTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGc                                    >  5:173575/1‑149 (MQ=255)
                                                                                                                   tgCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACa                               <  3:74/148‑1 (MQ=255)
                                                                                                                   tgCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCACGATAAGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAgg                             <  8:86714/149‑1 (MQ=255)
                                                                                                                       aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGa                          >  7:226727/1‑149 (MQ=255)
                                                                                                                       aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCAAACAGGCGa                          >  5:137567/1‑149 (MQ=255)
                                                                                                                       aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGACCAACAGGCGa                          <  5:1691/149‑1 (MQ=255)
                                                                                                                       aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGAGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGa                          >  1:87321/1‑149 (MQ=255)
                                                                                                                              ggCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgct                   <  2:87321/149‑1 (MQ=255)
                                                                                                                                 cAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCat                                                      >  4:56430/1‑111 (MQ=255)
                                                                                                                                 cAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCat                                                      <  3:56430/111‑1 (MQ=255)
                                                                                                                                 cAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgctggc                >  5:200902/1‑149 (MQ=255)
                                                                                                                                   gtttttccAGGCCATGGGCCTGGAAAAATCGGATGTTGATCGTGGTCGTCTTACGGCTCGCAATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgctggctg              >  3:287415/7‑149 (MQ=255)
                                                                                                                                     tCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTTGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGga            <  6:173575/149‑1 (MQ=255)
                                                                                                                                       gACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGga            <  8:226727/147‑1 (MQ=255)
                                                                                                                                        tccAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGgaga          >  4:53114/2‑148 (MQ=255)
                                                                                                                                          cAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACAt       >  5:226662/1‑149 (MQ=255)
                                                                                                                                              cccATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtgat   <  6:335345/149‑1 (MQ=255)
                                                                                                                                              cccATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtgat   <  6:330893/149‑1 (MQ=255)
                                                                                                                                               ccATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCTATAATGCTCGAGACGGTCAGAGTGGTATTGAGTTTAGGCACTTGGATGTCGATGTCGGTGATGGCCTTCAGGCGATCGCGGTGGCTGGAGACATtgatg  >  6:86865/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATG  >  NC_002947/5767207‑5767497

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: