Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,018,200 A→T 39.7% intergenic (+76/+55) prfC → / ← potF‑I peptide chain release factor 3/putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,018,2000AT39.7% 193.3 / 261.0 226intergenic (+76/+55)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (70/65);  new base T (48/41);  total (120/106)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-02

Only 100 of 240 total aligned reads displayed.

CCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCT  >  NC_002947/1018051‑1018349
                                                                                                                                                     |                                                                                                                                                     
 ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACCAGACCTGGAGGTCAt                                                                                                                                                      <  2:60900/150‑1 (MQ=255)
 ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGAAACACCACACCAGGAGGTCAt                                                                                                                                                      >  1:4187086/1‑150 (MQ=255)
  cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                                     <  1:4156672/150‑1 (MQ=255)
  cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                                     <  2:3659671/150‑1 (MQ=255)
   gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                                    <  1:933833/150‑1 (MQ=255)
   gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                                    <  1:1336853/150‑1 (MQ=255)
   gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                                    <  1:2024766/150‑1 (MQ=255)
     gTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACACGAAAAAACACAACAGGAGGTGAtgttt                                                                                                                                                  >  1:2500803/1‑148 (MQ=38)
      tCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGAt                                                                                                                                                 <  2:4550527/150‑1 (MQ=25)
        aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGAttc                                                                                                                                               <  1:2224136/150‑3 (MQ=25)
        aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                                               <  2:1140404/150‑1 (MQ=255)
         aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTcc                                                                                                                                              <  1:4234858/150‑1 (MQ=25)
          ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCTGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                             <  1:5726916/150‑1 (MQ=25)
          ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                             <  1:4465608/150‑1 (MQ=25)
          ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                             <  2:5276756/150‑1 (MQ=25)
          ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAAATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                             <  1:3575650/150‑1 (MQ=25)
            tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                                           <  1:6518948/150‑1 (MQ=21)
            tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                                           <  1:1058881/150‑1 (MQ=25)
            tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                                           <  1:3689300/150‑1 (MQ=25)
            tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                                           <  2:5984590/150‑1 (MQ=255)
            tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                                           <  2:6305568/150‑1 (MQ=255)
             ggCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATTCCtgtg                                                                                                                                          <  2:4987217/150‑1 (MQ=255)
              gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                                         <  1:3909089/150‑1 (MQ=25)
              gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                                         <  2:5391564/150‑1 (MQ=25)
              gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                                         <  1:5228149/150‑1 (MQ=25)
               cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTggg                                                                                                                                        <  2:4288593/150‑1 (MQ=255)
                gcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGa                                                                                                                                       <  1:1148110/150‑1 (MQ=21)
                gcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGa                                                                                                                                       <  2:5222437/150‑1 (MQ=21)
                 cTGATGGAAGAGCGCTGTCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAg                                                                                                                                      <  2:6215084/150‑1 (MQ=21)
                    aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGACATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCgg                                                                                                                                   <  2:1478025/150‑1 (MQ=255)
                    aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCgg                                                                                                                                   <  2:329153/150‑1 (MQ=21)
                    aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCGACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCgg                                                                                                                                   <  2:3750213/150‑1 (MQ=21)
                      ggAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGc                                                                                                                                 <  1:5837170/150‑1 (MQ=255)
                                    cGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                   >  2:3395042/1‑150 (MQ=255)
                                           cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAATAACCAGACCGGGAGGTCATGTGTTAGCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCTAg                                                                                                            >  1:5559072/1‑150 (MQ=38)
                                           cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                                         >  2:725028/1‑121 (MQ=18)
                                                tCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTcc                                                                                                                                              <  1:4013959/111‑1 (MQ=17)
                                                 ccGCGCTACCCGCGAGCACAATTACGGCCTTAACTCGGGGCTGCTTTGCAGCCCTTCGCGGGAACGCCCGCTCCCACCGGAATCACCTGACCTGGAGGTCATGTGATACGTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGcc                                                                                                      <  2:3200732/150‑1 (MQ=25)
                                                 ccGCGCGACCCGCGAGCACCATTAAGGCCTGACAGCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAAGCCCCCGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGcc                                                                                                      <  2:5448317/150‑1 (MQ=37)
                                                   gcgcTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCGGGAGGTCATGTGATACCGGTGGGAGCGGGCGGGCCCGCGCAGGGCCGCTAGGCGGcccc                                                                                                    >  2:5956015/1‑150 (MQ=255)
                                                               aGCACCCTTAAGGCCTGAACTCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGAATCACCAGACCTGCAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTcc                                                                                        <  2:411764/150‑1 (MQ=255)
                                                                gCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCACGAATCACCAGCCCGGGCGGGCCTGGGATACCGGGGGGGGCGGGCGTGCCCGCGAAGGGCCGCTCGGCGGCCCCATTTCAATTTCCg                                                                                       >  2:3667781/1‑150 (MQ=37)
                                                                     ctgccggCCTGACCGCGGGGCGGCTTGGCAGCCCTTCGCGGGCACGCCCGCTCCCCCCGGAATCCCCCGACCTGGAGGTCATGGGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCttt                                                                                  <  2:6655559/145‑1 (MQ=25)
                                                                     aTTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCttt                                                                                  <  2:3187837/150‑1 (MQ=255)
                                                                     aTTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCttt                                                                                  <  2:2161853/150‑1 (MQ=255)
                                                                            ccTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTg                                                                           >  2:3904847/1‑150 (MQ=255)
                                                                                      gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGTTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTc                                                                 <  2:3091510/150‑1 (MQ=255)
                                                                                       ggCTGCGTTGCAGCCCTTCGCGGGCACGCCCGCTCCCAGAGGAATCACCAGACCTGGAGGTCATGTGATACCTGGGGAGCGGGGGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCTATTTCCGGCTTTGATACTGGTCCAGATTcgc                                                               <  2:2577216/150‑1 (MQ=37)
                                                                                        gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                                               <  2:4736576/150‑1 (MQ=255)
                                                                                                         cgcgGGCACGCCCGCTCCCACAGTAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGAt                                              <  2:6542272/150‑1 (MQ=255)
                                                                                                          gcgGCCCCGCCCGCCCCCACCGCAATCACCTCACTTGTTGGTCTTGTTATACCTGGGGGCGGGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATc                                             <  1:4709895/150‑1 (MQ=255)
                                                                                                             ggCACGCCCGCTCCCACCGGAATCACCCGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTc                                          <  2:2463116/150‑1 (MQ=255)
                                                                                                                   cccGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGc                                    >  1:788946/1‑150 (MQ=21)
                                                                                                                      gCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATc                                 >  2:1990312/1‑150 (MQ=255)
                                                                                                                          ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCCTCGCCt                             >  1:4649991/1‑150 (MQ=21)
                                                                                                                          ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                             >  2:3145004/1‑150 (MQ=25)
                                                                                                                          ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                             >  2:3708169/1‑150 (MQ=255)
                                                                                                                          ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                             >  1:553338/1‑150 (MQ=255)
                                                                                                                           ccccggAAGCCCCAGACCTGGAGGTCATGGGATACCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGGACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            <  2:3004570/146‑1 (MQ=255)
                                                                                                                           caccGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            <  2:4032221/150‑1 (MQ=255)
                                                                                                                           cacaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCTAGGCGGCCCCATTTCTATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            >  2:6394296/1‑150 (MQ=21)
                                                                                                                           cacaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            >  2:5528956/1‑150 (MQ=25)
                                                                                                                           cacaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            >  1:6228264/1‑150 (MQ=25)
                                                                                                                           cacaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            >  2:26212/1‑150 (MQ=25)
                                                                                                                           cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                            >  1:5074419/1‑150 (MQ=255)
                                                                                                                            acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                                               <  1:3884252/114‑1 (MQ=17)
                                                                                                                            acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCTATCATCGCCGGAAACGCCTcc                           >  2:3120748/1‑150 (MQ=21)
                                                                                                                            acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGcc                                       >  1:2625232/1‑138 (MQ=21)
                                                                                                                            acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc                           >  1:4284924/1‑150 (MQ=25)
                                                                                                                            acaGGTATCACATGACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc                           >  2:1054742/1‑150 (MQ=25)
                                                                                                                            acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATTCGGTCGATCTTCGCCGGCATCGCCTcc                           >  1:6664161/1‑150 (MQ=255)
                                                                                                                            acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGCTCTTCGCCGGCATCGCCTcc                           >  2:2393872/1‑150 (MQ=255)
                                                                                                                              aGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCGTTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg                         >  2:4339061/1‑150 (MQ=18)
                                                                                                                              aGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg                         >  2:1939763/1‑150 (MQ=255)
                                                                                                                                 taTCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcgc                      >  1:5465222/2‑150 (MQ=25)
                                                                                                                                 aaTCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcgc                      <  2:4187086/150‑1 (MQ=255)
                                                                                                                                   tcacatgACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCAt                    >  2:5148439/7‑150 (MQ=37)
                                                                                                                                     acatgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCAACGCCTCCAGCGCATAc                  >  1:2721954/5‑150 (MQ=37)
                                                                                                                                     aCCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTcc                                                                                        >  1:1299966/1‑80 (MQ=255)
                                                                                                                                       atgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAg                >  1:6118112/3‑150 (MQ=37)
                                                                                                                                       atgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGc                                        <  2:486544/124‑1 (MQ=25)
                                                                                                                                       atgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAg                >  2:1282582/3‑150 (MQ=37)
                                                                                                                                       atgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAg                >  1:4840135/3‑150 (MQ=37)
                                                                                                                                       atgACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgcg                                                              >  2:421317/3‑104 (MQ=21)
                                                                                                                                        tgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGc               >  1:778173/2‑150 (MQ=38)
                                                                                                                                        aGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGc               <  1:4511191/150‑1 (MQ=255)
                                                                                                                                        aGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGc               >  2:2786105/1‑150 (MQ=255)
                                                                                                                                            cTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTg                                                                           <  2:3838077/86‑1 (MQ=255)
                                                                                                                                            cTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGc           >  1:6856562/1‑150 (MQ=255)
                                                                                                                                             tGGAGGGCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGcc          >  2:5617649/1‑150 (MQ=255)
                                                                                                                                               gAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCa        >  1:2981743/1‑150 (MQ=255)
                                                                                                                                               gAGGTCATGTGATACCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCa        >  2:6440439/1‑150 (MQ=255)
                                                                                                                                                aGGTCATGTGATCCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCAt       <  2:4467672/150‑1 (MQ=255)
                                                                                                                                                  gTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCa     >  2:1350791/1‑150 (MQ=255)
                                                                                                                                                  gTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCa     >  2:3880113/1‑150 (MQ=255)
                                                                                                                                                   tCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGACTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCAc    >  1:1769204/1‑150 (MQ=255)
                                                                                                                                                     aTGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGCATACGGTCGTTCGTCGCCGGCATCGCCTCCCGCGCGTACTGCGTGCCCCTCAccg  >  1:4904393/1‑149 (MQ=38)
                                                                                                                                                     aTGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCt  >  2:278296/1‑150 (MQ=255)
                                                                                                                                                     aTGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCt  >  1:1410027/1‑150 (MQ=255)
                                                                                                                                                     aTGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCAATTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCATGCCCATCACCt  >  1:3144582/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCT  >  NC_002947/1018051‑1018349

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