Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,585 A→T 5.9% intergenic (‑38/+264) hisZ ← / ← hflC ATP phosphoribosyltransferase regulatory subunit/integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,5850AT 5.9% 198.2 / 3.7 69intergenic (‑38/+264)hisZ/hflCATP phosphoribosyltransferase regulatory subunit/integral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base T (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

AACTCTGGAACAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGAC  >  NC_002947/5560436‑5560726
                                                                                                                                                     |                                                                                                                                             
aaCTCTGGAACAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCcagca                                                                                                                                               <  2:220820/150‑1 (MQ=255)
   tctGGAACAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGGCATTTCCCCTCTGGTCCGGGGGGCTTGCCAGCGGGc                                                                                                                                            >  1:68133/1‑150 (MQ=255)
   tctGGAACAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGc                                                                                                                                            <  1:300173/150‑1 (MQ=255)
     tGGAACAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCaa                                                                                                                                          >  1:957156/1‑150 (MQ=255)
        aaCAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCTGCTGGCAAGcc                                                                                                                                       <  1:153195/150‑1 (MQ=255)
        aaCAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCTGCTGGCAAGcc                                                                                                                                       >  2:153195/1‑150 (MQ=255)
        aaCAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAAcc                                                                                                                                       <  1:291344/150‑1 (MQ=255)
           aGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTc                                                                                                                                    <  1:83453/150‑1 (MQ=255)
            ggTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCTGCTGGCaagc                                                                                                                                        >  2:835220/1‑143 (MQ=255)
            ggTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCTGCTGGCaagc                                                                                                                                        <  1:835220/145‑3 (MQ=255)
            ggTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCACTACTTCCTCGATGCCATCTGGCCGCACCCAGCGTTGTACCGTTGCTATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTcc                                                                                                                                   >  1:879046/1‑149 (MQ=255)
               ccAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGaa                                                                                                                                >  1:527155/1‑150 (MQ=255)
                 aaCACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGaaaa                                                                                                                              <  2:957156/150‑1 (MQ=255)
                 aaCACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGaaaa                                                                                                                              >  2:1210440/1‑150 (MQ=255)
                  acacCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAg                                                                                                                             >  1:1261910/1‑150 (MQ=255)
                   cacCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGc                                                                                                                            <  1:900564/150‑1 (MQ=255)
                   cacCTGGCCGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGc                                                                                                                            <  1:465909/150‑1 (MQ=255)
                    acCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGcc                                                                                                                           >  1:927807/1‑150 (MQ=255)
                    acCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGcc                                                                                                                           <  2:927807/150‑1 (MQ=255)
                     ccTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCg                                                                                                                          <  2:89819/150‑1 (MQ=255)
                      cTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCgg                                                                                                                         <  1:1210440/150‑1 (MQ=255)
                      cTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCgg                                                                                                                         >  1:454089/1‑150 (MQ=255)
                             cgcgcgATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATcgcg                                                                                                                  <  2:454089/150‑1 (MQ=255)
                              gcgcgATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCgcgc                                                                                                                 >  2:430046/1‑150 (MQ=255)
                                    tctcGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCAcagca                                                                                                           >  2:256122/1‑150 (MQ=255)
                                       ggATGCCCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCTGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCg                                                                                                        <  1:380135/149‑1 (MQ=255)
                                        gATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCAAGCGGTCTACCGTTGCCATTACGTCATTTCCCCTCTGGTCCGAGCGGCGTGCCAGCATGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCgg                                                                                                       <  1:24392/150‑1 (MQ=255)
                                         aTGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCCGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCggg                                                                                                      <  1:1131250/150‑1 (MQ=255)
                                         aTGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCa                                                                                                                                   <  1:444342/121‑1 (MQ=255)
                                         aTGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCa                                                                                                                                   >  2:444342/1‑121 (MQ=255)
                                           gcgcgcAGCCTCAGTGGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGgca                                                                                                    <  2:879046/150‑1 (MQ=255)
                                           gcgcgcAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGgca                                                                                                    <  2:527155/150‑1 (MQ=255)
                                                     tCACGTGGCAGTACTTCCTCGATGCCATCTGGCCTCTTCCGTTGGTCTACCGTTGCCAGTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATcc                                                                                          <  2:622336/150‑1 (MQ=255)
                                                      cAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATccc                                                                                         <  2:521137/150‑1 (MQ=255)
                                                       aGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGt                                                                                                                        >  2:169406/1‑118 (MQ=255)
                                                       aGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGt                                                                                                                        <  1:169406/118‑1 (MQ=255)
                                                            gCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAgcgcc                                                                                   >  2:413887/1‑149 (MQ=255)
                                                            gCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCa                                                                                   >  2:413801/1‑150 (MQ=255)
                                                                  cTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACagcagc                                                                                                         >  1:238853/1‑122 (MQ=255)
                                                                  cTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACagcagc                                                                                                         <  2:238853/122‑1 (MQ=255)
                                                                            gCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGaca                                                                   >  1:765355/1‑150 (MQ=255)
                                                                              cATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACAcc                                                                 >  2:1173512/1‑150 (MQ=255)
                                                                               aTCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACAccc                                                                >  2:853352/1‑150 (MQ=255)
                                                                                tCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCt                                                               <  2:198786/150‑1 (MQ=255)
                                                                                   ggCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCaa                                                            <  2:765355/150‑1 (MQ=255)
                                                                                      agcagcCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAAtt                                                         <  1:853352/150‑1 (MQ=255)
                                                                                        cagccagcGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGc                                                       <  1:295259/150‑1 (MQ=255)
                                                                                                ggTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAg                                               <  1:413801/150‑1 (MQ=255)
                                                                                                ggTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGCATTGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAg                                               <  1:413887/150‑1 (MQ=255)
                                                                                                      ccGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGAcc                                         >  1:416167/1‑150 (MQ=255)
                                                                                                         ttGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCaaa                                      >  1:189939/1‑150 (MQ=255)
                                                                                                            ccATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCg                                   >  2:95996/1‑150 (MQ=255)
                                                                                                              aTTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGAGAGCTGACCAAACCGTc                                 >  2:990017/1‑150 (MQ=255)
                                                                                                               ttACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCcagca                                                                                                                                               >  2:266195/1‑39 (MQ=255)
                                                                                                               ttACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCcagca                                                                                                                                               <  1:266195/39‑1 (MQ=255)
                                                                                                                   gCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCg                                                                                                        >  1:907090/1‑74 (MQ=255)
                                                                                                                   gCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCg                                                                                                        <  2:907090/74‑1 (MQ=255)
                                                                                                                   gCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCg                            >  1:1173224/1‑150 (MQ=255)
                                                                                                                            ccTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGa                                                                     >  2:944651/1‑100 (MQ=255)
                                                                                                                            ccTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGa                                                                     <  1:944651/100‑1 (MQ=255)
                                                                                                                              tctGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGAc                                          <  1:1161732/125‑1 (MQ=255)
                                                                                                                              tctGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGAc                                          >  2:1161732/1‑125 (MQ=255)
                                                                                                                                tGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGGACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGGGAGCTGACCAAACCGTCAACCGGCGCTGAGCCacg               >  1:123360/1‑148 (MQ=255)
                                                                                                                                    ccGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAAc           <  1:256122/150‑1 (MQ=255)
                                                                                                                                    ccGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAAc           <  1:1173512/150‑1 (MQ=255)
                                                                                                                                      gggCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCTAATCAACCt         >  1:784732/1‑150 (MQ=255)
                                                                                                                                       ggcggcTTGCCATCAGGCAAACGGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGGGCCAATCAACCtt        <  2:68133/150‑1 (MQ=255)
                                                                                                                                        gcggcTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTg       >  2:89686/1‑150 (MQ=255)
                                                                                                                                        gcggcTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTg       <  2:189939/150‑1 (MQ=255)
                                                                                                                                          ggcTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCa     <  1:430046/150‑1 (MQ=255)
                                                                                                                                            cTTGCCGGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGa   <  2:491564/150‑1 (MQ=255)
                                                                                                                                             ttGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGAc  <  1:19278/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
AACTCTGGAACAGGTCCAACACCTGGCGGCGCGCGATCTCGATGCGCGCAGCCTCAGGTGGCAGTACTTCCTCGATGCCATCTGGCAGCAGCCAGCGGTCTACCGTTGCCATTACGCCATTTCCCCTCTGGTCCGGGCGGCTTGCCAGCAGGCAAACCGTCAGAAAAGCCGGTATCGCGCACAGCAGCGGGCAGCACTGATCCCAGCGCAGCCACCGTACCCGACACCCTCAAATTGCCTTGCGAGCTGACCAAACCGTCAACCGGCGCTGAGCCAATCAACCTTGCAGAC  >  NC_002947/5560436‑5560726

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