Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,575,647 C→G 7.7% M53I (ATG→ATC eutS ← putative ethanol utilization carboxysome structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,575,6470CG 7.7% 264.7 / 9.3 78M53I (ATG→ATCeutSputative ethanol utilization carboxysome structural protein
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base G (2/4);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ATTTAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGCGT  >  NC_000913/2575500‑2575794
                                                                                                                                                   |                                                                                                                                                   
aTTTAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGatc                                                                                                                                                     >  3:543149/1‑148 (MQ=255)
  ttAATAGGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGatcat                                                                                                                                                   >  4:359403/1‑148 (MQ=255)
  ttAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGatcat                                                                                                                                                   >  3:398812/1‑148 (MQ=255)
  ttAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGatcat                                                                                                                                                   <  3:916068/148‑1 (MQ=255)
  ttAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAACGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGatcat                                                                                                                                                   >  3:368222/1‑148 (MQ=255)
     ataaGCGACCTAAACCGCTGACGGCTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGc                                                                                                                                                >  3:452896/1‑148 (MQ=255)
      taaGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCg                                                                                                                                              >  2:1056113/1‑149 (MQ=255)
       aaGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCg                                                                                                                                              <  4:710263/148‑1 (MQ=255)
       aaGCGACCTAAACCGCTGACGGTTGGGGATAACGCTGCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCg                                                                                                                                              <  4:452896/148‑1 (MQ=255)
          cGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTc                                                                                                                                          <  1:165107/149‑1 (MQ=255)
           gACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggtga                                                                                                                                          <  4:225513/148‑5 (MQ=255)
               tAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCg                                                                                                                                      >  4:63460/1‑148 (MQ=255)
                aaaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCggg                                                                                                                                    <  1:1056113/149‑1 (MQ=255)
                 aaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCggg                                                                                                                                    <  4:368222/148‑1 (MQ=255)
                 aaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCggg                                                                                                                                    <  3:103854/148‑1 (MQ=255)
                 aaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCggg                                                                                                                                    <  3:727631/148‑1 (MQ=255)
                     gCTGACGGTTTGTGATAATGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGtt                                                                                                                                >  3:874631/1‑148 (MQ=255)
                       tGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGc                                                                                                                             >  1:131430/1‑149 (MQ=255)
                        gACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCg                                                                                                                            >  1:411316/1‑149 (MQ=255)
                             tttGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATa                                                                                                                       >  1:10693/1‑149 (MQ=255)
                              ttGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATaa                                                                                                                      <  2:624435/149‑1 (MQ=255)
                              ttGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATa                                                                                                                       <  4:705912/148‑1 (MQ=255)
                               tgtgATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATaa                                                                                                                      <  4:240209/148‑1 (MQ=255)
                                gtgATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAAt                                                                                                                     <  4:320572/148‑1 (MQ=255)
                                 tgATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGc                                                                                                                   >  1:815364/1‑149 (MQ=255)
                                  gATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATAGCCGTTTCGCCGGGAGTTAGCGTCATAATGc                                                                                                                   <  3:316368/148‑1 (MQ=255)
                                     aaCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAt                                                                                                               <  1:137564/149‑1 (MQ=255)
                                      aCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATTGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAt                                                                                                               >  3:571284/1‑148 (MQ=255)
                                       cGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAtt                                                                                                              >  4:658770/1‑148 (MQ=255)
                                        gCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTg                                                                                                             >  4:276694/1‑148 (MQ=255)
                                            cttctACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGcc                                                                                                         <  4:303439/148‑1 (MQ=255)
                                                 aCTGCGCCCACCGAGACATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATc                                                                                                   <  2:910799/149‑1 (MQ=255)
                                                       cccACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAAc                                                                                             <  2:1015333/149‑1 (MQ=255)
                                                         cACCGAGCCATAAATCACCAATGCACCGCTGAATTGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACg                                                                                            <  4:571284/148‑1 (MQ=255)
                                                         cACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACg                                                                                            >  3:77176/1‑148 (MQ=255)
                                                             gAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGAt                                                                                       <  2:804445/149‑1 (MQ=255)
                                                              aGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGAt                                                                                       >  3:961521/1‑148 (MQ=255)
                                                                 cATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGAtttt                                                                                    <  4:392967/148‑1 (MQ=255)
                                                                  aTAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttc                                                                                   >  4:780145/1‑148 (MQ=255)
                                                                   tAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttctt                                                                                 <  2:10693/149‑1 (MQ=255)
                                                                    aaaTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttcttc                                                                                <  2:751156/149‑1 (MQ=255)
                                                                      aTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGc                                                                              <  1:437766/149‑1 (MQ=255)
                                                                       tCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGc                                                                              >  3:371223/1‑148 (MQ=255)
                                                                        cACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCa                                                                            <  2:362810/149‑1 (MQ=255)
                                                                          ccAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAg                                                                           >  4:842078/1‑148 (MQ=255)
                                                                           cAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGtt                                                                         >  1:161794/1‑149 (MQ=255)
                                                                           cAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGtt                                                                         >  1:450945/1‑149 (MQ=255)
                                                                            aaTGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTc                                                                        >  1:675033/1‑149 (MQ=255)
                                                                               gCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCt                                                                      <  3:248511/148‑1 (MQ=255)
                                                                                  ccGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCAc                                                                   >  4:365285/1‑148 (MQ=255)
                                                                                  ccGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCAc                                                                   >  4:801266/1‑148 (MQ=255)
                                                                                     cTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCggg                                                               >  1:689257/1‑149 (MQ=255)
                                                                                             gTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGcgcg                                                        <  3:751043/148‑1 (MQ=255)
                                                                                              tCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGa                                                       >  3:589148/1‑148 (MQ=255)
                                                                                              tCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGAt                                                      <  2:131430/149‑1 (MQ=255)
                                                                                              tCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGAt                                                      >  1:800184/1‑149 (MQ=255)
                                                                                              tCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGAt                                                      <  2:417364/149‑1 (MQ=255)
                                                                                               cGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGAt                                                      >  3:104262/1‑148 (MQ=255)
                                                                                                gagaAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTTGCCAGTTCCTCACCGGGATGCGCGATga                                                    >  2:810137/1‑149 (MQ=255)
                                                                                                   aaaaCCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATgaga                                                  >  4:929921/1‑148 (MQ=255)
                                                                                                     aaCCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATgc                                               <  2:685818/149‑1 (MQ=255)
                                                                                                             tGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAgtgt                                       <  1:917353/149‑1 (MQ=255)
                                                                                                             tGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAgtg                                        <  4:874631/148‑1 (MQ=255)
                                                                                                              gCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGtgtg                                      >  2:153919/1‑149 (MQ=255)
                                                                                                                     aGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGtt                               <  2:1077111/149‑1 (MQ=255)
                                                                                                                      gCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGttt                              >  2:994719/1‑149 (MQ=255)
                                                                                                                         gCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGc                            <  4:961521/148‑1 (MQ=255)
                                                                                                                            ttCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCgg                         >  4:915148/1‑148 (MQ=255)
                                                                                                                                 cgcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGcaca                    >  4:518261/1‑148 (MQ=255)
                                                                                                                                  gcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGcc                                                                                          >  4:391594/1‑77 (MQ=255)
                                                                                                                                  gcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGcc                                                                                          <  3:391594/77‑1 (MQ=255)
                                                                                                                                   cTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACaaa                  >  3:507512/1‑148 (MQ=255)
                                                                                                                                   cTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGACATGCGCCAGTGTGACCTGTTTGCCGGGCACaaa                  >  3:791405/1‑148 (MQ=255)
                                                                                                                                       aTCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTcc              <  3:51198/148‑1 (MQ=255)
                                                                                                                                       aTCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCAATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTcc              <  3:915103/148‑1 (MQ=255)
                                                                                                                                        tCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCt             <  3:932425/148‑1 (MQ=255)
                                                                                                                                         cACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGa           <  2:800184/149‑1 (MQ=255)
                                                                                                                                          aCCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGa           <  4:507512/148‑1 (MQ=255)
                                                                                                                                          aCCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGa           <  4:371223/148‑1 (MQ=255)
                                                                                                                                          aCCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGa           <  4:77176/148‑1 (MQ=255)
                                                                                                                                           acGGCGATGATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATgaga                                                  >  4:180714/2‑108 (MQ=255)
                                                                                                                                           acGGCGATGATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATgaga                                                  <  3:180714/107‑1 (MQ=255)
                                                                                                                                              gCGATCATCGCCGTTTCGCCGGGGGTTAGCGTCATAATGCTAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAatgat      >  2:699848/1‑149 (MQ=255)
                                                                                                                                                  tcatcGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGCg   >  3:643669/1‑148 (MQ=255)
                                                                                                                                                  tcatcGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGCg   <  3:780145/148‑1 (MQ=255)
                                                                                                                                                  tcatcGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGCGt  >  1:174760/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTTAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGCGT  >  NC_000913/2575500‑2575794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: