Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 166,579 T→C 7.7% T110A (ACC→GCC)  pntAA ← pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947166,5790TC 7.7% 220.3 / 6.8 78T110A (ACC→GCC) pntAApyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base C (2/4);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GCAGCGGTCATCAGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGTACCAATTGGGC  >  NC_002947/166431‑166711
                                                                                                                                                    |                                                                                                                                     
gCAGCGGTCATCAGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGt                                                                                                                                       <  6:91524/149‑1 (MQ=255)
     ggTCATCAGCATTGGCATGAAGCGTGCGTAGTGATGGGCTGCCAGCAACACGGCCATGTAGCCGGCTATGTTGGCCAGCGACGACAGCACGTCCACGCTCTGCGCCCGCGAGGTTAGCGGCGCGCCTTCCAGGGAGAAGGCGGTAATGc                                                                                                                                  >  6:197597/1‑149 (MQ=255)
      gTCATCAGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGcc                                                                                                                                 >  2:76830/1‑149 (MQ=255)
         atcaGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTTCAGGGCGAAGGCGGTAATGCcgcg                                                                                                                              >  1:326677/1‑149 (MQ=255)
         atcaGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCcgcg                                                                                                                              >  4:151946/1‑149 (MQ=255)
          tcaGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGt                                                                                                                             >  3:16519/1‑149 (MQ=255)
            aGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTc                                                                                                                           <  3:183439/149‑1 (MQ=255)
             gCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCg                                                                                                                          <  1:316271/149‑1 (MQ=255)
              cATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCgg                                                                                                                         <  8:269920/149‑1 (MQ=255)
              cATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAACGCGGCATTACCGCCTTCgc                                                                                                                         <  5:168234/149‑2 (MQ=255)
               aTTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGc                                                                                                                       <  7:184730/149‑1 (MQ=255)
                      tGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGCCGAACGCGGCATTACCGCCTTCGccctggaag                                                                                                                <  3:277665/149‑7 (MQ=38)
                        aaGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGcc                                                                                                              <  3:35059/149‑1 (MQ=255)
                             tGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGCAACTCCCCCCCCGAGTTTCCCGCCCCGCATTCCAGGGCAAAGGCGTTAATGCCCCTTACGGCCAAATTGCCTATCa                                                                                                         >  4:225536/1‑149 (MQ=255)
                               ggTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAg                                                                                                        <  6:182500/148‑1 (MQ=255)
                                   gTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCa                                                                                                    <  6:284395/148‑1 (MQ=255)
                                        gggCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGtt                                                                                               >  5:109230/1‑148 (MQ=255)
                                         ggCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAAGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGtt                                                                                                      <  2:317688/140‑1 (MQ=255)
                                         ggCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAAGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGtt                                                                                                      >  1:317688/1‑140 (MQ=255)
                                           cTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAg                                                                                           <  7:247017/149‑1 (MQ=255)
                                           cTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAg                                                                                           >  8:247017/1‑149 (MQ=255)
                                            tGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAgg                                                                                          >  8:211655/1‑149 (MQ=255)
                                              ccAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGTCAGCACGTGCTGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGGCGCGTTCGTCCATCTTGCCGATCAGTTGACCGTTGAAgg                                                                                          <  5:197597/147‑1 (MQ=255)
                                              ccAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAggg                                                                                         <  6:141170/148‑1 (MQ=255)
                                               cAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGt                                                                                       >  4:186180/1‑149 (MQ=255)
                                                aGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGtt                                                                                      >  3:244774/1‑149 (MQ=255)
                                                aGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGtt                                                                                      <  5:327965/149‑1 (MQ=255)
                                                 gCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGtt                                                                                      >  3:61386/1‑148 (MQ=255)
                                                 gCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGtt                                                                                      >  5:208189/1‑148 (MQ=255)
                                                  cAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAACGCGGCATTACCGCCTTCgcc                                                                                                                        <  2:248618/114‑3 (MQ=255)
                                                  cAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAACGCGGCATTACCGCCTTCgcc                                                                                                                        >  1:248618/1‑112 (MQ=255)
                                                   aaCACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAg                                                                                           <  3:241831/141‑1 (MQ=255)
                                                   aaCACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAg                                                                                           >  4:241831/1‑141 (MQ=255)
                                                   aaCACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAg                                                                                   <  3:186180/149‑1 (MQ=255)
                                                    acacGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGa                                                                                    <  5:249391/147‑1 (MQ=255)
                                                    acacGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGa                                                                                    >  6:249391/1‑147 (MQ=255)
                                                       cGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCtt                                                                                                                 >  6:33178/1‑115 (MQ=255)
                                                       cGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCtt                                                                                                                 <  5:33178/115‑1 (MQ=255)
                                                          ccTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGccc                                                                            <  4:16519/149‑1 (MQ=255)
                                                          ccTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGcc                                                                             >  2:6999/1‑148 (MQ=255)
                                                          ccTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGcc                                                                             >  2:12553/1‑148 (MQ=255)
                                                            ttGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCccac                                                                          >  4:199108/1‑149 (MQ=255)
                                                              gTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCccac                                                                          <  1:76830/148‑1 (MQ=255)
                                                              gTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGa                                                                                    >  6:282384/1‑137 (MQ=255)
                                                              gTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGa                                                                                    <  5:282384/137‑1 (MQ=255)
                                                              gTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCccac                                                                          <  1:12553/147‑1 (MQ=255)
                                                              gTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCccac                                                                          <  1:6999/147‑1 (MQ=255)
                                                                aGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAg                                                                       >  1:55614/1‑148 (MQ=255)
                                                                 gCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGa                                                                      <  4:61386/148‑1 (MQ=255)
                                                                       gATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGc                                                               >  7:229778/1‑149 (MQ=255)
                                                                        aTGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGcc                                                              <  3:127953/149‑1 (MQ=255)
                                                                             ggCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAACGCGGCATTACCGCCTTCGccc                                                                                                                      >  2:187477/1‑88 (MQ=37)
                                                                             ggCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAACGCGGCATTACCGCCTTCGccc                                                                                                                      <  1:187477/88‑1 (MQ=37)
                                                                             ggCCTGCGACGACAGCACGGCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTCCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTg                                                          <  6:258487/148‑1 (MQ=255)
                                                                              gCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGt                                                         <  3:151946/148‑1 (MQ=255)
                                                                                     acgacAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGtt                                                        >  8:178349/1‑142 (MQ=255)
                                                                                     acgacAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGtt                                                        <  7:178349/142‑1 (MQ=255)
                                                                                     acgacAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGGTGATAAgg                                                 >  6:287056/1‑149 (MQ=255)
                                                                                        acAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGc                                               >  4:94410/1‑148 (MQ=255)
                                                                                           gCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCTCCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGt                                           >  5:106973/1‑149 (MQ=255)
                                                                                           gCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGt                                           >  7:35795/1‑149 (MQ=255)
                                                                                            cACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAgggg                                                                                        <  4:42131/103‑1 (MQ=255)
                                                                                             aCGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTc                                         >  6:238786/1‑149 (MQ=255)
                                                                                             aCGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTc                                         >  5:139410/1‑149 (MQ=255)
                                                                                                  cAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGaa                                                                     >  8:136816/1‑116 (MQ=255)
                                                                                                  cAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTgg                                     <  7:211655/148‑1 (MQ=255)
                                                                                                        ctGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTgg                              >  3:313036/1‑149 (MQ=255)
                                                                                                         tGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGc                                                               >  4:77470/1‑115 (MQ=255)
                                                                                                         tGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGc                                                               <  3:77470/115‑1 (MQ=255)
                                                                                                           cgcCCGCGAGGTTCGCGGCGCGGCTTTCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGCAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTggg                             <  2:326677/147‑1 (MQ=255)
                                                                                                           cgcCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTggggg                           >  6:240671/1‑149 (MQ=255)
                                                                                                           cgcCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTggg                             <  8:35795/147‑1 (MQ=255)
                                                                                                            gcCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAACGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGc                          <  6:208189/149‑1 (MQ=255)
                                                                                                                           ggcgcggcTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGTAc           <  4:244774/149‑1 (MQ=255)
                                                                                                                           ggcgcggcTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGTAc           <  5:240671/149‑1 (MQ=255)
                                                                                                                                  cTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGTACCAATTgg    <  2:55614/149‑1 (MQ=255)
                                                                                                                                     ccAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGTACCAATTGGgc  <  5:238786/148‑1 (MQ=255)
                                                                                                                                           cGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGt             <  7:282581/131‑1 (MQ=255)
                                                                                                                                           cGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGt             >  8:282581/1‑131 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
GCAGCGGTCATCAGCATTGGCATGAAGCGTGGGTAGTGATGGGCTGCCAGCAACACGGCCTTGTAGCCGGCGATGTTGGCCTGCGACGACAGCACGTCCAGGCTCTGCGCCCGCGAGGTTCGCGGCGCGGCTTCCAGGGCGAAGGCGGTAATGCCGCGTTCGGCCATCTTGCCGATCAGTTCACCGTTGAAGGGGTTGAGCATGCCCACCAGAAGGCTGCCGCTGTTGATCAGGGCCAGTTCCTGGTCGTTGGGGGCGACCACTTTGAGTACCAATTGGGC  >  NC_002947/166431‑166711

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: