Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,575,647 C→G 8.5% 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 8.5% 330.0 / 18.5 94M53I (ATG→ATCeutSputative ethanol utilization carboxysome structural protein
Reads supporting (aligned to +/- strand):  ref base C (46/40);  new base G (5/3);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

TTTAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGC  >  NC_000913/2575501‑2575792
                                                                                                                                                  |                                                                                                                                                 
tttAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGatcat                                                                                                                                                 <  2:554184/149‑1 (MQ=255)
          gACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCg                                                                                                                                       >  1:311951/1‑149 (MQ=255)
                aaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 <  2:1013040/149‑1 (MQ=255)
                aaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 >  2:878897/1‑149 (MQ=255)
                aaCCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 <  4:1284124/149‑1 (MQ=255)
                  ccGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggtga                                                                                                                                        <  3:1100422/140‑5 (MQ=255)
                  ccGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggtga                                                                                                                                        >  4:1100422/1‑136 (MQ=255)
                    gCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCAAAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTa                                                                                                                             <  4:650729/149‑1 (MQ=255)
                     cTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAg                                                                                                                            <  1:111662/149‑1 (MQ=255)
                          ggTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCa                                                                                                                       >  3:760391/1‑149 (MQ=255)
                            tttGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATa                                                                                                                     >  4:611783/1‑149 (MQ=255)
                            tttGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATa                                                                                                                     >  4:1084988/1‑149 (MQ=255)
                                 gATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTa                                                                                                                             <  4:319000/136‑1 (MQ=255)
                                 gATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTa                                                                                                                             >  3:319000/1‑136 (MQ=255)
                                   tAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCaa                                                                                                              <  3:1084988/149‑1 (MQ=255)
                                    aaCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggt                                                                                                                                          >  1:1234322/1‑118 (MQ=255)
                                    aaCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggt                                                                                                                                          <  2:1234322/120‑3 (MQ=255)
                                    aaCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 >  1:1298568/1‑129 (MQ=255)
                                    aaCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 <  2:1298568/129‑1 (MQ=255)
                                     aCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAtt                                                                                                            <  3:611783/149‑1 (MQ=255)
                                     aCGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAtt                                                                                                            >  4:1080071/1‑149 (MQ=255)
                                       gCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTgc                                                                                                          >  2:1031065/1‑149 (MQ=255)
                                        cttcttCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTgcg                                                                                                         <  4:760391/149‑1 (MQ=255)
                                            ttctACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggtg                                                                                                                                         >  2:1124389/1‑110 (MQ=255)
                                            ttctACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATGATCGCCggtg                                                                                                                                         <  1:1124389/113‑4 (MQ=255)
                                               tACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 <  2:1038986/118‑1 (MQ=255)
                                               tACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 >  1:1038986/1‑118 (MQ=255)
                                                   gcgcCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCgg                                                                                              >  1:683884/1‑149 (MQ=255)
                                                      cccACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAAc                                                                                           >  2:1299260/1‑149 (MQ=255)
                                                       ccACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGaa                                                                                            >  3:279906/1‑147 (MQ=255)
                                                       ccACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGaa                                                                                            <  4:279906/147‑1 (MQ=255)
                                                         aCCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGaa                                                                                            >  1:1123209/1‑145 (MQ=255)
                                                         aCCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGaa                                                                                            <  2:1123209/145‑1 (MQ=255)
                                                           cGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTc                                                                                                                                        >  1:1114043/1‑99 (MQ=255)
                                                           cGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGa                                                                                      <  1:878897/149‑1 (MQ=255)
                                                           cGAGCCATAAAACACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTc                                                                                                                                        <  2:1114043/99‑1 (MQ=255)
                                                               ccATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 <  4:1410850/102‑1 (MQ=255)
                                                               ccATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 >  3:1410850/1‑102 (MQ=255)
                                                                cATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttc                                                                                 >  1:669196/1‑149 (MQ=255)
                                                                cATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttc                                                                                 <  1:1299260/149‑1 (MQ=255)
                                                                cATAAATCACCAATGCACCGCTGAATAGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCTCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttc                                                                                 >  2:131463/1‑149 (MQ=255)
                                                                 aTAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 <  4:182582/100‑1 (MQ=255)
                                                                 aTAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGa                                                                                                                                 >  3:182582/1‑100 (MQ=255)
                                                                  tAAATCACCAATGCACCGCTGAATCGGTCGAGAACACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGCCATAATGCCAATTGCGCCCGCATCCGGAACGCCGa                                                                                      <  4:361513/142‑1 (MQ=255)
                                                                  tAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGa                                                                                      >  3:361513/1‑142 (MQ=255)
                                                                  tAAATCACCAATGCACCGCTGAATAGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttctt                                                                               >  1:867506/1‑149 (MQ=255)
                                                                   aaaTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttcttc                                                                              <  3:1080071/149‑1 (MQ=255)
                                                                     aTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttcttc                                                                              <  3:793628/147‑1 (MQ=255)
                                                                     aTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttcttc                                                                              >  4:793628/1‑147 (MQ=255)
                                                                     aTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGc                                                                            >  3:860752/1‑149 (MQ=255)
                                                                     aTCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGc                                                                            <  2:311951/149‑1 (MQ=255)
                                                                      tCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGcc                                                                           >  2:572127/1‑149 (MQ=255)
                                                                      tCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGcc                                                                           <  1:1031065/149‑1 (MQ=255)
                                                                        aCCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAg                                                                         >  3:586822/1‑149 (MQ=255)
                                                                           aaTGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTc                                                                      <  4:860752/149‑1 (MQ=255)
                                                                               cACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAt                                                                                                             <  1:546146/106‑1 (MQ=255)
                                                                               cACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAAt                                                                                                             >  2:546146/1‑106 (MQ=255)
                                                                                aCCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATaa                                                                                                                    >  3:1032297/1‑98 (MQ=255)
                                                                                aCCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATaa                                                                                                                    <  4:1032297/98‑1 (MQ=255)
                                                                                       aaTCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGtt                                                                                                                              <  4:1414595/81‑1 (MQ=255)
                                                                                       aaTCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGtt                                                                                                                              >  3:1414595/1‑81 (MQ=255)
                                                                                        aTCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATgc                                                         <  2:669196/149‑1 (MQ=255)
                                                                                          cGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATgcgc                                                       <  1:572127/149‑1 (MQ=255)
                                                                                                 gaAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATgaga                                                >  3:847580/1‑149 (MQ=255)
                                                                                                    aaCCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATgc                                             >  2:57175/1‑149 (MQ=255)
                                                                                                     aCCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATgcg                                            >  3:293990/1‑149 (MQ=255)
                                                                                                        gATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCa                                         <  2:683884/149‑1 (MQ=255)
                                                                                                               acaTCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGAc                                  >  3:1168218/1‑149 (MQ=255)
                                                                                                                  tCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTg                               >  1:1178639/1‑149 (MQ=255)
                                                                                                                     gCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGc                                                                                                                           <  3:343927/54‑1 (MQ=255)
                                                                                                                     gCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGc                                                                                                                           >  4:343927/1‑54 (MQ=255)
                                                                                                                           ttCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCggg                      >  1:972462/1‑149 (MQ=255)
                                                                                                                               gcgcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGc                                                                            >  2:454483/1‑91 (MQ=255)
                                                                                                                               gcgcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGc                                                                            <  1:454483/91‑1 (MQ=255)
                                                                                                                               gcgcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGcaca                  >  4:708753/1‑149 (MQ=255)
                                                                                                                                 gcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACaaa                >  1:502516/1‑149 (MQ=255)
                                                                                                                                 gcTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACaaa                >  1:1395500/1‑149 (MQ=255)
                                                                                                                                  cTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAAt               <  2:502516/149‑1 (MQ=255)
                                                                                                                                  cTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAAt               <  2:867506/149‑1 (MQ=255)
                                                                                                                                    aaaTCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTc             <  4:586822/149‑1 (MQ=255)
                                                                                                                                     aaTCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTcc            >  3:1205648/1‑149 (MQ=255)
                                                                                                                                      aTCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCt           >  3:1341960/1‑149 (MQ=255)
                                                                                                                                      aTCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCt           <  4:293990/149‑1 (MQ=255)
                                                                                                                                       tCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTg          >  3:635687/1‑149 (MQ=255)
                                                                                                                                        cACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttcttc                                                                              >  1:931634/1‑80 (MQ=255)
                                                                                                                                        cACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTttcttc                                                                              <  2:931634/80‑1 (MQ=255)
                                                                                                                                        cACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTcc            >  2:1106631/1‑146 (MQ=255)
                                                                                                                                        cACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTcc            <  1:1106631/146‑1 (MQ=255)
                                                                                                                                         aCCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGaa        <  3:708753/149‑1 (MQ=255)
                                                                                                                                          acGGCGATGATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAat       >  3:982652/2‑149 (MQ=255)
                                                                                                                                          acGGCGATGATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAat       >  1:284626/2‑149 (MQ=255)
                                                                                                                                             gCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAatgat    >  1:828679/1‑149 (MQ=255)
                                                                                                                                             gCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAatgat    <  2:828679/149‑1 (MQ=255)
                                                                                                                                             gCGATCATCGCAGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAatgat    >  1:1085982/1‑149 (MQ=255)
                                                                                                                                               gATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGc  <  4:1168218/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTTAATAAGCGACCTAAACCGCTGACGGTTTGTGATAACGCTTCTTCTACTGCGCCCACCGAGCCATAAATCACCAATGCACCGCTGAATCGGTCGAGAAAACCGATATGCACATCAGCGGCTTTCAGCGCTAAATCACCGGCGATCATCGCCGTTTCGCCGGGAGTTAGCGTCATAATGCCAATTGCGCCCGCATCCGGAACGCCGATTTTCTTCGCCAGTTCCTCACCGGGATGCGCGATGAGATGCGCCAGTGTGACCTGTTTGCCGGGCACAAATTCCTGAATGATGC  >  NC_000913/2575501‑2575792

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