Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,540,103 A→G 6.0% C76R (TGC→CGC)  nudC ← NADH pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,540,1030AG 6.0% 202.9 / 2.4 67C76R (TGC→CGC) nudCNADH pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base G (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CCCAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTGGCGTCCACCAGAAA  >  NC_002947/4539956‑4540250
                                                                                                                                                    |                                                                                                                                                    
cccAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGc                                                                                                                                                    <  2:405581/150‑1 (MQ=255)
  cAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCCCCCCCCCCCCCCCCCCCCCCCCGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCt                                                                                                                                                  >  2:133296/1‑150 (MQ=255)
    gCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCATGCAGGTGGTGGGAGGTGCAGCCTgg                                                                                                                                                <  2:677472/150‑1 (MQ=255)
         tttGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACAt                                                                                                                                           >  1:421397/1‑150 (MQ=255)
         tttGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCAACGTGCAGATCTTTTCACCGCATCCCAGCACCTTGTAGTTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACAt                                                                                                                                           <  1:778904/150‑1 (MQ=255)
               tCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTgg                                                                                                                                     >  1:1015835/1‑150 (MQ=255)
                     gCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTgc                                                                                                                               >  2:157641/1‑150 (MQ=255)
                     gCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTgc                                                                                                                               >  1:1503454/1‑150 (MQ=255)
                         gCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAg                                                                                                                           <  1:157641/150‑1 (MQ=255)
                          ccGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCACCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGt                                                                                                                          >  1:578266/1‑150 (MQ=255)
                              aaCTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCca                                                                                                                      >  1:355305/1‑150 (MQ=255)
                                 tGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCcagca                                                                                                                   >  2:174507/1‑150 (MQ=255)
                                  gCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCcagcag                                                                                                                  <  2:1503454/150‑1 (MQ=255)
                                      gaGAAGCAATTCTCGCAGGTCCAGGAGCCGATCCGTGCGGGATACCCCAGCACCTTGCAGGTGTCGAAGTCGCCTGCGAGCAAGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAac                                                                                                              <  2:91796/149‑1 (MQ=255)
                                      cAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGTCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAac                                                                                                              >  1:1184045/1‑150 (MQ=255)
                                       caggacctcgacgtgctgtgcgagcatggcatcggccactttgacggccagccggtgttcctgctggaactgcgcagcgccactgatgtgcCAGGCCGCGCCAGCTGCAGCCt                                                                                                                                                  >  2:1623332‑M1/92‑113 (MQ=255)
                                       caggacctcgacgtgctgtgcgagcatggcatcggccactttgacggccagccggtgttcctgctggaactgcgcagcgccactgatgtgcCAGGCCGCGCCAGCTGCAGCCt                                                                                                                                                  <  1:1623332‑M1/22‑1 (MQ=255)
                                              aTGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTg                                                                                                      <  2:421397/150‑1 (MQ=255)
                                                            aGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGc                                                                                       <  2:365935/150‑1 (MQ=255)
                                                            aGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGc                                                                                       <  1:978239/150‑1 (MQ=255)
                                                             ggTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGcc                                                                                      >  1:1560314/1‑150 (MQ=255)
                                                                 ccGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATg                                                                                  >  1:152212/1‑150 (MQ=255)
                                                                        gtgCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGACCGCACGCAGGCCGCGCCACCTGCAGTCTGGCACATCAGTGGAGCTGCGCAGGTACTGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCt                                                                           >  1:108803/1‑150 (MQ=255)
                                                                           cGGCATAGCCCAGCAGCTTGCAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCCAGCCGTCAAAGGTGGCCATGCCATGATCGc                                                                        >  2:192294/1‑150 (MQ=255)
                                                                               aTAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCa                                                                                                                                         >  1:737159/1‑82 (MQ=255)
                                                                               aTAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCa                                                                                                                                         <  2:737159/82‑1 (MQ=255)
                                                                                   cccAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACg                                                                >  2:332769/1‑150 (MQ=255)
                                                                                    ccAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGt                                                               >  2:1589964/1‑150 (MQ=255)
                                                                                     cAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGGGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAg                                                                                           >  1:616115/1‑121 (MQ=255)
                                                                                     cAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGcc                                                                                                                                                   >  2:1036204/1‑66 (MQ=255)
                                                                                     cAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGcc                                                                                                                                                   <  1:1036204/66‑1 (MQ=255)
                                                                                     cAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAg                                                                                           <  2:616115/121‑1 (MQ=255)
                                                                                        cACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAg                                                           <  2:1500150/150‑1 (MQ=255)
                                                                                        cACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTACAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGAACAGCACGTCGAg                                                           <  1:1587102/150‑1 (MQ=255)
                                                                                           cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAac                                                                                                              >  2:1581860/1‑97 (MQ=255)
                                                                                           cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAac                                                                                                              <  1:1581860/97‑1 (MQ=255)
                                                                                           cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATg                                                                                  <  1:1353238/124‑1 (MQ=255)
                                                                                           cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATg                                                                                  >  2:1353238/1‑124 (MQ=255)
                                                                                             tGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAacac                                                                                                            <  1:1561858/97‑1 (MQ=255)
                                                                                             tGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAacac                                                                                                            >  2:1561858/1‑97 (MQ=255)
                                                                                                aGGTGTCGAAGTCGCCTTCGAGCATCAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGgc                                                   <  1:354658/150‑1 (MQ=255)
                                                                                                      cGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTa                                             <  2:1560314/150‑1 (MQ=255)
                                                                                                       gAAGTCGCCTTCGAGCATGAAGGCACGCATGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAg                                            >  1:1111407/1‑150 (MQ=255)
                                                                                                       gAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAg                                            >  2:929781/1‑150 (MQ=255)
                                                                                                          gTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGc                                                                                 <  2:1528237/110‑1 (MQ=255)
                                                                                                          gTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGc                                                                                 >  1:1528237/1‑110 (MQ=255)
                                                                                                             gCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATg                                                                             <  1:23643/111‑1 (MQ=255)
                                                                                                             gCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATg                                                                             >  2:23643/1‑111 (MQ=255)
                                                                                                               cTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGcaca                                                                     >  1:678703/1‑117 (MQ=255)
                                                                                                               cTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGcaca                                                                     <  2:678703/117‑1 (MQ=255)
                                                                                                                ttCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTcgcg                                   >  1:865697/1‑150 (MQ=255)
                                                                                                                     gccactgatgtgcCAGGCTGCAGCTGGCGCGGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTcc                                                       <  2:498833‑M1/112‑1 (MQ=255)
                                                                                                                     gccactgatgtgcCAGGCTGCAGCTGGCGCGGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTcc                                                       >  1:498833‑M1/14‑125 (MQ=255)
                                                                                                                     gCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGCGGCGCTGCGCAGTTCCAGCAGGAACACCGTCTGGCCGTCAAATTGGCCGATGCCATGCTCGCACAGCACGTCGCGGTCCTGGCGCTTTAGCCAGTCGCGGGTGa                              >  2:1569735/1‑150 (MQ=255)
                                                                                                                       aTGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCca                                                                                                                      >  1:1214725/1‑61 (MQ=255)
                                                                                                                       aTGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCca                                                                                                                      <  2:1214725/61‑1 (MQ=255)
                                                                                                                        actgatgtgcCAGGCTGCAGCTGGCGCGGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGca                                                                       <  1:454301‑M1/96‑1 (MQ=255)
                                                                                                                        actgatgtgcCAGGCTGCAGCTGGCGCGGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGca                                                                       >  2:454301‑M1/11‑106 (MQ=255)
                                                                                                                           acgcacgcAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAgcg                        >  1:1034101/1‑150 (MQ=255)
                                                                                                                              cacgcaGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAgcgcgc                     <  1:1589964/150‑1 (MQ=255)
                                                                                                                                   AGGCCGCGTCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGGTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGACATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTg                >  2:1184045‑M1/1‑150 (MQ=255)
                                                                                                                                        gcgcCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCTGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTGGCGTc           <  2:952527/150‑1 (MQ=255)
                                                                                                                                           gcCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTGGCGTcca         >  2:1086649/1‑150 (MQ=255)
                                                                                                                                            ccAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTGGCGTccac        <  2:578266/150‑1 (MQ=255)
                                                                                                                                                  gCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGACGTCAAGGAGGCAGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTGGCGTCCACCAGaaa  <  1:1605237/150‑1 (MQ=255)
                                                                                                                                                   cAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCa                                         >  2:477333/1‑110 (MQ=255)
                                                                                                                                                   cAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCa                                         <  1:477333/110‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTGGCGTCCACCAGAAA  >  NC_002947/4539956‑4540250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: