Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,673 G→C 10.1% intergenic (+21/+42) 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,6730GC10.1% 312.1 / 19.9 106intergenic (+21/+42)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 = 9.99e-01

Only 100 of 110 total aligned reads displayed.

GCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTC  >  NC_000913/3530525‑3530822
                                                                                                                                                      |                                                                                                                                                     
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                                                                                                                                              <  2:881521/150‑1 (MQ=255)
        gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCg                                                                                                                                              >  1:345768/1‑150 (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:127096/150‑1 (MQ=255)
                    gAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACt                                                                                                                                  <  2:1697237/150‑1 (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)
                                   tgtgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATacg                                                                                                                   <  2:1216978/150‑1 (MQ=255)
                                         tACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGt                                                                                                             <  2:1368587/150‑1 (MQ=255)
                                             gCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACt                                                                                                         <  1:511072/150‑1 (MQ=255)
                                             gCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACt                                                                                                         <  2:646688/150‑1 (MQ=255)
                                              cGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCTTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCTTGATTTTACGACTAATACACCTGATGCGCTTTGTCACTc                                                                                                        <  1:258305/150‑1 (MQ=255)
                                                 tAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                                     >  2:681598/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                            >  1:1419029/1‑138 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                                >  2:1117980/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                                >  1:952245/1‑150 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                                <  2:1381604/150‑1 (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                                                                                         <  2:442486/150‑1 (MQ=255)
                                                             gTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAtt                                                                                         >  1:1703189/1‑150 (MQ=255)
                                                             gTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAtt                                                                                         >  2:1664842/1‑150 (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)
                                                                                        ccGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATc                                                              >  1:52338/1‑150 (MQ=255)
                                                                                           caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGCGGGCTCTGCCGGATGTGACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgcg                                                           >  1:130653/1‑150 (MQ=255)
                                                                                           caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgcg                                                           >  2:936715/1‑150 (MQ=255)
                                                                                              caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCa                                                        <  2:1081909/150‑1 (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                                                <  2:52338/150‑1 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                                >  2:550550/1‑150 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                                >  1:552068/1‑150 (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)
                                                                                                                                     aggTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGc                 <  2:1161390/148‑1 (MQ=255)
                                                                                                                                     aaatCCGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGc                 >  1:1457009/4‑150 (MQ=255)
                                                                                                                                      aatCCGGC‑‑AGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCt                >  1:1605543/3‑150 (MQ=255)
                                                                                                                                       gTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtg               <  1:177696/150‑1 (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)
                                                                                                                                                      ggCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGCTCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTtttc  >  1:1266189/1‑150 (MQ=255)
                                                                                                                                                      ggCGCTGCCGGATTTTAAGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCGTCGTCCCCAGTTTTCCATCCTTCGCAGCTGGGAGTAATTtttc  <  2:255819/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGC‑‑AGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTC  >  NC_000913/3530525‑3530822

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