Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,749,313 C→T 6.2% G339G (GGC→GGT PP_4203 → ubiquinone oxidoreductase electron transfer flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,749,3130CT 6.2% 194.7 / 2.0 64G339G (GGC→GGTPP_4203ubiquinone oxidoreductase electron transfer flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/32)
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.67e-01

CGCCTCAAGCACCACCCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAGCGCCTTCAAGGCCAGCTGGC  >  NC_002947/4749167‑4749461
                                                                                                                                                  |                                                                                                                                                     
cgcCTCAAGCACCACCCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCAc                                                                                                                                                     <  1:195759/149‑1 (MQ=255)
  ccTCAAGCACCACCCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCa                                                                                                                                                   <  2:123225/149‑1 (MQ=255)
          accaccCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAAAGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTg                                                                                                                                           >  3:45178/1‑149 (MQ=255)
           ccaccCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTgg                                                                                                                                          <  6:159840/149‑1 (MQ=255)
            caccCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGc                                                                                                                                         >  4:310884/1‑149 (MQ=255)
              cccGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGcc                                                                                                                                        >  2:1156/1‑148 (MQ=255)
                     aTCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCaa                                                                                                                                <  6:142258/149‑1 (MQ=255)
                      tCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGTACCATGAACGTGGCCAAGATCAAg                                                                                                                               >  4:238242/1‑149 (MQ=255)
                         gCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAggg                                                                                                                             >  7:12916/1‑148 (MQ=255)
                           cAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAg                                                                                                                          <  1:1620/149‑1 (MQ=255)
                             gTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGc                                                                                                                         >  6:209891/1‑148 (MQ=255)
                                ccTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGTCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATa                                                                                                                    <  4:236808/149‑1 (MQ=255)
                                 cTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATAc                                                                                                                   >  2:246382/1‑149 (MQ=255)
                                 cTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATAc                                                                                                                   >  1:231534/1‑149 (MQ=255)
                                     aaGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAgg                                                                                                                              >  6:70464/1‑135 (MQ=255)
                                     aaGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACcgcg                                                                                                               <  5:209891/149‑1 (MQ=255)
                                      aGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGa                                                                                                              <  2:79752/149‑1 (MQ=255)
                                       ggcggcAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGAt                                                                                                             >  8:171450/1‑149 (MQ=255)
                                        gcggcAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATg                                                                                                            <  7:171450/149‑1 (MQ=255)
                                          ggcAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGaa                                                                                                          <  1:20344/149‑1 (MQ=255)
                                           gcAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAg                                                                                                         <  2:231534/149‑1 (MQ=255)
                                              aGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTcc                                                                                                      >  6:240870/1‑149 (MQ=255)
                                                 gcATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGc                                                                                                   <  5:240870/149‑1 (MQ=255)
                                                  cATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAgg                                                                                                                              >  6:88461/1‑122 (MQ=255)
                                                  cATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAgg                                                                                                                              <  5:88461/122‑1 (MQ=255)
                                                           cGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTcgccg                                                                                         >  2:88827/1‑149 (MQ=255)
                                                           cGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTcgcc                                                                                          <  7:166490/148‑1 (MQ=255)
                                                               gcCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGccgccg                                                                                      >  8:57814/1‑148 (MQ=255)
                                                                 ccGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAg                                                                                   <  1:246382/149‑1 (MQ=255)
                                                                             cAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTc                                                                                             >  1:81212/1‑127 (MQ=255)
                                                                             cAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTc                                                                                             <  2:81212/127‑1 (MQ=255)
                                                                              aaGGGCGGCATCATCTCGCTGCCGAAAATTGTCTTCAATGGTGGCTCGCTGTTCGGTGGCGACGTGGTCACCATGAACGTGGCCAAGATCAAGGGCACCCATACCGCGATGAAGTCCGGCATGCTCGCCTCCGTTGCAGTGGCGGACGc                                                                      <  4:139043/149‑1 (MQ=255)
                                                                                gggCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGTACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCtt                                                                    <  3:238242/149‑1 (MQ=255)
                                                                                   cggcATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGa                                                                 <  8:12916/149‑1 (MQ=255)
                                                                                      cATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCg                                                              >  2:39909/1‑149 (MQ=255)
                                                                                        tCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATCTCCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGc                                                             <  1:88827/149‑1 (MQ=255)
                                                                                         taaCTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCg                                                           <  3:127271/148‑1 (MQ=255)
                                                                                           aCTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCgg                                                          >  2:87327/1‑148 (MQ=255)
                                                                                           aCTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCg                                                           >  4:127271/1‑147 (MQ=255)
                                                                                            cTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTa                                                        <  1:87327/149‑1 (MQ=255)
                                                                                                          tGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGTACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGccgccg                                                                                      <  8:287556/105‑1 (MQ=255)
                                                                                                          tGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGTACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGccgccg                                                                                      >  7:287556/1‑105 (MQ=255)
                                                                                                               ttCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAAtt                                     >  5:156517/1‑149 (MQ=255)
                                                                                                                tCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTg                                    <  7:236436/149‑1 (MQ=255)
                                                                                                                 cAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTg                                    >  8:236436/1‑148 (MQ=255)
                                                                                                                  aaTGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAgg                                                 >  5:279475/1‑134 (MQ=255)
                                                                                                                  aaTGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAgg                                                 <  6:279475/134‑1 (MQ=255)
                                                                                                                         gcgcgcTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCTGTAGCGAAGGCGGCGATCAATTGAACAGCTAt                           <  1:39909/149‑1 (MQ=255)
                                                                                                                         gcgcgcTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCGAt                           <  3:310884/149‑1 (MQ=255)
                                                                                                                         gcgcgcTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGACGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTAt                           >  6:116289/1‑149 (MQ=255)
                                                                                                                               tGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCa                                        >  6:165764/1‑130 (MQ=255)
                                                                                                                               tGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCa                                        <  5:165764/130‑1 (MQ=255)
                                                                                                                               tGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAAtt                                     >  2:12277/1‑133 (MQ=255)
                                                                                                                               tGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAAtt                                     <  1:12277/133‑1 (MQ=255)
                                                                                                                                    ggTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAgc                     <  1:319789/144‑1 (MQ=255)
                                                                                                                                    ggTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAgc                     >  2:319789/1‑144 (MQ=255)
                                                                                                                                        gCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAGCGCCTTCAAg            >  6:311807/1‑149 (MQ=255)
                                                                                                                                            ccTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTc                        >  1:59070/1‑133 (MQ=255)
                                                                                                                                            ccTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTc                        <  2:59070/133‑1 (MQ=255)
                                                                                                                                             cTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATAACGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATAGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTAAGCGCCTTCAAGGCCAg       >  6:82034/1‑149 (MQ=255)
                                                                                                                                               gggCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAAt                                      <  7:51425/116‑1 (MQ=255)
                                                                                                                                               gggCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAAt                                      >  8:51425/1‑116 (MQ=255)
                                                                                                                                               gggCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAGCGCCTTCAAGGCCAgct     >  7:99053/1‑149 (MQ=255)
                                                                                                                                                ggCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAATCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAGCGCCTTCAAGGCCAgctg    >  4:271939/1‑149 (MQ=255)
                                                                                                                                                  cACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAGCGCCTTCAAGGCCAgctggc  >  2:118685/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGCCTCAAGCACCACCCGGTGATCAGCCAGTACCTCGAAGGCGGCAAGCGCATCAGCTACGGCGCCCGCGCCTTGGCCAAGGGCGGCATCAACTCGCTGCCGAAAATGGTCTTCAATGGTGGCGCGCTGATCGGTTGCGACCTGGGCACCATGAACGTGGCCAAGATCAAGGGCAGCCATACCGCGATGAAGTCCGGCATGCTCGCCGCCGAAGCAGTGGCCGACGCTTTGATCGCCGGTAGCGAAGGCGGCGATCAATTGAACAGCTATGTCAGCGCCTTCAAGGCCAGCTGGC  >  NC_002947/4749167‑4749461

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: