Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,204,390 T→C 5.7% intergenic (+101/‑206) yehK → / → yehL uncharacterized protein/putative hexameric AAA+ MoxR family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,204,3900TC 5.7% 214.4 / 3.3 70intergenic (+101/‑206)yehK/yehLuncharacterized protein/putative hexameric AAA+ MoxR family ATPase
Reads supporting (aligned to +/- strand):  ref base T (40/26);  new base C (2/2);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

TTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGC  >  NC_000913/2204246‑2204535
                                                                                                                                                |                                                                                                                                                 
ttAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTtata                                                                                                                                                <  2:217988/148‑1 (MQ=255)
 tAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATaa                                                                                                                                               <  2:253405/148‑1 (MQ=255)
         aaCTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGTGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATtcgttcgt                                                                                                                                      >  4:232923/1‑149 (MQ=255)
               aGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGaa                                                                                                                                 >  2:261684/1‑148 (MQ=255)
                       ttCAACCAAGAACTAGTATAATATACTCTAAATAATTCAATTTGGTCAGCTCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCAGTGGTTATAAttcgttc                                                                                                                                        <  2:294864/133‑1 (MQ=255)
                         cAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGGCCCAATGGCCTGGATTATAACTCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTg                                                                                                                       <  7:324780/148‑1 (MQ=255)
                          aaCCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGca                                                                                                                     <  1:237490/149‑1 (MQ=255)
                            ccAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGca                                                                                                                     >  7:180706/1‑147 (MQ=255)
                               aGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGc                                                                                                                <  4:305536/149‑1 (MQ=255)
                               aGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCAGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGc                                                                                                                >  1:232435/1‑149 (MQ=255)
                                   aTAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCtgtg                                                                                                            <  1:261684/149‑1 (MQ=255)
                                      gAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTgtggtg                                                                                                         >  6:257847/1‑149 (MQ=255)
                                      gAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTgtggtg                                                                                                         >  1:355829/1‑149 (MQ=255)
                                       aataatATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGgtgt                                                                                                        <  5:257847/149‑1 (MQ=255)
                                         taataTACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGTCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTa                                                                                                      <  4:25292/149‑1 (MQ=255)
                                            tataCTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGc                                                                                                   >  2:145803/1‑149 (MQ=255)
                                                 tctAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGa                                                                                               >  6:254033/1‑148 (MQ=255)
                                                     aataatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGaaaa                                                                                            <  7:191819/147‑1 (MQ=255)
                                                     aataatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCGGaaaa                                                                                            >  8:111865/1‑147 (MQ=255)
                                                       taatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGt                                                                                        <  8:267145/149‑1 (MQ=255)
                                                        aatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGt                                                                                        >  1:247891/1‑148 (MQ=255)
                                                             aaaTTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGtaaata                                                                                   >  2:45913/1‑148 (MQ=255)
                                                              aaTTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGtaaataa                                                                                  >  7:114449/1‑148 (MQ=255)
                                                                   gTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAAc                                                                            <  1:145803/149‑1 (MQ=255)
                                                                      cGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCCCAGGCTGTGGTGTTAGGCCTAAAAATAGTAAATAAAAGAACCTg                                                                         >  5:259885/1‑149 (MQ=255)
                                                                      cGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTg                                                                         <  8:180706/149‑1 (MQ=255)
                                                                       gATCCCGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAAAAGTAAATAAAAGAACCTga                                                                        >  2:177212/1‑148 (MQ=255)
                                                                             gcgcAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCAc                                                                   >  8:43904/1‑148 (MQ=255)
                                                                                 aaCGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACcag                                                                >  3:252525/1‑147 (MQ=255)
                                                                                 aaCGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTCATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACcag                                                                >  1:305740/1‑147 (MQ=255)
                                                                                 aaCGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTCATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACcag                                                                >  7:261548/1‑147 (MQ=255)
                                                                                    gTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGt                                                           <  6:71586/149‑1 (MQ=255)
                                                                                    gTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGt                                                           >  8:155308/1‑149 (MQ=255)
                                                                                       ccAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGtttt                                                        >  3:19311/1‑149 (MQ=255)
                                                                                         aaTGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGtttttt                                                      >  4:127484/1‑149 (MQ=255)
                                                                                             gCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGTAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttatt                                                  <  6:237591/149‑1 (MQ=255)
                                                                                             gCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttat                                                   <  3:232923/148‑1 (MQ=255)
                                                                                             gCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTAGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttatt                                                  >  3:240301/1‑149 (MQ=255)
                                                                                               cTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTCATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttattta                                                <  2:305740/149‑1 (MQ=255)
                                                                                               cTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTCATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttattta                                                <  8:261548/149‑1 (MQ=255)
                                                                                               cTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCAGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttattta                                                >  2:222377/1‑149 (MQ=255)
                                                                                                tGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTAc                                               >  7:130191/1‑149 (MQ=255)
                                                                                                  gATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACt                                              <  2:304827/148‑1 (MQ=255)
                                                                                                  gATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACt                                              >  1:304827/1‑148 (MQ=255)
                                                                                                   aTTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATTATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACtgt                                            >  1:83101/1‑149 (MQ=255)
                                                                                                       taAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGAt                                         <  7:43446/148‑1 (MQ=255)
                                                                                                       taAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGAt                                         <  7:111865/148‑1 (MQ=255)
                                                                                                        aaaTCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCt                                       >  1:371672/1‑149 (MQ=255)
                                                                                                         aaTCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCt                                       <  4:252525/148‑1 (MQ=255)
                                                                                                          aTCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGc                                     >  8:132467/1‑149 (MQ=255)
                                                                                                          aTCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGc                                     >  5:166086/1‑149 (MQ=255)
                                                                                                           tctcATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCt                                    >  5:9508/1‑149 (MQ=255)
                                                                                                           tctcATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGACTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCt                                    >  7:150951/1‑149 (MQ=255)
                                                                                                              cATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTTGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTc                                 >  5:237004/1‑149 (MQ=255)
                                                                                                               aTTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTc                                 >  4:159139/1‑148 (MQ=255)
                                                                                                                ttATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCATTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTc                                 >  2:200394/1‑147 (MQ=255)
                                                                                                                  aTCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAg                              <  1:45913/148‑1 (MQ=255)
                                                                                                                    cTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTc                                 <  7:192222/143‑1 (MQ=255)
                                                                                                                    cTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTc                                 >  8:192222/1‑143 (MQ=255)
                                                                                                                        aTTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATAttt                        <  2:355829/148‑1 (MQ=255)
                                                                                                                         ttGCAACTGGGTCAAGCCGTGGTTATAATTCGATCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTATAATATACTGTGATCTGCTTTCCAGATAttttt                      >  8:212735/1‑149 (MQ=255)
                                                                                                                         ttGCAACGGGGGCCAGGCGTGTTTATTATTTTTTCGTTGGAACACGATGTTGCACAGGCTGTTGTGTTAGGCCTTAAAATAGTAAATAAAAGAACCTGCGTCACCATGAGGTTTTTTTATTTACTGTGATGTGCTTTCCAGATAttttt                      <  4:240301/149‑1 (MQ=255)
                                                                                                                           gCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCg                    >  5:131418/1‑149 (MQ=255)
                                                                                                                             aaCGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGc                   >  8:86578/1‑148 (MQ=255)
                                                                                                                             aaCGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCt                  >  2:119761/1‑149 (MQ=255)
                                                                                                                                ggggTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCAc                                                                   >  2:348508/1‑97 (MQ=255)
                                                                                                                                ggggTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCAc                                                                   <  1:348508/97‑1 (MQ=255)
                                                                                                                                    tCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAaacaa           <  7:155308/149‑1 (MQ=255)
                                                                                                                                         ccGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAAt      >  2:268422/1‑149 (MQ=255)
                                                                                                                                           gTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgc    <  1:268422/149‑1 (MQ=255)
                                                                                                                                              gTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCGGGTTTCCAGATATTTTTCGGTCAAACAACTAATgcgc  >  2:120133/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGC  >  NC_000913/2204246‑2204535

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