Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,509 A→T 10.8% G248G (GGT→GGA purA ← adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,558,5090AT10.8% 178.7 / 11.5 65G248G (GGT→GGApurAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base T (3/4);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

GTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGCATCATC  >  NC_002947/5558365‑5558653
                                                                                                                                                 |                                                                                                                                                
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                                                                                            <  2:175659/149‑1 (MQ=255)
                                                   ggAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATg                                                                                            <  5:14757/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)
                                                                                                                                                 |                                                                                                                                                
GTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGCATCATC  >  NC_002947/5558365‑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: