Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,673 G→C 5.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,6730GC 5.1% 224.1 / 3.9 79intergenic (+21/+42)hslO/yhgEheat shock protein Hsp33/DUF4153 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/46);  new base C (2/2);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

AGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTC  >  NC_000913/3530524‑3530822
                                                                                                                                                     |                                                                                                                                                     
aGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAgg                                                                                                                                                       <  1:664390/150‑1 (MQ=255)
 gCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAggg                                                                                                                                                      <  1:104242/150‑1 (MQ=255)
      cTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGGCTCCGGAAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTg                                                                                                                                                 <  1:789692/150‑1 (MQ=255)
      cTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAAAGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTg                                                                                                                                                 >  2:1181023/1‑150 (MQ=255)
         gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                              <  2:586516/150‑1 (MQ=255)
         gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCg                                                                                                                                              >  2:364149/1‑150 (MQ=255)
         gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTGCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCg                                                                                                                                              <  2:388475/150‑1 (MQ=255)
          cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCgg                                                                                                                                             <  1:1165607/150‑1 (MQ=255)
             aaGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACACCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGAtt                                                                                                                                          >  2:59180/1‑150 (MQ=255)
                 tGGCGAAATTGCCATGCATTGTGCTTACTTCGGAAACCACTATCTGCTCAATGCGATGGATATTGCGGAAATCCGCAACAACGCGGCTCCGGCAGATCCGCAAGTTCATAAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                                      <  1:950972/150‑1 (MQ=255)
                  ggCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCTCAACAACGCGTCTCAGGCAGATCCGCAAGTTCATTAATGAGTATTTCCGGCAGAGACGCCACTGCCGGATTTTaag                                                                                                                                     >  2:969784/1‑148 (MQ=255)
                    cGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATAAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGAc                                                                                                                                   <  1:806625/150‑1 (MQ=255)
                     gaaagggACATGCATTGTGATTACTGCGGTAACCACTCTCTGTTCAATCCGATGGATATTGCTGAAAGCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACt                                                                                                                                  <  1:185803/144‑1 (MQ=255)
                       aaTTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTaa                                                                                                                                <  2:1130157/150‑1 (MQ=255)
                         ttGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGCAAATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATa                                                                                                                              >  2:814795/1‑150 (MQ=255)
                           gACATGCATTGTGATTACTGCGGGAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATaca                                                                                                                            <  1:681987/150‑1 (MQ=255)
                                gCATTGTGATTACTGCGGTTACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGa                                                                                                                       <  2:671988/150‑1 (MQ=255)
                                    tgtgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGTCTAATACACCTGATgcg                                                                                                                   <  1:699146/150‑1 (MQ=255)
                                     gtgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCTCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATgcgc                                                                                                                  <  2:921493/150‑1 (MQ=255)
                                       taTTAATGCGGTAAGATCTATGTGTTAAATGCGATGTATTTTTCTTAATTCCGCTACAACGCGTCTCCGGCAGATCCTCAAGTTAATTAATGAGTATGTCCGGAAGAGAGGGCTCTGCAGGATTTAACGACTATTACACCTGATGCGCtt                                                                                                                <  2:221488/149‑1 (MQ=255)
                                        aTTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCttt                                                                                                               <  1:846199/150‑1 (MQ=255)
                                          tACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGt                                                                                                             <  1:1188180/150‑1 (MQ=255)
                                           aCTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                            <  2:1055093/150‑1 (MQ=255)
                                           aCTGCGGTAAACACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGATACAACGCGTCTCCTGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGTGCTTTGTc                                                                                                            <  2:542613/150‑1 (MQ=255)
                                            cTGCGGTAACCACTATCTGTTCATTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCa                                                                                                           <  2:629301/150‑1 (MQ=255)
                                             tGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCAc                                                                                                          <  1:491321/150‑1 (MQ=255)
                                             tGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCAc                                                                                                          >  2:821258/1‑150 (MQ=255)
                                              gCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACt                                                                                                         <  2:541484/150‑1 (MQ=255)
                                               cGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTc                                                                                                        <  1:1279579/150‑1 (MQ=255)
                                                 gTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTATACGACTAATACACCTGATGCGCTATGTCACTCAt                                                                                                      >  1:1209195/1‑150 (MQ=255)
                                                    aCCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCa                                                                                                   <  2:774502/150‑1 (MQ=255)
                                                    aCCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGAAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCa                                                                                                   <  1:749091/150‑1 (MQ=255)
                                                    aCCACTATCTGTTCAATGCGATGGATATTCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAg                                                                                                  >  1:1183969/1‑150 (MQ=255)
                                                     ccACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTATTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAg                                                                                                  >  2:318220/1‑150 (MQ=255)
                                                       aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTGATTAATGAGTATCTCCGGCAGAGAGTGCTCTGCCGGATTTTACGAATTATACACCTGATGCGGTTTGTCACTCATTCAGAt                                                                                                >  2:795548/1‑150 (MQ=255)
                                                       aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACATGATGCGCTATTTCAAACATTCatat                                                                                                >  2:978085/1‑150 (MQ=255)
                                                           tCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTCATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTcac                                                                                            <  1:621554/150‑1 (MQ=255)
                                                           tCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTcac                                                                                            <  1:896486/150‑1 (MQ=255)
                                                             tGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAt                                                                                          >  2:343283/1‑150 (MQ=255)
                                                               ttCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTc                                                                                        <  1:1124098/150‑1 (MQ=255)
                                                                tCAATGCTATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCa                                                                                       <  2:182181/150‑1 (MQ=255)
                                                                   aTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGcc                                                                                    >  2:1238231/1‑150 (MQ=255)
                                                                    tGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCAGTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCt                                                                                   >  2:1082970/1‑150 (MQ=255)
                                                                      cGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCttt                                                                                 >  1:912814/1‑150 (MQ=255)
                                                                       gATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTAAGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTc                                                                                <  2:611867/150‑1 (MQ=255)
                                                                         tGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCAc                                                                              <  1:1140207/150‑1 (MQ=255)
                                                                         tGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCAc                                                                              >  2:1154599/1‑150 (MQ=255)
                                                                              atTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATcc                                                                         <  2:111695/150‑1 (MQ=255)
                                                                                  cTGAAATCCGCAACAACGCGTCTCCGGCAGCTCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGAGGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCGTCCACTa                                                                     >  1:615806/1‑150 (MQ=255)
                                                                                      aaTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAg                                                                 >  2:7788/1‑150 (MQ=255)
                                                                                         ccGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACTACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATc                                                              <  1:364149/150‑1 (MQ=255)
                                                                                         ccGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTAAGCCTTTCACCATCCACTAAAAGATc                                                              >  2:1063475/1‑150 (MQ=255)
                                                                                          cGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcg                                                             <  2:522707/150‑1 (MQ=255)
                                                                                          cGCAACAACGCGTATCCTGCATATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGATCTGCCGGATTTTACGACTAATACACCTGATGCTCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcg                                                             <  1:1154291/150‑1 (MQ=255)
                                                                                            caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTTTGTCCGGGAGAGAGGGCTCTGGCGGATTTTACGACTATTACACCTGCGGCGCTTTGTCGCTCTTTCAGATTGACTGTCGGCCGTTCCCCGTCGACTACATGATagcg                                                           >  2:788981/1‑150 (MQ=255)
                                                                                               caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCCGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCa                                                        <  2:708091/150‑1 (MQ=255)
                                                                                               caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGACCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCa                                                        <  1:79056/150‑1 (MQ=255)
                                                                                                 ccgcgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACg                                                      <  2:1172827/149‑1 (MQ=255)
                                                                                                      gctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCGCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTa                                                 <  2:111787/149‑1 (MQ=255)
                                                                                                        tcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAgg                                               <  2:487119/150‑1 (MQ=255)
                                                                                                           ggCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCtt                                            <  1:1181023/150‑1 (MQ=255)
                                                                                                             cAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCg                                          <  2:310361/150‑1 (MQ=255)
                                                                                                                aTCCGCAAGTTGATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCGTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTcc                                       <  2:826055/150‑1 (MQ=255)
                                                                                                                aTCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTcc                                       <  1:1251670/150‑1 (MQ=255)
                                                                                                                      aaGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGtt                                 >  2:1118452/1‑150 (MQ=255)
                                                                                                                         ttCATTAATGAGTATGTCCGGCAGAGAGGGCGCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTc                              <  2:657819/150‑1 (MQ=255)
                                                                                                                           cATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCCCGCTTTAGGCTTCGTCCCCAGTTTTCCa                            >  2:546578/1‑150 (MQ=255)
                                                                                                                           cATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCa                            <  2:566999/150‑1 (MQ=255)
                                                                                                                               aaTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt                        >  1:1037152/1‑150 (MQ=255)
                                                                                                                               aaTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt                        <  2:1150335/150‑1 (MQ=255)
                                                                                                                                aTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCAATGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGTTCGCGCCACGCTTTAGGCTTCGGCCCCAGTTTTCCATCCtg                       >  1:28202/1‑149 (MQ=255)
                                                                                                                                  gAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCCCTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCg                     >  1:250166/1‑150 (MQ=255)
                                                                                                                                  gAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCGCCAGTTTTCCATCCTTCg                     >  1:190180/1‑150 (MQ=255)
                                                                                                                                      aaatCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCAGATGCTCTTTGTCACTCATTAAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCTCGACACGCTTTAGGCTTCGTCCACAGGTTTCCATCCTTCGCAGc                 >  1:493430/4‑150 (MQ=255)
                                                                                                                                       tGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCt                >  1:609658/1‑150 (MQ=255)
                                                                                                                                          ccGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCCCCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg             >  2:873841/1‑150 (MQ=255)
                                                                                                                                          ccGGCAGAGAGGTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCATCtgtg             <  1:59180/150‑1 (MQ=255)
                                                                                                                                           cGGCAGAGCGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCCCTCATTCAGATTCACATTCAGGCTTTCACCATCCACTACAAGCTCGCGCCACGCTTTAGGCGTCGTCCCCAGGTGTCCATCCTTCGCAGCtgggg            >  1:349368/1‑147 (MQ=255)
                                                                                                                                           cGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGa            <  2:896939/150‑1 (MQ=255)
                                                                                                                                            ggCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAg           <  1:446670/150‑1 (MQ=255)
                                                                                                                                                 agagGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAAtt      >  2:538347/1‑150 (MQ=255)
                                                                                                                                                    gggCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAAttttt   <  2:152988/150‑1 (MQ=255)
                                                                                                                                                     ggCTCTGCCGGATTTTACGACTAATTCACCTGATGCGCGTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTtttc  <  2:808921/150‑1 (MQ=255)
                                                                                                                                                     ggCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAAttttta  >  2:294613/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTC  >  NC_000913/3530524‑3530822

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