Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,070,882 T→A 6.9% E25V (GAG→GTG)  rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,8820TA 6.9% 203.0 / 5.2 72E25V (GAG→GTG) rutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (2/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

TTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTT  >  NC_000913/1070735‑1071031
                                                                                                                                                   |                                                                                                                                                     
ttGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGACCAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCAGACGGGTGTCTcc                                                                                                                                                     >  2:537456/1‑150 (MQ=255)
  gCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACCCGATCCGTGCCGGAGAACAGTTTGCTGATGTCGGCCCCCAGTTCATCAGGGCATAAATCGCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCg                                                                                                                                                   >  1:821271/1‑150 (MQ=255)
     gCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGAGGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCa                                                                                                                                                <  1:410125/150‑1 (MQ=255)
       gggAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCGGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGGCTCATCGCTGGCGGGAGACACCCGTCAGc                                                                                                                                              <  1:215772/149‑1 (MQ=255)
        ggAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGcc                                                                                                                                             <  2:1629853/150‑1 (MQ=255)
         gAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAGcga                                                                                                                                            <  1:1352766/150‑3 (MQ=255)
          aaaGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGGCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGt                                                                                                                                           <  1:657827/150‑1 (MQ=255)
           aaGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAGcgatg                                                                                                                                          <  1:1164305/150‑5 (MQ=255)
            aGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTtg                                                                                                                                         <  2:587541/150‑1 (MQ=255)
              tGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGGtgtg                                                                                                                                       >  1:1077170/1‑150 (MQ=255)
               gCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGa                                                                                                                                      >  1:971490/1‑150 (MQ=255)
                ccGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAg                                                                                                                                     <  2:1150911/150‑1 (MQ=255)
                   ggCGCAGACGGTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGACCAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGc                                                                                                                                  <  1:1386837/150‑1 (MQ=255)
                     cgcAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGGCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCgg                                                                                                                                <  1:1388018/150‑1 (MQ=255)
                           cGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCg                                                                                                                          >  2:1205929/1‑150 (MQ=255)
                           cGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCACCCGCCAGCCGTTGTGAGTGCGGGCATCg                                                                                                                          >  1:435307/1‑150 (MQ=255)
                             ttCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGt                                                                                                                        <  1:580312/150‑1 (MQ=255)
                              tCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTg                                                                                                                       <  2:1471293/150‑1 (MQ=255)
                                     cTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAggg                                                                                                                <  2:69979/150‑1 (MQ=255)
                                      ttCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGt                                                                                                               >  1:616138/1‑150 (MQ=255)
                                        cTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGc                                                                                                             >  1:1410187/1‑150 (MQ=255)
                                         tGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTTTGAGTGCGGGCATCGGTGAACAGGGTGCt                                                                                                            <  1:1555128/150‑1 (MQ=255)
                                         tGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCt                                                                                                            <  2:479384/150‑1 (MQ=255)
                                         tGCCGTGCGGGT‑‑ACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAACTCTCCCGTAACGTCTCAACGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTaa                                                                                                          >  1:1338378/1‑150 (MQ=255)
                                            cGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCGGGTGTGCGGGCGGGCATCGGTGAACCGGGTggttag                                                                                                         >  2:166495/1‑145 (MQ=255)
                                            cGTGAGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAg                                                                                                         <  2:1554427/150‑1 (MQ=255)
                                                  ggTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCAcc                                                                                                   >  1:802420/1‑150 (MQ=255)
                                                    tAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTg                                                                                                 <  2:1455803/150‑1 (MQ=255)
                                                     aaaCACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTgg                                                                                                >  2:1538041/1‑150 (MQ=255)
                                                     aaaCACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTgg                                                                                                >  2:1288592/1‑150 (MQ=255)
                                                     aaaCACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTgg                                                                                                >  2:759444/1‑150 (MQ=255)
                                                     aaaCACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTgg                                                                                                >  2:557662/1‑150 (MQ=255)
                                                     aaaCACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCAgg                                                                                                >  2:1617797/1‑150 (MQ=255)
                                                        cacGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCt                                                                                             <  1:1270987/150‑1 (MQ=255)
                                                        cacGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCACCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCt                                                                                             <  1:1415570/150‑1 (MQ=255)
                                                        cacGATCCGTGCCGGAGAACAGTTAGCTGAAGTATGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCACATCGCTGACGGGTGTCTCCAGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCt                                                                                             >  1:1163247/1‑150 (MQ=255)
                                                                  gCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCg                                                                                   <  2:455817/150‑1 (MQ=255)
                                                                      gagaACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCa                                                                               <  1:256545/150‑1 (MQ=255)
                                                                       agaACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCAt                                                                              >  1:982430/1‑150 (MQ=255)
                                                                       agaACAGTTAGCTGACGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCAt                                                                              >  1:529578/1‑150 (MQ=255)
                                                                         aaCAGTTCGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGGCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAg                                                                            <  1:943945/150‑1 (MQ=255)
                                                                         aaCAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCGCATCGCTGACGGGTGTCGCCCGCCGGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAg                                                                            <  2:1350437/150‑1 (MQ=255)
                                                                                  cTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCtt                                                                   <  1:557662/150‑1 (MQ=255)
                                                                                  cTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCtt                                                                   <  1:1538041/150‑1 (MQ=255)
                                                                                  cTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCtt                                                                   >  1:443302/1‑150 (MQ=255)
                                                                                  cTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGGCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTTAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCtt                                                                   <  2:653213/150‑1 (MQ=255)
                                                                                    gATGTCGGCCCCCCTTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGAGGTCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTAc                                                                 <  1:888255/150‑1 (MQ=255)
                                                                                       gTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAgg                                                              >  1:464981/1‑150 (MQ=255)
                                                                                               ccATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTc                                                      <  2:560218/150‑1 (MQ=255)
                                                                                                  tttCTTCAGGGCATAAATGTGCCGTATGGTCTCATCGCTGACGGGGGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTTAACAGGGTGCTAATCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACg                                                   <  1:1155966/150‑1 (MQ=255)
                                                                                                  tttCATCAGGGCATAATTCTGCAGTTACGTCTCATCGCTGCCTGGTGTCTCACGCCAGACGTTGTGAGTGCGGGCATCTGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACg                                                   <  1:36845/150‑1 (MQ=255)
                                                                                                     catcaGGGCATAATTCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGGGTTTGTGAGTGCGGGAATCGTTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATTGCAATTCCTTACAGGGCGGCTTCACGGTg                                                <  2:1163247/150‑1 (MQ=255)
                                                                                                      atccGGGCATAAATCTCCCGTATCGTCTCAGCGCTGACGGGTGTCTCCCGCCAGACGTTGTGAGTGCGGGCATCGGGGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGa                                               <  2:1064109/150‑1 (MQ=255)
                                                                                                        caGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATg                                             <  2:1077170/150‑1 (MQ=255)
                                                                                                             cATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACGAGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCGTCACGGTGATGTAACa                                        >  2:128215/1‑150 (MQ=255)
                                                                                                                aaaTCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGc                                     >  2:1212688/1‑150 (MQ=255)
                                                                                                                   tctcCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTgg                                  <  2:1267655/150‑1 (MQ=255)
                                                                                                                   tctcCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTgg                                  >  1:1544636/1‑150 (MQ=255)
                                                                                                                    ctcCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGc                                 >  2:460584/1‑150 (MQ=255)
                                                                                                                    ctcCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGc                                 >  2:7309/1‑150 (MQ=255)
                                                                                                                         gTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGcc                            <  2:529578/150‑1 (MQ=255)
                                                                                                                         ccgcactctcATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGcc                            >  1:744729/7‑150 (MQ=255)
                                                                                                                           aaCGTCGCATCGCTGACGGGTGGCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCg                          <  2:268080/150‑1 (MQ=255)
                                                                                                                             cGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTACTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGtt                        >  2:1519356/1‑150 (MQ=255)
                                                                                                                                ctcATCGCTGAGCGGGGTCTCCCGCCAGCCGGTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGAGGTATCAGGCTGGCAAGCCCGTTGAg                     <  1:482829/150‑1 (MQ=255)
                                                                                                                                ctcATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAg                     >  2:987408/1‑150 (MQ=255)
                                                                                                                                  cATCGCTGACGGGTGTCTCCCGCCTGCCGTTGTGAGTGCGGGCATCGGTTAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTa                   <  1:1617797/150‑1 (MQ=255)
                                                                                                                                      gCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACag               >  1:517480/1‑150 (MQ=255)
                                                                                                                                       cTGGCGGGAGACACCCGTCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACaga              >  1:621912/1‑150 (MQ=255)
                                                                                                                                           cGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGCACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCAt          >  1:1162418/1‑150 (MQ=255)
                                                                                                                                             ggTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTa        >  1:1620135/1‑150 (MQ=255)
                                                                                                                                             ggTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTa        >  1:1522579/1‑150 (MQ=255)
                                                                                                                                                 tctcCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATTGCAATTCCGTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACg    >  1:1210886/1‑150 (MQ=255)
                                                                                                                                                  ctcCCGCCAGCCGTTGTGAGTGTGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCCCGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGtt  >  2:575239/1‑150 (MQ=255)
                                                                                                                                                   tcCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGtt  >  1:1724222/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTT  >  NC_000913/1070735‑1071031

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: