Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,090,047 A→G 6.0% L297P (CTG→CCG)  astA‑II ← arginine N‑succinyltransferase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,090,0470AG 6.0% 197.4 / 2.3 68L297P (CTG→CCG) astA‑IIarginine N‑succinyltransferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base G (2/2);  total (33/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 = 3.36e-01

CGTCTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCT  >  NC_002947/5089901‑5090194
                                                                                                                                                   |                                                                                                                                                   
cGTCTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAgg                                                                                                                                                   <  3:252150/149‑1 (MQ=255)
  tatGTTCCCAGACCTCAAACCGCGACCACGCGCACGCTGGCGCCCACACCCACCCCCAGCGCGTCGGCTGCTGCCACGCTCACGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGACCAGAAACACCCCGCGCTACTCCTGCAAAATGcc                                                                                                                                                 >  2:81229/3‑149 (MQ=255)
   cTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGATGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCg                                                                                                                                                >  7:126842/1‑149 (MQ=255)
    tGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGt                                                                                                                                               >  5:164666/1‑149 (MQ=255)
      ttGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTg                                                                                                                                             <  5:147652/149‑1 (MQ=255)
        gCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACGGCCTGTTGc                                                                                                                                            <  2:278858/148‑1 (MQ=255)
             aCCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCa                                                                                                                                      <  7:25041/149‑1 (MQ=255)
               cTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAgg                                                                                                                                    <  4:301227/149‑1 (MQ=255)
                 cAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCTCCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTa                                                                                                                                  >  4:264451/1‑149 (MQ=255)
                 cAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTa                                                                                                                                  >  8:112042/1‑149 (MQ=255)
                  aaaCCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTAc                                                                                                                                 <  1:228997/149‑1 (MQ=255)
                       gcgACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGgcg                                                                                                                            >  3:256553/1‑149 (MQ=255)
                        cgACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGgcgc                                                                                                                           <  8:279569/149‑1 (MQ=255)
                         gACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGcgcc                                                                                                                          <  8:126842/149‑1 (MQ=255)
                          aCCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGcgccc                                                                                                                         <  4:92026/149‑1 (MQ=255)
                              ggCGCACGCTGGCGCCCTCCCCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCt                                                                                                                     >  3:190189/1‑149 (MQ=255)
                              ggCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCt                                                                                                                     <  4:218888/149‑1 (MQ=255)
                              ggCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCt                                                                                                                     >  3:193989/1‑149 (MQ=255)
                                cgcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAGCACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCttt                                                                                                                   <  3:47324/149‑1 (MQ=255)
                                cgcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCttt                                                                                                                   <  3:121792/149‑1 (MQ=255)
                                 gcACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTg                                                                                                                  <  8:31277/149‑1 (MQ=255)
                                         gcgcCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGc                                                                                                          <  8:141401/149‑1 (MQ=255)
                                               tCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGcaca                                                                                                                                         <  7:190749/112‑1 (MQ=255)
                                               tCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGcaca                                                                                                                                         >  8:190749/1‑112 (MQ=255)
                                                  cccTCGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACGGGGGCCTCTTTGAt                                                                                                 <  1:81229/149‑1 (MQ=255)
                                                   ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACCGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                >  2:61394/1‑149 (MQ=255)
                                                   ccACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATc                                                                                                >  4:188978/1‑149 (MQ=255)
                                                      cGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTc                                                                                             <  1:182924/149‑1 (MQ=255)
                                                        cccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCAc                                                                                           <  4:190189/149‑1 (MQ=255)
                                                        cccAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCAc                                                                                           <  4:193989/149‑1 (MQ=255)
                                                           aGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCgg                                                                                        >  1:303369/1‑149 (MQ=255)
                                                                           ccACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTgg                                                                        <  6:164666/149‑1 (MQ=255)
                                                                                   aGGCTGACCGGCTTGCCAGGCTCCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATc                                                                  <  3:264451/147‑1 (MQ=255)
                                                                                     gCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGc                                                              >  2:19970/1‑149 (MQ=255)
                                                                                       tGACCGGCTTGCCAGGCACCCAGTCCAGGGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACCCCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCg                                                             >  1:170817/1‑149 (MQ=255)
                                                                                       tGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCg                                                             >  8:106714/1‑148 (MQ=255)
                                                                                        gACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGa                                                           >  2:52134/1‑149 (MQ=255)
                                                                                          ccGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGAt                                                          <  8:28526/148‑1 (MQ=255)
                                                                                           cGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGa                                                           >  5:308314/1‑146 (MQ=255)
                                                                                           cGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGa                                                           <  6:308314/146‑1 (MQ=255)
                                                                                           cGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGc                                                        >  5:136990/1‑149 (MQ=255)
                                                                                           cGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGCGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGc                                                        <  1:80683/149‑1 (MQ=255)
                                                                                           cGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGAAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGc                                                        >  2:80683/1‑149 (MQ=255)
                                                                                             gCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAAGACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCa                                                      <  8:127285/149‑1 (MQ=255)
                                                                                                 gcgaGGCTCCCAGTGCA‑GTTCCAGCAACACGGCGCGGTAGTCCTGCAACAGGCCGTTGTAGACCAGGTACGGGCGCCCGCCTTTGACCggg                                                                                                            <  4:105441/88‑1 (MQ=255)
                                                                                                 gCCAGGCACCCAGTCCA‑GGTGCAGAAACACCGCGCGGTAGTCCTGCAACAGGCCGTTCTACACCACGTACGGGCGCCCCCCATAGACCgcg                                                                                                            >  3:105441/1‑89 (MQ=255)
                                                                                                 gCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCagag                                                   <  1:52134/148‑1 (MQ=255)
                                                                                                   cAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTg                                                >  2:124884/1‑149 (MQ=255)
                                                                                                        aCCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCAGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTgcgcgc                                           <  7:112042/149‑1 (MQ=255)
                                                                                                         cccAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCa                                          >  5:140783/1‑149 (MQ=255)
                                                                                                            aGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATg                                        >  6:69589/1‑148 (MQ=255)
                                                                                                                 a‑ggtccagCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGt                                  >  6:255662/1‑149 (MQ=255)
                                                                                                                       cagCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGgc                             >  1:115774/1‑149 (MQ=255)
                                                                                                                               ccGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCg                     <  6:136990/149‑1 (MQ=255)
                                                                                                                                cgcgcgGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCg                     >  5:162783/1‑148 (MQ=255)
                                                                                                                                 gcgcgGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAg                                     <  4:195694/131‑1 (MQ=255)
                                                                                                                                 gcgcgGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAg                                     >  3:195694/1‑131 (MQ=255)
                                                                                                                                 gcgcgGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGaa                   <  4:256553/149‑1 (MQ=255)
                                                                                                                                         gTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCg            <  5:69589/148‑1 (MQ=255)
                                                                                                                                          gtgtGCAACGGCCTGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGAt          >  4:82341/4‑149 (MQ=255)
                                                                                                                                           tgtGCAACGGCCTGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTc                      >  7:288029/3‑136 (MQ=255)
                                                                                                                                           tgtGCAACGGCCTGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTc                      <  8:288029/134‑1 (MQ=255)
                                                                                                                                             tGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg         >  6:145241/1‑147 (MQ=255)
                                                                                                                                             tGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATg         >  5:216853/1‑147 (MQ=255)
                                                                                                                                             tGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTa       >  8:23772/1‑149 (MQ=255)
                                                                                                                                               cAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTa       <  6:216853/147‑1 (MQ=255)
                                                                                                                                               cAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGt     >  8:288437/1‑149 (MQ=255)
                                                                                                                                                  cAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCt  <  2:303369/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGTCTGTTGCCAGACCTCAAACCGCGACCAGGCGCACGCTGGCGCCCTCACCCACGCCCAGCGCGTCGGCTGCTGCCACGCTCAGGCTGACCGGCTTGCCAGGCACCCAGTCCA‑GGTCCAGCAACACCGCGCGGTAGTCCTGCAACAGGCCGTTGCACACCAGGTACGGGCGCCCGCCTTTGACCGGGGCCTCTTCGATCTTCACCGGCACCACCCGGCTCTGGGCGATCGAGCGGATGCCAGAGGTGCGCGCATGCAGGGTCGGGCCGCCGTCGAAGATGTCGATGTAGTGCT  >  NC_002947/5089901‑5090194

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: ‑