Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,045,725 T→G 7.9% Q129P (CAG→CCG)  PP_5298 ← glutamine amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,045,7250TG 7.9% 199.1 / 7.7 75Q129P (CAG→CCG) PP_5298glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base G (3/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

TGACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACCCGCCACACTGGG  >  NC_002947/6045579‑6045873
                                                                                                                                                  |                                                                                                                                                    
tGACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCCGTGGTCAAGGTAGCCAGGCAGTTCCTGCCCACGCTgg                                                                                                                                                    >  7:159396/1‑149 (MQ=255)
 gACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTgg                                                                                                                                                    <  4:178117/148‑1 (MQ=255)
 gACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTgg                                                                                                                                                    <  2:361591/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                                                                                                                 >  6:65587/1‑149 (MQ=255)
                                   aacaacCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTg                                                                                                                 >  7:298226/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                                                                                                                                              <  8:262319/115‑1 (MQ=255)
                                        cccGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACCAGCGTGTGCAgg                                                                                                                                              >  7:262319/1‑115 (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                                                                                                                     <  3:147798/127‑1 (MQ=255)
                                                     aTGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCt                                                                                                                     >  4:147798/1‑127 (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                                                                 <  6:417314/146‑1 (MQ=255)
                                                                                      tCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGctgc                                                                 >  5:417314/1‑146 (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                                                      >  6:185025/1‑145 (MQ=255)
                                                                                                  tCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCg                                                      <  5:185025/145‑1 (MQ=255)
                                                                                                    cTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAggg                                                >  3:379717/1‑149 (MQ=255)
                                                                                                    cTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAggg                                                >  3:162893/1‑149 (MQ=255)
                                                                                                    cTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAgg                                                 <  4:341417/148‑1 (MQ=255)
                                                                                                     tCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGc                                                       >  3:86066/1‑141 (MQ=255)
                                                                                                     tCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGc                                                       <  4:86066/141‑1 (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)
                                                                                                                                                gCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACCCGCCACACTg    >  1:176242/1‑149 (MQ=255)
                                                                                                                                                 cTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACCCGCCACACTgg   <  2:238845/149‑1 (MQ=255)
                                                                                                                                                  tGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACCCGCCACACTggg  >  3:349390/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGACCTCGAACTGCGCCACCAACCCCAGGCGCTCGAACAACCCGCCAGGCTCGATGCTGACCGAGTGGCGAGGGCCGTATTGCACCTCCAGGGGGGCGTCCTCAGGTTCGCGGTGGTCAAGGTAGCCAGGCAGTTCCTGCACACGCTGGTGCAGGCTGCCGCCAAGGGCCACGTTCAGCTCCTGGTAGCCTCGGCAGATGCAGAATACGGGCACACCGACGGCAATGGCTGCCTGCAACAGCGGCAGGGTCAGGCGATCACGCGCCAGATCATGCCGGGTACCCGCCACACTGGG  >  NC_002947/6045579‑6045873

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