Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,006,733 C→A 6.0% A381A (GCG→GCT)  PP_3535 ← gamma‑glutamyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,006,7330CA 6.0% 214.7 / 3.3 67A381A (GCG→GCT) PP_3535gamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base A (2/2);  total (33/34)
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.27e-01

TCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGG  >  NC_002947/4006585‑4006880
                                                                                                                                                    |                                                                                                                                                   
caacggCCAGGCCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGcgcc                                                                                                                                                     >  3:139373/5‑149 (MQ=255)
 cATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGcgccg                                                                                                                                                    >  8:155276/1‑149 (MQ=255)
     gCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCcgcgc                                                                                                                                                 <  2:189100/148‑1 (MQ=255)
         gggCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTg                                                                                                                                            >  2:9112/1‑149 (MQ=255)
           gCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCGGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGc                                                                                                                                          <  4:331070/149‑1 (MQ=255)
           gCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGc                                                                                                                                          <  5:86008/149‑1 (MQ=255)
            ccTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGcc                                                                                                                                         <  7:155276/149‑1 (MQ=255)
             cTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCa                                                                                                                                        <  7:289042/149‑1 (MQ=255)
             cTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGCCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCCATCAGCCCCGCGCGGATGGc                                                                                                                                          >  3:230520/1‑147 (MQ=255)
                  cGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTaa                                                                                                                                   <  4:139373/149‑1 (MQ=255)
                      gTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCa                                                                                                                                <  6:309702/148‑1 (MQ=255)
                      gTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCa                                                                                                                                <  5:154849/148‑1 (MQ=255)
                        cATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGaa                                                                                                                             >  7:327165/1‑149 (MQ=255)
                         aTCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAAt                                                                                                                            >  4:191947/1‑149 (MQ=255)
                                        gTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGa                                                                                                             >  2:236199/1‑149 (MQ=255)
                                         tGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGAc                                                                                                            <  3:191947/149‑1 (MQ=255)
                                           agagGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGAcc                                                                                                           <  1:23681/148‑1 (MQ=255)
                                                tGGAGATGCGCAGTGGGGTGCGGTCGGGCGCAAGCGCCATGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACg                                                                                                     <  4:230520/149‑1 (MQ=255)
                                                 ggAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAg                                                                                                                                  <  7:254706/119‑1 (MQ=255)
                                                 ggAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAg                                                                                                                                  >  8:254706/1‑119 (MQ=255)
                                                  gagaTGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGTGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACggg                                                                                                   >  2:19864/1‑149 (MQ=255)
                                                   agaTGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGc                                                                                                  >  4:62415/1‑149 (MQ=255)
                                                       gcgcAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGTCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCgg                                                                                                                                               >  8:126472/1‑100 (MQ=255)
                                                       gcgcAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCgg                                                                                                                                               <  7:126472/100‑1 (MQ=255)
                                                        cgcAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACg                                                                                               >  8:21490/1‑147 (MQ=255)
                                                        cgcAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTa                                                                                             <  1:9112/149‑1 (MQ=255)
                                                         gcAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAg                                                                                            >  7:135487/1‑149 (MQ=255)
                                                         gcAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAg                                                                                            >  5:204465/1‑149 (MQ=255)
                                                             tGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAAt                                                                                                                            >  8:244454/1‑113 (MQ=255)
                                                             tGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAAt                                                                                                                            <  7:244454/113‑1 (MQ=255)
                                                                cGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGt                                                                                     <  8:135487/149‑1 (MQ=255)
                                                                  aGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCg                                                                                   <  6:204465/149‑1 (MQ=255)
                                                                   gCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTc                                                                                    <  1:188827/147‑1 (MQ=255)
                                                                   gCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTc                                                                                    >  2:188827/1‑147 (MQ=255)
                                                                    cGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACgg                                                                                                    >  8:62309/1‑130 (MQ=255)
                                                                    cGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACgg                                                                                                    <  7:62309/130‑1 (MQ=255)
                                                                      acctggCGCTTGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGTCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGc                                                                                 >  4:264288/5‑147 (MQ=255)
                                                                       tCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCNGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAg                                                                                            <  1:321764/135‑1 (MQ=255)
                                                                       tCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGNGCNTTGACCGGGACGGGCACNTAg                                                                                            >  2:321764/1‑135 (MQ=255)
                                                                       tCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAg                                                                                            >  2:316675/1‑135 (MQ=255)
                                                                       tCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAg                                                                                            <  1:316675/135‑1 (MQ=255)
                                                                       tCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAg                                                                              >  3:292862/1‑149 (MQ=255)
                                                                         gggCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCg                                                                                         <  7:327600/136‑1 (MQ=255)
                                                                         gggCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCg                                                                                         >  8:327600/1‑136 (MQ=255)
                                                                           gcgcAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGtt                                                                                       <  7:69170/136‑1 (MQ=255)
                                                                           gcgcAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGtt                                                                                       >  8:69170/1‑136 (MQ=255)
                                                                             gcAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGaa                                                                                                                             >  1:327014/1‑96 (MQ=255)
                                                                             gcAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGaa                                                                                                                             <  2:327014/96‑1 (MQ=255)
                                                                               aaGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGAc                                                                                                            >  8:243297/1‑111 (MQ=255)
                                                                               aaGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGAc                                                                                                            <  7:243297/111‑1 (MQ=255)
                                                                                  cgcCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAgcgc                                                                   <  1:236199/149‑1 (MQ=255)
                                                                                   gcCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGc                                                                                                  >  7:10535/1‑117 (MQ=255)
                                                                                   gcCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGc                                                                                                  <  8:10535/117‑1 (MQ=255)
                                                                                                cTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAggg                                                                                                                   >  8:28515/1‑87 (MQ=255)
                                                                                                cTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAggg                                                                                                                   <  7:28515/87‑1 (MQ=255)
                                                                                                          gCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAg                                                        >  4:153350/1‑136 (MQ=255)
                                                                                                          gCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAg                                                        <  3:153350/136‑1 (MQ=255)
                                                                                                            ccGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGcc                                                                                                                 <  6:205512/77‑1 (MQ=255)
                                                                                                            ccGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGcc                                                                                                                 >  5:205512/1‑77 (MQ=255)
                                                                                                                     cgcgcTTCATGCTGTGCTCGCCGATCAGCGCAGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGAt                                 >  6:297518/1‑148 (MQ=255)
                                                                                                                     cgcgcTTCATGCTGTGCTCGCCGATCAGCGCAGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGAt                                 >  6:292627/1‑148 (MQ=255)
                                                                                                                      gcgcTTCCTGCTGGGCGCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAAGCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTGGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCa                               <  4:218310/149‑1 (MQ=255)
                                                                                                                         cTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGt                            <  8:327165/149‑1 (MQ=255)
                                                                                                                                cTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGc                     >  3:150386/1‑149 (MQ=255)
                                                                                                                                 tgtgCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGc                     >  1:330230/1‑148 (MQ=255)
                                                                                                                                        gCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGctcggcctcgg                                     <  8:88085/125‑1 (MQ=255)
                                                                                                                                         ccGATCAGCGCAGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAAcggcgg            <  5:297518/149‑1 (MQ=255)
                                                                                                                                         ccGATCAGCGCAGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAAcggcgg            <  5:292627/149‑1 (MQ=255)
                                                                                                                                           gATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCgggg          >  6:357598/1‑149 (MQ=255)
                                                                                                                                               aGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTc      >  5:271078/1‑149 (MQ=255)
                                                                                                                                                   ccgcGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAgg  <  2:330230/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGG  >  NC_002947/4006585‑4006880

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