Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,671 A→C 5.7% intergenic (+19/+44) 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,6710AC 5.7% 218.3 / 2.2 80intergenic (+19/+44)hslO/yhgEheat shock protein Hsp33/DUF4153 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/48);  new base C (3/2);  total (30/50)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ATAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATT  >  NC_000913/3530522‑3530818
                                                                                                                                                     |                                                                                                                                                   
aTAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCagaga                                                                                                                                                     >  1:492840/1‑150 (MQ=255)
aTAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCAGCAAGTTCATTATTGAGTATGTCCGGCagaga                                                                                                                                                     >  2:319107/1‑150 (MQ=255)
 tAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCagagag                                                                                                                                                    >  1:1188903/1‑150 (MQ=255)
  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                                                                          >  2:1154599/1‑150 (MQ=255)
                                                                           tGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCAc                                                                          <  1:1140207/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
ATAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATT  >  NC_000913/3530522‑3530818

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: