Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,089,078:1 +C 5.8% coding (761/1029 nt) astA‑I ← arginine N‑succinyltransferase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,089,0781.C 5.8% 276.7 / 10.7 69coding (761/1029 nt)astA‑Iarginine N‑succinyltransferase subunit beta
Reads supporting (aligned to +/- strand):  ref base . (32/33);  new base C (2/2);  total (34/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.06e-01

TGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAACTTGGGCATCAG  >  NC_002947/5088930‑5089226
                                                                                                                                                     |                                                                                                                                                    
tGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGGCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGgccgggccgg                                                                                                                                                       <  4:75364/149‑1 (MQ=255)
   ggCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCAt                                                                                                                                                   >  1:129119/1‑149 (MQ=255)
      agggATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGCCGGCCCGGCCATCGa                                                                                                                                                <  6:152039/146‑1 (MQ=255)
       gTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGaa                                                                                                                                               <  2:52593/149‑1 (MQ=255)
       gTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGCCGGCCCGGCCATCGag                                                                                                                                               >  5:283072/1‑148 (MQ=255)
       gTGATGCGACAGTCGTCGCGCTTGCGGTTGGGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGaa                                                                                                                                               <  7:9766/149‑1 (MQ=255)
        tGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCg                                                                                                                                                 >  6:300360/1‑146 (MQ=255)
        tGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAg                                                                                                                                              <  6:266490/149‑1 (MQ=255)
          aTGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCTCTGGCCGGGCCGGCATCGAAGAt                                                                                                                                            <  7:47447/149‑1 (MQ=255)
          aTGCGACAGTCGGCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGa                                                                                                                                             <  1:285455/148‑1 (MQ=255)
             cGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTc                                                                                                                                         <  4:269843/149‑1 (MQ=255)
             cGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTc                                                                                                                                         >  5:172697/1‑149 (MQ=255)
             cGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTc                                                                                                                                         <  6:172697/149‑1 (MQ=255)
             cGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGCCGGCCCGGCCATCGagtgtgac                                                                                                                                         <  6:283072/149‑8 (MQ=255)
               aCAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGa                                                                                                                                       >  4:221446/1‑149 (MQ=255)
               aCAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGa                                                                                                                                       <  2:64967/149‑1 (MQ=255)
               aCAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGa                                                                                                                                       <  5:300360/149‑1 (MQ=255)
                   tcgtcgCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACg                                                                                                                                     >  6:287651/1‑147 (MQ=255)
                    cgtcgCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAg                                                                                                                                  >  3:171145/1‑149 (MQ=255)
                      tcgCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGt                                                                                                                                    >  5:189627/1‑145 (MQ=255)
                      tcgCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGt                                                                                                                                    <  6:189627/145‑1 (MQ=255)
                              gCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTa                                                                                                                          <  3:271579/147‑1 (MQ=255)
                               cGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATTGAAGATGTCGACGTAGCCCTGGTAGtt                                                                                                                       <  1:182561/149‑1 (MQ=255)
                                ggTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGtt                                                                                                                       >  3:48247/1‑148 (MQ=255)
                                  ttGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGaa                                                                                                                    <  4:171145/149‑1 (MQ=255)
                                      ggATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCt                                                                                                                >  2:229228/1‑149 (MQ=255)
                                             gaAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTc                                                                                                              <  5:89918/144‑1 (MQ=255)
                                             gaAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTc                                                                                                              >  6:89918/1‑144 (MQ=255)
                                             gaAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCt                                                                                                         >  3:247395/1‑149 (MQ=255)
                                                aGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGa                                                                                                      <  3:221446/149‑1 (MQ=255)
                                                    gTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCAt                                                                                                  <  4:247395/149‑1 (MQ=255)
                                                       gcgTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGc                                                                                               >  5:150686/1‑149 (MQ=255)
                                                            aTCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAgcg                                                                                          >  4:26223/1‑149 (MQ=255)
                                                             tCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCt                                                                                                         >  4:48254/1‑133 (MQ=255)
                                                                  tGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCCGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCt                                                                                    >  5:72074/1‑149 (MQ=255)
                                                                      gTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGgt                                                                                <  2:63415/149‑1 (MQ=255)
                                                                          ctaccGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCg                                                                                                                                        >  2:123799/1‑89 (MQ=255)
                                                                          ctaccGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCg                                                                                                                                        <  1:123799/89‑1 (MQ=255)
                                                                            accGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCttt                                                                          >  4:210342/1‑149 (MQ=255)
                                                                             ccGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCttt                                                                          <  4:48247/148‑1 (MQ=255)
                                                                                     aCCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTtgt                                                                        <  4:17433/142‑1 (MQ=255)
                                                                                     aCCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTtgt                                                                        >  3:17433/1‑142 (MQ=255)
                                                                                      ccAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACgg                                                                >  7:36735/1‑149 (MQ=255)
                                                                                      ccAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACgg                                                                >  7:41571/1‑149 (MQ=255)
                                                                                          ggTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGa                                                            <  7:33026/149‑1 (MQ=255)
                                                                                                    tcGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTgcg                                                  >  8:277038/1‑149 (MQ=255)
                                                                                                    tcGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTgcg                                                  >  1:192036/1‑149 (MQ=255)
                                                                                                       cgTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTgcgcgc                                               >  3:281144/1‑149 (MQ=255)
                                                                                                              cgcgGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCg                                        >  1:145351/1‑149 (MQ=255)
                                                                                                              cgcgGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCg                                        <  3:210342/149‑1 (MQ=255)
                                                                                                               gcgGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCg                                        >  8:122243/1‑148 (MQ=255)
                                                                                                               gcgGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCg                                        >  6:299990/1‑148 (MQ=255)
                                                                                                                cgGATCTTGGCAGTGTCACACTCGATGGCCGGGCCAGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGa                                      <  6:72074/149‑1 (MQ=255)
                                                                                                                  gATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACa                                    <  8:36735/149‑1 (MQ=255)
                                                                                                                  gATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACa                                    <  8:41571/149‑1 (MQ=255)
                                                                                                                        ggCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAg                               >  2:153527/1‑148 (MQ=255)
                                                                                                                        ggCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAg                               <  1:229228/148‑1 (MQ=255)
                                                                                                                         gCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGc                              <  4:281144/148‑1 (MQ=255)
                                                                                                                          cAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAg                            <  2:129119/149‑1 (MQ=255)
                                                                                                                            gtgtCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGgt                          <  7:122243/149‑1 (MQ=255)
                                                                                                                                 acacTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATg                                                                                                 <  5:46667/73‑1 (MQ=255)
                                                                                                                                 acacTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATg                                                                                                 >  6:46667/1‑73 (MQ=255)
                                                                                                                                 acacTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAg                               >  6:91932/1‑139 (MQ=255)
                                                                                                                                 acacTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAg                               <  5:91932/139‑1 (MQ=255)
                                                                                                                                 acacTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTa                       >  1:23806/1‑147 (MQ=255)
                                                                                                                                 acacTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTa                       >  6:264694/1‑147 (MQ=255)
                                                                                                                                    ttCGATGCCGGCCCGGCCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACATCg                  >  8:245046/2‑149 (MQ=255)
                                                                                                                                       gATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCgg                 <  2:192036/147‑1 (MQ=255)
                                                                                                                                         tGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAAc             <  3:26223/149‑1 (MQ=255)
                                                                                                                                         tGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAAc             >  4:221705/1‑149 (MQ=255)
                                                                                                                                           gccgggccggCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAACtt           <  7:277038/149‑1 (MQ=255)
                                                                                                                                            ccgggccggCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAACTTg          <  5:299990/149‑1 (MQ=255)
                                                                                                                                             cgggccggCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAACTTgg         <  2:23806/149‑1 (MQ=255)
                                                                                                                                              gggccggCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCg                                                                                  >  8:115989/1‑75 (MQ=255)
                                                                                                                                              gggccggCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCg                                                                                  <  7:115989/75‑1 (MQ=255)
                                                                                                                                                    gCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAGGCAGGTGTACAGCGGGAACTTGGGCATCAg  >  6:267347/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTACGGCGCGGATCTTGGCAGTGTCACACTCGATGGCCGGGCCGGCATCGAAGATGTCGACGTAGCCCTGGTAGTTGAAGCCTTCCTGCTTGAGCATGGCCAGCGCCGGCTCGGTGTCTTTGTGCACACGGCCGATCACACTGCGCGCCGCTTCGGACAGGAAGCAGGTGTACAGCGGGAACTTGGGCATCAG  >  NC_002947/5088930‑5089226

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: