Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,667 T→A 6.6% Q294L (CAG→CTG)  hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,6670TA 6.6% 218.9 / 6.5 77Q294L (CAG→CTG) hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base A (3/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCC  >  NC_002947/5559518‑5559816
                                                                                                                                                     |                                                                                                                                                     
cATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCt                                                                                                                                                       <  2:706664/150‑1 (MQ=255)
    gCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGc                                                                                                                                                   <  2:579343/150‑1 (MQ=255)
      cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGACAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcga                                                                                                                                                 <  2:987201/150‑1 (MQ=255)
       aCTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat                                                                                                                                                <  2:221718/150‑1 (MQ=255)
       aCTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat                                                                                                                                                <  1:683334/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                                                                                                                                >  1:9689/1‑130 (MQ=255)
                                           gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCt                                                                                                                                <  2:9689/130‑1 (MQ=255)
                                             ccTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAAcacca                                                                                                          >  2:611999/1‑150 (MQ=255)
                                               tCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGaa                                                                                                                        >  1:1085963/1‑134 (MQ=255)
                                               tCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGaa                                                                                                                        <  2:1085963/134‑1 (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:7019/150‑1 (MQ=255)
                                                                     gggTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCAc                                                                                  >  1:1164319/1‑150 (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                                                                                    >  2:592729/1‑123 (MQ=255)
                                                                                              aGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAAcc                                                                                    <  1:592729/123‑1 (MQ=255)
                                                                                               gTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCgtgt                                                                                                 >  1:26416/1‑109 (MQ=255)
                                                                                               gTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCgtgt                                                                                                 <  2:26416/109‑1 (MQ=255)
                                                                                               gTGGCCGGGCGCGCCAGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGc                                                        >  2:784210/1‑150 (MQ=255)
                                                                                                tGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCg                                                       <  1:851181/150‑1 (MQ=255)
                                                                                                   ccGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCa                                                    <  2:1157840/150‑1 (MQ=255)
                                                                                                   ccGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCa                                                    <  1:833319/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)
                                                                                                                                                     |                                                                                                                                                     
CATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCC  >  NC_002947/5559518‑5559816

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: ‑