Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,844,774 G→C 8.3% V466L (GTG→CTG) ‡ PP_5122 → GMC family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,844,7740GC 8.3% 215.8 / 7.1 70V466L (GTG→CTG) ‡PP_5122GMC family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base C (4/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

ACACGTCCAGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCGGTGTAT  >  NC_002947/5844632‑5844920
                                                                                                                                                 |                                                                                                                                                   
acacGTCCAGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATCTGATgg                                                                                                                                               >  3:200845/1‑149 (MQ=255)
  acGTCCAGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGggc                                                                                                                                             <  4:196766/149‑1 (MQ=255)
    gTCCAGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGggcgg                                                                                                                                           >  1:211457/1‑149 (MQ=255)
        aGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCCGCCCATCACATGggcgctgc                                                                                                                                        <  3:353757/149‑5 (MQ=255)
        aGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTgc                                                                                                                                       <  5:231047/149‑1 (MQ=255)
        aGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTgc                                                                                                                                       <  6:345860/149‑1 (MQ=255)
        aGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTgc                                                                                                                                       <  7:166899/149‑1 (MQ=255)
        aGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTgc                                                                                                                                       <  8:185845/149‑1 (MQ=255)
               cGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCg                                                                                                                                >  2:4575/1‑149 (MQ=255)
               cGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATCACATGggcgctgc                                                                                                                                        >  2:248594/1‑137 (MQ=255)
               cGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATCACATGggcgctgc                                                                                                                                        <  1:248594/141‑5 (MQ=255)
                         aCCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGAccc                                                                                                                      >  6:334971/1‑149 (MQ=255)
                           ccATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCcgc                                                                                                                    >  1:255664/1‑149 (MQ=255)
                           ccATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCcgc                                                                                                                    >  3:55248/1‑149 (MQ=255)
                                   ccccGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCAcgcg                                                                                                            >  2:6018/1‑149 (MQ=255)
                                    cccGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCgg                                                                                                           >  5:354309/1‑149 (MQ=255)
                                     ccGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGgt                                                                                                          >  8:59467/1‑149 (MQ=255)
                                     ccGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGgt                                                                                                          <  3:362789/149‑1 (MQ=255)
                                     ccGGTGCATAGCCACGCCAGCGCCGGCTCAAGCCTGGAACAGGCCCTGACCATGAACGACAAACAGCGCCTGGATCCCTCATCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCAcgcg                                                                                                            <  1:299381/149‑1 (MQ=255)
                                           cATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCAcgcg                                                                                                            >  7:284228/1‑141 (MQ=255)
                                           cATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCAcgcg                                                                                                            <  8:284228/141‑1 (MQ=255)
                                           cATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGa                                                                                                    >  7:240820/1‑149 (MQ=255)
                                            aTAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGa                                                                                                    >  5:155303/1‑148 (MQ=255)
                                             tAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATTGGCGGCTGCGCCCTCGGCGAATACGCGCGCCACGCGGTGTGCGACa                                                                                                  <  6:152935/149‑1 (MQ=255)
                                               gCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGc                                                                                                >  8:254016/1‑149 (MQ=255)
                                                 gACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCt                                                                                              <  2:50383/149‑1 (MQ=255)
                                                      cAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGgcc                                                                                         >  6:132741/1‑149 (MQ=255)
                                                        gcgcCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGgccgc                                                                                       <  7:254016/149‑1 (MQ=255)
                                                           ccgccTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCAt                                                                                    <  6:253855/149‑1 (MQ=255)
                                                               ctctAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTgcg                                                                                                                                      <  1:229416/95‑1 (MQ=255)
                                                               ctctAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTgcg                                                                                                                                      >  2:229416/1‑95 (MQ=255)
                                                                  tAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATggg                                                                                                                                              <  3:123534/84‑1 (MQ=255)
                                                                  tAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATggg                                                                                                                                              >  4:123534/1‑84 (MQ=255)
                                                                      cTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTgg                                                                         >  3:321682/1‑148 (MQ=255)
                                                                       tGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGa                                                                        >  2:103221/1‑148 (MQ=255)
                                                                          aaCAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAAc                                                                    <  2:255664/149‑1 (MQ=255)
                                                                               gCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTc                                                               <  5:334971/149‑1 (MQ=255)
                                                                                 ccTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGa                                                             <  2:211457/149‑1 (MQ=255)
                                                                                 ccTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGa                                                             >  4:116969/1‑149 (MQ=255)
                                                                                 ccTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGa                                                             >  5:80101/1‑149 (MQ=255)
                                                                                     aCCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCa                                                         <  6:354309/149‑1 (MQ=255)
                                                                                       cATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCCCGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCac                                                        >  4:315122/1‑148 (MQ=255)
                                                                                       cATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCacg                                                       >  3:42988/1‑149 (MQ=255)
                                                                                       cATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCac                                                        >  1:32758/1‑148 (MQ=255)
                                                                                        aTGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGaaa                                                                      <  4:226060/133‑1 (MQ=255)
                                                                                        aTGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGaaa                                                                      >  3:226060/1‑133 (MQ=255)
                                                                                         tGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCCAGCTGGAAAACCTGTCGATTCacga                                                      <  1:4575/149‑1 (MQ=255)
                                                                                            tCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGt                                                  <  6:186719/149‑1 (MQ=255)
                                                                                               aCAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCGGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCa                                               <  5:206117/149‑1 (MQ=255)
                                                                                                cAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCAc                                              >  5:134251/1‑149 (MQ=255)
                                                                                                 aaaCTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACt                                             <  6:80101/149‑1 (MQ=255)
                                                                                                       cgcCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTc                                         <  1:6018/147‑1 (MQ=255)
                                                                                                        gcCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACCGTTccc                                       >  8:330660/1‑148 (MQ=255)
                                                                                                          cTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCccac                                    <  6:134251/149‑1 (MQ=255)
                                                                                                            ggAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCccacca                                  <  5:10021/149‑1 (MQ=255)
                                                                                                                   t‑‑tCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCgg                           >  2:219357/1‑149 (MQ=255)
                                                                                                                           aCGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCaa                     >  6:300838/1‑149 (MQ=255)
                                                                                                                              cgcTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAAccc                  <  8:240820/149‑1 (MQ=255)
                                                                                                                               gcTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGgccgc                                                                                       >  3:39614/1‑80 (MQ=255)
                                                                                                                               gcTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGgccgc                                                                                       <  4:39614/80‑1 (MQ=255)
                                                                                                                                  tGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCTCGCGGTGTGCGACTGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCGGCACCAGCATCGGCGCCTACGGGCaa              <  4:42988/149‑1 (MQ=255)
                                                                                                                                   gggCAGCGCCCATGTGATGGGCGGCTGCGCCGTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTAGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCaa              <  5:132741/148‑1 (MQ=255)
                                                                                                                                       agccgcCCATCACATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGt          >  3:363228/4‑149 (MQ=255)
                                                                                                                                       agccgcCCATCACATGGGCGGCTGCGCCCACGGCGAAGACCCGCGCCACCCGGTGTGCGACCCCCTGGGACGCCATCACCAGCTGGAAAACCTGTCGATTCACGACCGTTCACTGTTCCCCACACGCATCCGCGCCAACCCCCAATTGt          >  2:284465/4‑149 (MQ=255)
                                                                                                                                       aGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTc         <  1:103221/149‑1 (MQ=255)
                                                                                                                                        gcgcCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCg        <  3:315122/149‑1 (MQ=255)
                                                                                                                                          cgcCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCgg       <  2:32758/149‑1 (MQ=255)
                                                                                                                                            cccATGTGATGGTCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCGgtg     <  4:244813/149‑1 (MQ=255)
                                                                                                                                             ccATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCGgtgt    >  3:200914/1‑149 (MQ=255)
                                                                                                                                               aTGTGAGGCGCGGCTGCGCCCGCGGCGAAGTCCCGCGCCACGCGGTGTGCGACAGCCTGTGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCGGTGTAt  <  4:200845/149‑1 (MQ=255)
                                                                                                                                                 gtgATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTc         >  5:282762/1‑140 (MQ=255)
                                                                                                                                                 gtgATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTCTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTc         <  6:282762/140‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACACGTCCAGTTTGCCGCTGGCGCCACCCATGTAACCCCGGTGCATAGCGACGCCAGCGCCGCCTCTAGCCTGGAACAGGCCCTGACCATGATCGACAAACTGCGCCTGGAGCCCT‑‑TCCGCACGCGCTTGGGCAGCGCCCATGTGATGGGCGGCTGCGCCCTCGGCGAAGACCCGCGCCACGCGGTGTGCGACAGCCTGGGCCGCCATCACCAGCTGGAAAACCTGTCGATTCACGACGGTTCACTGTTCCCCACCAGCATCGGCGCCAACCCGCAATTGTCGGTGTAT  >  NC_002947/5844632‑5844920

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: