Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,672 G→C 10.1% intergenic (+20/+43) hslO → / ← yhgE heat shock protein Hsp33/DUF4153 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,530,6720GC10.1% 314.3 / 19.8 106intergenic (+20/+43)hslO/yhgEheat shock protein Hsp33/DUF4153 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/53);  new base C (6/5);  total (48/58)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 108 total aligned reads displayed.

GCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTT  >  NC_000913/3530525‑3530821
                                                                                                                                                     |                                                                                                                                                     
gCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAggg                                                                                                                                                     <  2:400174/150‑1 (MQ=255)
  aTCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGct                                                                                                                                                   <  1:880180/150‑1 (MQ=255)
     cTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCAATAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTg                                                                                                                                                >  1:738719/1‑150 (MQ=255)
       ggCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGcc                                                                                                                                              >  2:177696/1‑150 (MQ=255)
        gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGCCCTCTCTGCCg                                                                                                                                             <  2:1663960/150‑1 (MQ=255)
        gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCg                                                                                                                                             >  1:345768/1‑150 (MQ=255)
        gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCg                                                                                                                                             <  2:881521/150‑1 (MQ=255)
         cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCTACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCgg                                                                                                                                            <  2:468571/150‑1 (MQ=255)
         cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGCCCTCTCTGCCgg                                                                                                                                            <  1:796772/150‑1 (MQ=255)
         cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCgg                                                                                                                                            <  2:1345045/150‑1 (MQ=255)
         cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCgg                                                                                                                                            >  1:1725294/1‑150 (MQ=255)
         cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGCAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCgg                                                                                                                                            <  1:1009587/150‑1 (MQ=255)
           gAAGATGGCTTAATTGACATGCATTTTGATTTCTGGGTTAACCAGTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGGCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGCCCGCTCTGCCGGac                                                                                                                                          <  2:1157557/150‑2 (MQ=255)
           gAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGCCCTCTCTGCCGGac                                                                                                                                          <  1:1652435/150‑2 (MQ=255)
                    gAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACt                                                                                                                                 <  2:1697237/150‑1 (MQ=255)
                    gAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACt                                                                                                                                 <  2:127096/150‑1 (MQ=255)
                        ttGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGTTTTTACGACTAATa                                                                                                                             >  2:1141099/1‑150 (MQ=255)
                          gACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATaca                                                                                                                           >  2:1720359/1‑150 (MQ=255)
                           aCATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATacac                                                                                                                          <  1:415563/150‑1 (MQ=255)
                            aaTGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACAcc                                                                                                                         <  2:1286838/149‑1 (MQ=255)
                              tGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTg                                                                                                                       <  1:1529861/150‑1 (MQ=255)
                                gatagtgGTTACCGCAGTAATCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAAAGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTATGCCGGATTTTACGACTAATACACCTGAt                                                                                                                     >  2:335714/5‑150 (MQ=255)
                                 aTTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATg                                                                                                                    <  1:1304946/150‑1 (MQ=255)
                                  ttGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATgc                                                                                                                   <  1:161651/150‑1 (MQ=255)
                                   tgtgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATacg                                                                                                                  <  2:1216978/150‑1 (MQ=255)
                                         tACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGt                                                                                                            <  2:1368587/150‑1 (MQ=255)
                                             gCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACt                                                                                                        <  2:646688/150‑1 (MQ=255)
                                             gCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACt                                                                                                        <  1:511072/150‑1 (MQ=255)
                                              cGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCTTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCTTGATTTTACGACTAATACACCTGATGCGCTTTGTCACTc                                                                                                       <  1:258305/150‑1 (MQ=255)
                                                 tAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                                    >  2:681598/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                           <  2:1419029/138‑1 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                           >  1:1419029/1‑138 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                               >  1:952245/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                               >  1:604986/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                               <  2:1381604/150‑1 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                               >  2:1117980/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCGCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                               <  2:29599/150‑1 (MQ=255)
                                                          tCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTcac                                                                                           <  1:1669405/150‑1 (MQ=255)
                                                           cTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTcaca                                                                                          <  2:748088/150‑1 (MQ=255)
                                                            tGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAt                                                                                         >  2:1141554/1‑150 (MQ=255)
                                                            tGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGGCTCCGGCAGATCCGCTAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAt                                                                                         <  1:315073/150‑1 (MQ=255)
                                                             gTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAtt                                                                                        >  1:1703189/1‑150 (MQ=255)
                                                             gTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAtt                                                                                        <  2:442486/150‑1 (MQ=255)
                                                                 aaTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGc                                                                                    <  1:915949/150‑1 (MQ=255)
                                                                 aaTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGc                                                                                    <  2:1113900/150‑1 (MQ=255)
                                                                   tGCGTTGGATATTGCTGATATCCGCAACAACGCGTCTCCGGCAGATCCGAATGTTCATTAATGAGTTTGTGCGGCAGAGAGAGGGCTCTGCCGGATTTTACGACTAATACACATGATGCGCTTTGTCACTCATTCAGATTCACTTTCatg                                                                                    <  2:1511824/150‑3 (MQ=255)
                                                                      gATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGATTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTc                                                                               >  1:726129/1‑150 (MQ=255)
                                                                         ggATATTGCTGAAATCCTCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCAcc                                                                            <  1:407340/150‑1 (MQ=255)
                                                                           atatTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCAt                                                                          <  1:743502/150‑1 (MQ=255)
                                                                           atatTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCAt                                                                          <  1:1488493/150‑1 (MQ=255)
                                                                                gCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCCGCCTTTCACCATCCACt                                                                     >  1:256443/1‑150 (MQ=255)
                                                                                 cTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATctacta                                                                    >  1:938183/1‑150 (MQ=255)
                                                                                  tGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTAc                                                                   >  1:1018522/1‑150 (MQ=255)
                                                                                  tGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTAc                                                                   >  1:809926/1‑150 (MQ=255)
                                                                                     aaTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAg                                                                >  2:35257/1‑150 (MQ=255)
                                                                                        ccGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCTACTACAAGATc                                                             <  2:938183/150‑1 (MQ=255)
                                                                                           caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGCGGGCTCTGCCGGATGTGACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgcg                                                          >  1:130653/1‑150 (MQ=255)
                                                                                           caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgcg                                                          >  2:936715/1‑150 (MQ=255)
                                                                                            aacaacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCgcgc                                                         <  2:218818/150‑1 (MQ=255)
                                                                                              caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCa                                                       <  2:1081909/150‑1 (MQ=255)
                                                                                              caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACAAACCACTACAAGATCGCGCCa                                                       >  2:565475/1‑150 (MQ=255)
                                                                                                   cgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCtt                                                  <  1:1720359/150‑1 (MQ=255)
                                                                                                   cgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCtt                                                  >  1:1019537/1‑150 (MQ=255)
                                                                                                   cgTCTCCGGCAGATCCGCAACTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCtt                                                  >  1:1562463/1‑150 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                               >  1:552068/1‑150 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                               >  2:550550/1‑150 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                               <  1:1308142/150‑1 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                               <  1:335714/150‑1 (MQ=255)
                                                                                                        ccGGCATATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGc                                             <  1:48175/150‑1 (MQ=255)
                                                                                                        ccGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGc                                             <  2:604986/150‑1 (MQ=255)
                                                                                                        ccGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGc                                             <  1:1436663/150‑1 (MQ=255)
                                                                                                            cAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCg                                         <  1:442562/150‑1 (MQ=255)
                                                                                                            cAGATCCGAATGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCg                                         <  2:350065/150‑1 (MQ=255)
                                                                                                              gATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTc                                       <  1:1112507/150‑1 (MQ=255)
                                                                                                                tCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTccc                                     <  2:1677883/150‑1 (MQ=255)
                                                                                                                tCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTccc                                     >  1:729675/1‑150 (MQ=255)
                                                                                                                 ccGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCCTCCACTACAAGATCGCGCCCCGCTTTAGGCTTCGTcccc                                    >  1:1053789/1‑150 (MQ=255)
                                                                                                                 ccGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTcccc                                    >  1:742186/1‑150 (MQ=255)
                                                                                                                   gCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAg                                  <  2:1457430/150‑1 (MQ=255)
                                                                                                                   gCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAg                                  <  1:681598/150‑1 (MQ=255)
                                                                                                                      aGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCATCCCCAGatt                               >  2:1042090/1‑150 (MQ=255)
                                                                                                                        ttCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTc                             <  1:1263109/150‑1 (MQ=255)
                                                                                                                         tCATTAATGAGTATGGCCGGC‑‑AGAGAGGGCGCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGACTTcc                            <  2:1071440/150‑1 (MQ=255)
                                                                                                                          cATTAATGAGTATGTCCGGC‑‑AGTGAGGGCTCTGCCGGATTTTACGACTAATTCAGCTGATGCGCTTTGTCACTCTTTCAGATTCACTTTCAGGCTTTCTCCATGCACTATAATATCGCTCCTCGCTTTAGGCTTCGTCCCCAGTTTTGCa                           <  2:1182678/150‑1 (MQ=255)
                                                                                                                          cATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCa                           >  1:1055556/1‑150 (MQ=255)
                                                                                                                              aaTGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCCCCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt                       >  1:1573284/1‑150 (MQ=255)
                                                                                                                              aaTGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt                       >  2:202329/1‑150 (MQ=255)
                                                                                                                                      aatCCGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCt               >  1:1605543/3‑150 (MQ=255)
                                                                                                                                       gTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCGTTGTCCCTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtg              >  1:334788/1‑150 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg            >  1:766997/1‑150 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCt               <  1:996173/147‑1 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCggtg            >  1:490973/1‑150 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCGTCGCAGCtggg            >  1:484501/1‑148 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg            >  1:1432202/1‑150 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg            <  2:738719/150‑1 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg            >  2:1656687/1‑150 (MQ=255)
                                                                                                                                         ccGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg            >  1:784749/1‑150 (MQ=255)
                                                                                                                                           ggC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAg          <  1:1588958/150‑1 (MQ=255)
                                                                                                                                            gC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGt         >  2:650590/1‑150 (MQ=255)
                                                                                                                                                     gggCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAAttttt  >  2:750946/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTT  >  NC_000913/3530525‑3530821

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