Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,045,722 C→A 6.9% R130L (CGT→CTT) ‡ PP_5298 ← glutamine amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,045,7220CA 6.9% 221.1 / 5.7 73R130L (CGT→CTT) ‡PP_5298glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base A (3/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

TGACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACC  >  NC_002947/6045579‑6045861
                                                                                                                                               |                                                                                                                                           
tGACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCCGTGGTCAAGGTAGCCAGGCAGTTCCTGCCCACGCTgg                                                                                                                                        >  7:159396/1‑149 (MQ=255)
 gACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTgg                                                                                                                                        <  2:361591/148‑1 (MQ=255)
 gACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTgg                                                                                                                                        <  4:178117/148‑1 (MQ=255)
    cTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTg                                                                                                                                      >  4:128752/1‑147 (MQ=255)
     tCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAg                                                                                                                                   >  3:341417/1‑149 (MQ=255)
       gAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAg                                                                                                                                   <  1:150069/147‑1 (MQ=255)
       gAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAg                                                                                                                                   <  1:211706/147‑1 (MQ=255)
       gAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAg                                                                                                                                   <  8:149298/147‑1 (MQ=255)
        aaCTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACCAGCGTGTGCAggaa                                                                                                                                <  3:282255/149‑3 (MQ=255)
                cacgaACCCCAGGGGCTGGAACAACCCGCCAGGGTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAATGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCGaa                                                                                                                        <  7:62244/149‑1 (MQ=255)
                 accaACCCCAGGCGCTCGACCAAGCCGCCATGCTCGTTGGTGACCTCGTGGCGAGGGCCTTATTGCACCTCCAGGGGGGCGTCCTAAGGTTCGCGGTGTGCAAGGTAGCCAGGCTGTTCCCGCACACGCTGGTGCAGGCTgcc                                                                                                                             <  7:61829/143‑1 (MQ=255)
                     aCCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTTCACACGCTGGTGCAGGCTGCCGCCAAGGGcc                                                                                                                   >  4:69408/1‑149 (MQ=255)
                         cAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGt                                                                                                               >  3:350808/1‑149 (MQ=255)
                           ggCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGtt                                                                                                              >  2:128511/1‑148 (MQ=255)
                            gcgcTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAgg                                                                                                                      <  8:136607/139‑1 (MQ=255)
                            gcgcTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAgg                                                                                                                      >  7:136607/1‑139 (MQ=255)
                                  gAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCt                                                                                                      >  5:392346/1‑149 (MQ=255)
                                   aacaacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTg                                                                                                     >  7:298226/1‑149 (MQ=255)
                                   aacaacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTg                                                                                                     >  6:65587/1‑149 (MQ=255)
                                      aacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGt                                                                                                                                       <  4:376127/112‑1 (MQ=255)
                                      aacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGt                                                                                                                                       >  3:376127/1‑112 (MQ=255)
                                      aacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTgg                                                                                                    >  2:382004/1‑147 (MQ=255)
                                      aacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTa                                                                                                  <  7:180381/149‑1 (MQ=255)
                                        cccGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACCAGCGTGTGCAgg                                                                                                                                  >  7:262319/1‑115 (MQ=255)
                                        cccGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACCAGCGTGTGCAgg                                                                                                                                  <  8:262319/115‑1 (MQ=255)
                                        cccGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGc                                                                                                >  7:101511/1‑149 (MQ=255)
                                          cGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCt                                                                                              <  5:379465/149‑1 (MQ=255)
                                            cccgggtCGATGCTGACCGTGTGGCGAGGGCCGTATTGCACCTCCTGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCGGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCg                                                                                            <  6:392346/143‑1 (MQ=255)
                                             cAGGCTCGATTCTGACCGTGTGGCGATGGCCGTATTGCTCCTCCAGGGGGGCGTCGTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCgg                                                                                           <  8:159396/149‑1 (MQ=255)
                                              aGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGc                                                                                          <  5:65587/149‑1 (MQ=255)
                                                  tCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGAt                                                                                      >  7:133076/1‑149 (MQ=255)
                                                     aTGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCt                                                                                                         >  4:147798/1‑127 (MQ=255)
                                                     aTGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCt                                                                                                         <  3:147798/127‑1 (MQ=255)
                                                     aTGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCa                                                                                   <  8:133076/149‑1 (MQ=255)
                                                           aCCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATAc                                                                             <  1:128511/149‑1 (MQ=255)
                                                           aCCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATAc                                                                             >  3:125413/1‑149 (MQ=255)
                                                              gAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACggg                                                                          <  1:382004/149‑1 (MQ=255)
                                                              gAGCGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTACACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACggg                                                                          <  7:57228/149‑1 (MQ=255)
                                                                   gCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATTCAGAATACGGGGacac                                                                     >  5:105889/1‑149 (MQ=255)
                                                                     gAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCg                                                                   >  3:352611/1‑149 (MQ=255)
                                                                     gAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCAGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACAcc                                                                    >  1:129211/1‑148 (MQ=255)
                                                                               ttGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATg                                                          <  4:350808/148‑1 (MQ=255)
                                                                                  cACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGct                                                       <  3:133773/148‑1 (MQ=255)
                                                                                      tCCCGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgcct                                                   >  5:387792/1‑148 (MQ=255)
                                                                                      tCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgcctg                                                  >  2:37825/1‑149 (MQ=255)
                                                                                      tCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgcctg                                                  >  2:353169/1‑149 (MQ=255)
                                                                                      tCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgc                                                     >  5:417314/1‑146 (MQ=255)
                                                                                      tCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgc                                                     <  6:417314/146‑1 (MQ=255)
                                                                                      tCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCAGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgcctg                                                  >  8:349104/1‑149 (MQ=255)
                                                                                        ccggggggCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgcct                                                   <  6:387792/144‑1 (MQ=255)
                                                                                        cAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCa                                                <  8:101511/149‑1 (MQ=255)
                                                                                         aGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGcc                                                                                                                   <  3:1513/81‑1 (MQ=255)
                                                                                         aGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGcc                                                                                                                   >  4:1513/1‑81 (MQ=255)
                                                                                         aGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCaa                                               <  2:278555/149‑1 (MQ=255)
                                                                                                 gTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGAACACCGACGGCAATGGCTGCCTGCAACAGCgg                                         >  7:68692/1‑147 (MQ=255)
                                                                                                  tCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCg                                          <  5:185025/145‑1 (MQ=255)
                                                                                                  tCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCg                                          >  6:185025/1‑145 (MQ=255)
                                                                                                    cTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAggg                                    >  3:379717/1‑149 (MQ=255)
                                                                                                    cTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAggg                                    >  3:162893/1‑149 (MQ=255)
                                                                                                    cTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAgg                                     <  4:341417/148‑1 (MQ=255)
                                                                                                     tCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGc                                           <  4:86066/141‑1 (MQ=255)
                                                                                                     tCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGc                                           >  3:86066/1‑141 (MQ=255)
                                                                                                       aGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCa                                 <  3:128752/149‑1 (MQ=255)
                                                                                                              cgGTGGTCAAGGGAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATc                          <  3:69408/149‑1 (MQ=255)
                                                                                                                gtggtCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCAc                        <  4:352611/149‑1 (MQ=255)
                                                                                                                     cAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACgcgc                    >  8:163480/1‑148 (MQ=255)
                                                                                                                          tAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATc              <  4:125413/149‑1 (MQ=255)
                                                                                                                            gCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCAt            <  1:353169/149‑1 (MQ=255)
                                                                                                                               aGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCAcgc                      <  2:154459/136‑1 (MQ=255)
                                                                                                                               aGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCAcgc                      >  1:154459/1‑136 (MQ=255)
                                                                                                                                 gCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCgg       <  7:349104/149‑1 (MQ=255)
                                                                                                                                      gccTGCACCAGCGTGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTAcc  >  7:21729/2‑149 (MQ=255)
                                                                                                                                      gccTGCACCAGCGTGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTAcc  >  4:377067/2‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                           
TGACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACC  >  NC_002947/6045579‑6045861

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: