Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,540,091 C→T 6.2% G80S (GGC→AGC)  nudC ← NADH pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,540,0910CT 6.2% 187.8 / 2.8 66G80S (GGC→AGC) nudCNADH pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (2/2);  total (37/29)
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.73e-01

CCCGCTCCCAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTG  >  NC_002947/4539950‑4540236
                                                                                                                                             |                                                                                                                                                   
cccGCTCCCAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGc                                                                                                                                            >  1:952527/1‑150 (MQ=255)
      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)
                     tCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCa                                                                                                                                              >  1:1448085/1‑127 (MQ=255)
                     tCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCa                                                                                                                                              <  2:1448085/127‑1 (MQ=255)
                     tCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTgg                                                                                                                       >  1:1015835/1‑150 (MQ=255)
                          ggCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGc                                                                                                                                         >  1:651926/1‑127 (MQ=255)
                          ggCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGc                                                                                                                                         <  2:651926/127‑1 (MQ=255)
                           gCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTgc                                                                                                                 >  1:1503454/1‑150 (MQ=255)
                           gCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTgc                                                                                                                 >  2:157641/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                                                                                                                                    <  1:1623332‑M1/22‑1 (MQ=255)
                                             caggacctcgacgtgctgtgcgagcatggcatcggccactttgacggccagccggtgttcctgctggaactgcgcagcgccactgatgtgcCAGGCCGCGCCAGCTGCAGCCt                                                                                                                                    >  2:1623332‑M1/92‑113 (MQ=255)
                                                    aTGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTg                                                                                        <  2:421397/150‑1 (MQ=255)
                                                                  aGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGc                                                                         <  1:978239/150‑1 (MQ=255)
                                                                  aGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGc                                                                         <  2:365935/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                                                                                                                           <  2:737159/82‑1 (MQ=255)
                                                                                     aTAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCa                                                                                                                           >  1:737159/1‑82 (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                                                                                                <  1:1581860/97‑1 (MQ=255)
                                                                                                 cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAac                                                                                                >  2:1581860/1‑97 (MQ=255)
                                                                                                 cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATg                                                                    >  2:1353238/1‑124 (MQ=255)
                                                                                                 cTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATg                                                                    <  1:1353238/124‑1 (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                                                               >  2:23643/1‑111 (MQ=255)
                                                                                                                   gCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATg                                                               <  1:23643/111‑1 (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                                                                                                        <  2:1214725/61‑1 (MQ=255)
                                                                                                                             aTGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCca                                                                                                        >  1:1214725/1‑61 (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)
                                                                                                                                             |                                                                                                                                                   
CCCGCTCCCAGCGAATTTGCGTCATCGGCTGGCCGCAACTGCCGCAGAAGCGATGCTCGCGGGCCCAGGTGCCGATCTGTGCGGCATAGCCCAGCACCTTGTAGGTGTCGAAGTCGCCTTCGAGCATGAACGCACGCAGGCCGCGCCAGCTGCAGCCTGGCACATCAGTGGCGCTGCGCAGTTCCAGCAGGAACACCGGCTGGCCGTCAAAGTGGCCGATGCCATGCTCGCACAGCACGTCGAGGTCCTGGCGCTTTAGCCAGTCGCGGGGGAACAGCGCGCCGTTG  >  NC_002947/4539950‑4540236

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