Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,673 A→G 7.0% I292T (ATC→ACC)  hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,6730AG 7.0% 201.7 / 4.2 71I292T (ATC→ACC) hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

CACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCG  >  NC_002947/5559524‑5559820
                                                                                                                                                     |                                                                                                                                                   
cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGACAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcga                                                                                                                                                     <  2:987201/150‑1 (MQ=255)
 aCTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat                                                                                                                                                    <  1:683334/150‑1 (MQ=255)
 aCTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat                                                                                                                                                    <  2:221718/150‑1 (MQ=255)
    tggcGGCATCTGGATGCCGCCAGTTGGCAAGACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgatcga                                                                                                                                                 >  2:1013197/4‑150 (MQ=255)
    gTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgatcga                                                                                                                                                 <  1:621707/150‑1 (MQ=255)
     tCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgatcgat                                                                                                                                                >  2:35328/1‑150 (MQ=255)
        ggCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTgg                                                                                                                                             >  2:492445/1‑150 (MQ=255)
        cgCATCCAGATGCCGCCAGTTGGGCATACGACCCCCCCTCCCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCCGGGGCCGGGCGCGCCCCGCCAAAGTCGGCCCCGATATCGGCGTAGCGACCGCCCTGGGCGATCGATTgg                                                                                                                                             <  1:1082498/149‑1 (MQ=255)
         gCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGc                                                                                                                                            >  1:851209/1‑150 (MQ=255)
             ccAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTc                                                                                                                                        <  1:492445/150‑1 (MQ=255)
                     cgccAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCa                                                                                                                                >  2:833319/1‑150 (MQ=255)
                      gccAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAACGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGACCGAGTGGCCTTCGCCTGGCAc                                                                                                                               >  1:399597/1‑150 (MQ=255)
                      gccAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCAc                                                                                                                               >  2:775897/1‑150 (MQ=255)
                              gCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTg                                                                                                                                              >  2:498419/1‑127 (MQ=255)
                              gCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTg                                                                                                                                              <  1:498419/127‑1 (MQ=255)
                              gCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACgg                                                                                                                       >  1:596502/1‑150 (MQ=255)
                                     gcccTCCGCTGGCACCAGCGTGGCCAGGGTGTTCAAAGCCGCGGGGCATTCAGGGGGCGGGTGCGTGCGGCTAAAGTCGGTGCCGGTATCCTCGTCGCGACCGCGCTGGGCGATCGATTGGCGTTCGCCTGGCACGAACACGGCGAacac                                                                                                                <  2:316529/148‑1 (MQ=255)
                                     gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCt                                                                                                                                    <  2:9689/130‑1 (MQ=255)
                                     gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCt                                                                                                                                    >  1:9689/1‑130 (MQ=255)
                                       ccTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAAcacca                                                                                                              >  2:611999/1‑150 (MQ=255)
                                         tCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGaa                                                                                                                            <  2:1085963/134‑1 (MQ=255)
                                         tCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGaa                                                                                                                            >  1:1085963/1‑134 (MQ=255)
                                                         tgACCAGGGCCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCATCGTAGCGACCGCCATGGGCAAAAGATTAAACTTCGAATGGCACAAAAGCGGCCAACACCACGCCCGTgtggtagtggg                                                                                            >  2:851181/1‑149 (MQ=255)
                                                           aCCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTaa                                                                                          <  1:1133304/150‑1 (MQ=255)
                                                             cAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGc                                                                                                                      >  1:518789/1‑120 (MQ=255)
                                                             cAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGc                                                                                                                      <  2:518789/120‑1 (MQ=255)
                                                             cAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAAcc                                                                                        >  2:562408/1‑150 (MQ=255)
                                                               gggTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCAc                                                                                      >  1:1164319/1‑150 (MQ=255)
                                                               gggTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCAc                                                                                      <  1:7019/150‑1 (MQ=255)
                                                                ggTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACg                                                                                     >  1:1250944/1‑150 (MQ=255)
                                                                    ttCAAATCCGTGGAGAAACCAGTGACCTAGCGCGCCCGGCAAAAGTGCACACCGAGATCGGCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCGGTCACGAACACGGCGAACACCACGCCCGTGTGGGAGTGGTGACCACGCAGc                                                                                 >  2:180665/1‑150 (MQ=255)
                                                                            cGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAg                                                                         <  1:112165/150‑1 (MQ=255)
                                                                             gTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACg                                                                                                                        <  2:1153062/102‑1 (MQ=255)
                                                                             gTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACg                                                                                                                        >  1:1153062/1‑102 (MQ=255)
                                                                             gTGGAGAAACCAGTGGCCGGGCGCGCCCCGCCAAAGACGGCGCCGATATCGTCGTAGCGACCGCGCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACTGCGAACACCACGCCTATGAGGTAAAGGAAACCACGCAGCTCGCCGAgg                                                                        >  2:515905/1‑150 (MQ=255)
                                                                                  gaAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                                   <  1:35328/150‑1 (MQ=255)
                                                                                  gaAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGAACGAACACGGCGAACACCACGCCCGTGTGGTAGAATTAACCACGCAGCTCGCCGAGGTCGaa                                                                   <  1:359312/150‑1 (MQ=255)
                                                                                    aaCCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGt                                                                 >  2:814806/1‑150 (MQ=255)
                                                                                     aCCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTa                                                                <  2:1164319/150‑1 (MQ=255)
                                                                                      ccAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTAc                                                               <  1:1013197/150‑1 (MQ=255)
                                                                                       cAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTAGTAACCACGCAGCTCGCCGACGTCGAAGTACa                                                              >  1:488342/1‑150 (MQ=255)
                                                                                        aGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAAcc                                                                                        <  1:592729/123‑1 (MQ=255)
                                                                                        aGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAAcc                                                                                        >  2:592729/1‑123 (MQ=255)
                                                                                         gTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCgtgt                                                                                                     <  2:26416/109‑1 (MQ=255)
                                                                                         gTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCgtgt                                                                                                     >  1:26416/1‑109 (MQ=255)
                                                                                         gTGGCCGGGCGCGCCAGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGc                                                            >  2:784210/1‑150 (MQ=255)
                                                                                          tGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCg                                                           <  1:851181/150‑1 (MQ=255)
                                                                                             ccGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCa                                                        <  1:833319/150‑1 (MQ=255)
                                                                                             ccGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCa                                                        <  2:1157840/150‑1 (MQ=255)
                                                                                                    cgcCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCgg                                                 <  1:562408/150‑1 (MQ=255)
                                                                                                    cccccGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCgg                                                 <  2:809841/148‑1 (MQ=255)
                                                                                                         gtcaaacgTCTGCGCCTATATCGTCGTAGCAGCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAgc                                            <  1:180665/143‑1 (MQ=255)
                                                                                                         ggCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAgc                                            <  1:814806/150‑1 (MQ=255)
                                                                                                          gCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAgcg                                           <  1:775897/150‑1 (MQ=255)
                                                                                                           ccAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAgcgg                                          >  2:306165/1‑149 (MQ=255)
                                                                                                           ccAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAgcgc                                          <  2:399597/150‑1 (MQ=255)
                                                                                                             aaaGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGAGCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGa                                        >  1:802827/1‑150 (MQ=255)
                                                                                                              aaGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGAc                                       <  2:1250944/150‑1 (MQ=255)
                                                                                                                        cGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGACCACCACCCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGCCCGGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGAt                             <  1:1066683/150‑1 (MQ=255)
                                                                                                                        cGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGAt                             >  2:314006/1‑150 (MQ=255)
                                                                                                                         gATATCGTCGGAGGGACCGGGCGGGGAAAAAGATTGGGCTTTGCGTGGGACGGAGAAGGCGAACATCAGGCCCGTGTGGTAGTGGTAACCACGCAGGACGGCGAGGTCGAAGTAGAGCGGCAGGTCGGGGGAGGGCGGCGCCAGGGGAtg                            >  2:41563/1‑149 (MQ=255)
                                                                                                                                   aatcgatcGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGGACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGcc                  >  1:402534/8‑150 (MQ=255)
                                                                                                                                    atcgatcGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAGCCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGAc                                       >  1:862814/7‑128 (MQ=255)
                                                                                                                                    atcgatcGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGAc                                       <  2:862814/122‑1 (MQ=255)
                                                                                                                                     tcGATCGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcag                >  1:1061304/2‑150 (MQ=255)
                                                                                                                                     gCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcag                <  1:784210/150‑1 (MQ=255)
                                                                                                                                       gACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagca              >  1:710494/1‑150 (MQ=255)
                                                                                                                                        aCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                                   >  1:1284439/1‑96 (MQ=255)
                                                                                                                                        aCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                                   <  2:1284439/96‑1 (MQ=255)
                                                                                                                                        aCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagcag             <  2:380438/150‑1 (MQ=255)
                                                                                                                                         tcGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAgg            <  2:402534/149‑1 (MQ=255)
                                                                                                                                          cGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGt           >  1:175646/1‑150 (MQ=255)
                                                                                                                                           gCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTc          >  2:489547/1‑150 (MQ=255)
                                                                                                                                            cccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCa         <  2:175646/150‑1 (MQ=255)
                                                                                                                                            cccTGGGCGATCGATTGGCCTTCGACTGGTACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCa         <  2:802827/150‑1 (MQ=255)
                                                                                                                                             ccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCAt        >  1:760414/1‑150 (MQ=255)
                                                                                                                                               tGGGCGATCGATTGGCTTTCGCTTAGGACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCAtac      <  2:488342/150‑3 (MQ=255)
                                                                                                                                                   cgatcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAgcg  <  1:37587/150‑1 (MQ=255)
                                                                                                                                                    gatcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGAc                                       <  2:707914/112‑1 (MQ=255)
                                                                                                                                                    gatcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGAc                                       >  1:707914/1‑112 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCG  >  NC_002947/5559524‑5559820

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