Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,507 A→T 10.6% I249N (ATT→AAT)  purA ← adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,558,5070AT10.6% 181.4 / 11.2 66I249N (ATT→AAT) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base T (3/4);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

GAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGCATCATC  >  NC_002947/5558359‑5558653
                                                                                                                                                     |                                                                                                                                                  
gAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGaaa                                                                                                                                                    >  8:138826/1‑149 (MQ=255)
      gTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccac                                                                                                                                              >  7:225552/1‑149 (MQ=255)
       tGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccac                                                                                                                                              >  1:296924/1‑149 (MQ=255)
         gCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCgg                                                                                                                                           <  4:242823/149‑1 (MQ=255)
            aCGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCAGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggttc                                                                                                                                        <  8:20836/149‑4 (MQ=255)
             cGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggttcc                                                                                                                                       <  6:286322/149‑5 (MQ=255)
             cGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggttcc                                                                                                                                       <  2:298685/149‑5 (MQ=255)
                   gCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGtt                                                                                                                                 >  5:23552/1‑149 (MQ=255)
                   gCCAAGGGAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGAAGTGgtgt                                                                                                                                  <  8:36848/148‑1 (MQ=255)
                    cctatgTAGCACCGGTCTCGTCGAACAGTTCGGTCGTGAACGGAGAGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGATGATGTAGTCGAGGTACATCGGGCCGACGCCGTAACCGGTGGAAATACCACCGGCAGTGGTGTTg                                                                                                                                <  8:141122/144‑1 (MQ=255)
                    ccAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTg                                                                                                                                <  2:259314/149‑1 (MQ=255)
                    ccAAGGTAGCACCGGTCGCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTTGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTg                                                                                                                                <  2:259318/149‑1 (MQ=255)
                     cAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACGGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTg                                                                                                                                >  1:85683/1‑148 (MQ=255)
                     cAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTg                                                                                                                                >  2:257420/1‑148 (MQ=255)
                       aGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAg                                                                                                                             <  2:32342/149‑1 (MQ=255)
                        ggTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGc                                                                                                                            <  6:156254/149‑1 (MQ=255)
                            gCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTgg                                                                                                                         <  7:138826/148‑1 (MQ=255)
                              aCCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAgtggt                                                                                                                                    >  7:60804/1‑135 (MQ=255)
                              aCCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAgtggt                                                                                                                                    <  8:60804/135‑1 (MQ=255)
                                 ggTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACg                                                                                                                    >  6:14757/1‑148 (MQ=255)
                                   tctcGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTAc                                                                                                                 >  1:175659/1‑149 (MQ=255)
                                   tctcGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGCGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTAc                                                                                                                 >  3:151545/1‑149 (MQ=255)
                                     tcgtcgAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACgg                                                                                                               >  2:332847/1‑149 (MQ=255)
                                          gAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAg                                                                                                           <  8:225552/148‑1 (MQ=255)
                                              aGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGcc                                                                                                      <  1:257420/149‑1 (MQ=255)
                                                    ggtcggGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGt                                                                                                 <  6:23552/149‑1 (MQ=255)
                                                        gggAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGCGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGAt                                                                                             <  4:151545/149‑1 (MQ=255)
                                                         ggAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATg                                                                                            <  5:14757/149‑1 (MQ=255)
                                                         ggAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATg                                                                                            <  2:175659/149‑1 (MQ=255)
                                                            aCGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTcc                                                                                         <  2:296924/149‑1 (MQ=255)
                                                                acgggAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGt                                                                                                             >  5:229937/4‑125 (MQ=255)
                                                                acgggAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGt                                                                                                             <  6:229937/122‑1 (MQ=255)
                                                                       cccACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGcc                                                                              >  2:174295/1‑149 (MQ=255)
                                                                        ccACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGccc                                                                             >  5:207944/1‑149 (MQ=255)
                                                                        ccACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGccc                                                                             >  8:350060/1‑149 (MQ=255)
                                                                          aCGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTg                                                                           >  5:62568/1‑149 (MQ=255)
                                                                           cgcgAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTgg                                                                          <  8:278748/149‑1 (MQ=255)
                                                                                   gtgtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTc                                                                  <  2:189675/149‑1 (MQ=255)
                                                                                     gtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCg                                                                 >  7:312596/1‑148 (MQ=255)
                                                                                     gtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCg                                                                 >  6:168163/1‑148 (MQ=255)
                                                                                      tAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGaa                                                               <  1:174295/149‑1 (MQ=255)
                                                                                      tAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGaa                                                               <  5:168163/149‑1 (MQ=255)
                                                                                              ggTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACatgat                                                         <  8:312596/147‑1 (MQ=255)
                                                                                               gTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAAc                                                              <  4:198277/141‑1 (MQ=255)
                                                                                               gTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAAc                                                              >  3:198277/1‑141 (MQ=255)
                                                                                                  aTACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                                   >  3:336268/1‑149 (MQ=255)
                                                                                                  aTACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                                   >  4:277879/1‑149 (MQ=255)
                                                                                                  aTACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCt                                                    <  6:207944/148‑1 (MQ=255)
                                                                                                      cGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGAt                                                                                             <  5:318522/103‑1 (MQ=255)
                                                                                                      cGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGAt                                                                                             >  6:318522/1‑103 (MQ=255)
                                                                                                      cGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTc                                                      <  7:35294/142‑1 (MQ=255)
                                                                                                      cGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTc                                                      >  8:35294/1‑142 (MQ=255)
                                                                                                      cGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCg                                               >  8:169604/1‑149 (MQ=255)
                                                                                                                    aGGTACATCGGTCCGACGCCGGAACCGGTGGAAATAACACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGAGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGtt                                 >  8:132565/1‑149 (MQ=255)
                                                                                                                    aGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGtt                                 <  3:277879/149‑1 (MQ=255)
                                                                                                                    aGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGtt                                 >  6:404899/1‑149 (MQ=255)
                                                                                                                       tACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTtgtg                              <  6:62568/149‑1 (MQ=255)
                                                                                                                        tgaTCGGGCCGACGCCGGAACCGTTGGAAATACCAGCGGCAGTGTTGTTGGAGCTGGTGACGTTCGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCGCTGGGCGCCTTCGAACATGATGTCCTTGCCGGGGCGACGCAGGTtgtg                              <  4:336268/146‑1 (MQ=255)
                                                                                                                         cATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGACCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTtgtg                              >  7:275703/1‑147 (MQ=255)
                                                                                                                                 cGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTcggcg                    <  2:85683/149‑1 (MQ=255)
                                                                                                                                  gACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTcggcgg                   <  5:404899/149‑1 (MQ=255)
                                                                                                                                      ccactgccGGTGGTATTTCCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGAc               >  2:226420/7‑149 (MQ=255)
                                                                                                                                       cactgccGGTGGTATTTCCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACg              >  3:114983/6‑149 (MQ=255)
                                                                                                                                       cGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCg                                               >  4:289193/1‑116 (MQ=255)
                                                                                                                                       cGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCg                                               <  3:289193/116‑1 (MQ=255)
                                                                                                                                         ctgccGGTGGTATTTCCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTc            >  4:299016/4‑149 (MQ=255)
                                                                                                                                              ggTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGGAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGca       >  3:301164/1‑149 (MQ=255)
                                                                                                                                                   aaaTACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGCatcatc  <  7:169604/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGCATCATC  >  NC_002947/5558359‑5558653

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: