Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,269,032 G→T 11.5% intergenic (+13/+136) kdsA → / ← ldrA 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase/toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,269,0320GT11.5% 187.0 / 4.1 70intergenic (+13/+136)kdsA/ldrA3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base T (7/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCGGATCCGGAACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTA  >  NC_000913/1268886‑1269180
                                                                                                                                                  |                                                                                                                                                     
atccggatccggAACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGaa                                                                                                                                                     >  2:239849/1‑149 (MQ=255)
   cggatccggaACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGaaaaa                                                                                                                                                  >  2:268445/1‑149 (MQ=255)
     gatccggaACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGCACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATa                                                                                                                                                <  2:1224324/149‑1 (MQ=255)
                tGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTcgcg                                                                                                                                     >  1:595566/1‑149 (MQ=255)
                tGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCg                                                                                                                                     >  1:595343/1‑149 (MQ=255)
                 gCGAAAGGTGATGGTCCATCCGCGCGGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGACAGCGATTGTTGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAAGAAAGTATTAGCGt                                                                                                                                    >  1:525256/1‑149 (MQ=255)
                   gAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAACAAGTAATCTTTGTTGCTTGAAACAGAAAGTATTAgggggc                                                                                                                                  >  2:845507/1‑144 (MQ=255)
                    aaaTGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCt                                                                                                                                 <  2:208942/149‑1 (MQ=255)
                     aaTGTGGTGGTCCCTCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCTGAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTg                                                                                                                                <  1:1126496/149‑1 (MQ=255)
                       tgtgATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCGTCAAGAAGAAAAAAGCGATTGATGAACTGAAGGACGGTTTTGAGGAAATGGATTCCAGCAAGGAATTTTTTCTTCTAGAAAAATAAAGGATCATCGTTCTGCg                                                                                                                              >  2:788872/1‑149 (MQ=255)
                        gtgATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGGTGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAACGTATTAGCGTTCTGCGt                                                                                                                             <  2:998594/149‑1 (MQ=255)
                          gATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCGGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAGAAAAAAAGTATTAGCGTTTTGCATTa                                                                                                                           >  2:277845/1‑149 (MQ=255)
                              gTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGAc                                                                                                                       >  1:388279/1‑149 (MQ=255)
                                ccATCCGGGGTGCCGCTGGCTAACCTGGAACCGTTCCGCAAGCAGATGAAAGCGATTGATGAGCTGGTGAAAGGCTTCCAAGAACTTGATACCCGCAAGTAATGTTTTTTGCTTGAAAAACAAAGTATTTGCGTTCTGCGTTAAGACtt                                                                                                                     <  2:615152/149‑1 (MQ=255)
                                       ctcTGCCGCTGGCTAACCTGGAGCCGTTCCTCAAGCAGATGACAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCAt                                                                                                              <  2:435431/147‑1 (MQ=255)
                                       cgcTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGATCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCAt                                                                                                              <  2:324012/149‑1 (MQ=255)
                                            ccGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTTTAGTATAAAGTATTAGCGTTCTGCGGTAAGACTTTTTTCGtgtgtg                                                                                                         >  2:1104673/1‑149 (MQ=255)
                                              gCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGcc                                                                                                       >  1:23499/1‑149 (MQ=255)
                                               cTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCg                                                                                                      >  1:1204352/1‑149 (MQ=255)
                                               cTGGCTAAACTGGAACCGGTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATCCCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCg                                                                                                      <  2:781608/149‑1 (MQ=255)
                                                 ggCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTCTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGa                                                                                                   <  1:1055781/149‑1 (MQ=255)
                                                   cTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATa                                                                                                 <  1:409860/149‑1 (MQ=255)
                                                    tAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATAc                                                                                                >  1:1222425/1‑149 (MQ=255)
                                                    tAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATAc                                                                                                >  1:433699/1‑149 (MQ=255)
                                                    tAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATGATCGTTCTGCGTTAAGACTTTTTTCATGGGTGACGGATAc                                                                                                >  1:1223184/1‑149 (MQ=255)
                                                          tgAACCGGTCCTCAAGCCGATGAAAGCGATTGATGATCTGGTGAAAGGGTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACaaaaaa                                                                                          <  2:194032/148‑1 (MQ=255)
                                                               cGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCg                                                                                     <  1:784862/149‑1 (MQ=255)
                                                                     tCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCt                                                                               >  2:496600/1‑149 (MQ=255)
                                                                           agagGAAAGCGATTGATGACCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTAcc                                                                         <  1:416581/149‑1 (MQ=255)
                                                                           aGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAGCTGGATACCAGCAAGGATTCTTTTGTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTAcc                                                                         <  1:1118864/149‑1 (MQ=255)
                                                                           aGATGAAAGCGATTGATGATCTGGTGAAAGGTATCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTAcc                                                                         >  2:627299/1‑149 (MQ=255)
                                                                                 aaGCGATGGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCg                                                                   <  2:525256/149‑1 (MQ=255)
                                                                                      aTTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAACTAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAATGGCCGCAGACTGTTACGCCTGCGGCCgg                                                              >  2:1346279/1‑149 (MQ=255)
                                                                                      aTTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCgg                                                              >  2:1346261/1‑149 (MQ=255)
                                                                                          atgatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATTGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTc                                                          >  2:462254/1‑149 (MQ=255)
                                                                                            gatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCgg                                                        >  2:965091/1‑149 (MQ=255)
                                                                                            gatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCAACGGAAAGACTTTTTTATGAGTCAAGAAACTCTAAAAAGGGAGCAGTCACTTACTCCTGCGGCCGGTTTCgg                                                        >  1:642458/1‑149 (MQ=255)
                                                                                             gacTGGGGAAAGTTTTACAAACACCGGCTACCAGCCAGTAATCTCTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                                                       <  2:229358/147‑1 (MQ=37)
                                                                                             atCTGGTGAAAGGTTTCGAAGAACCGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                                                       >  1:1004466/1‑149 (MQ=255)
                                                                                             atCTAGTGAAAGTTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                                                       <  2:140601/149‑1 (MQ=255)
                                                                                                 ggTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCa                                                   >  2:189352/1‑149 (MQ=255)
                                                                                                   tGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCGCCGGCCGGTTGCGGGCGCata                                                 >  1:368748/1‑149 (MQ=255)
                                                                                                       aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTTTTTTTATTTTAAAGTATTAGCGGTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGc                                             >  2:454254/1‑149 (MQ=37)
                                                                                                       aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTAGCCCTGCGGCCGGTTTCGGGCGCATATTGc                                             <  1:836528/149‑1 (MQ=255)
                                                                                                       aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGc                                             <  2:97614/149‑1 (MQ=255)
                                                                                                         gTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTTTTTTTAAAATAAAGTATTAGCGTTCCGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGGGGCCGGTTTCGGGCGCATATTGCCa                                           >  1:732202/1‑149 (MQ=38)
                                                                                                           ttCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGCATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATc                                         <  2:388279/149‑1 (MQ=255)
                                                                                                               aagaaCTGGATACCAGCAAGTAATCTTTTTTTCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGACCGGACTCGAGCGCATATTGCCGTCGCgg                                     >  2:715853/1‑149 (MQ=38)
                                                                                                               aagaaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACgg                                     <  1:32016/149‑1 (MQ=255)
                                                                                                                    cTGGATACCAGCAAGTAATCGTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCCGATACAAACAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGcc                                <  1:289577/149‑1 (MQ=255)
                                                                                                                     tGGATACCAGCAAGTAATCTTTTTTTTTTTTTTTTTAAAGTATTAGCGGTCTGCGTTAAGACTTTTTTCATGGGTGGCGGATACAAAAAAGGCCTCCGGCTTTTACCGCTGCGGCCGGTTTCGGGCGCTTTTTTCCATTACCGCAGctt                               >  2:276577/1‑147 (MQ=19)
                                                                                                                            cAGCAAGTAATTTTTTTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGccc                        <  1:496600/149‑1 (MQ=255)
                                                                                                                            cAGCAAGTAATCTTTTTTTTTTTTTTTTTAAAGTATTAACGTCCTGCGTTAAGACTTTTTTCATGGGGGCCGAATACAGGAAAGGCCGCAGGCTGTTACCCCTGAGGCCGTTTTCGTTCGCTTATTGCCATAACGGCCGCCTGACGccc                        >  2:674642/1‑149 (MQ=16)
                                                                                                                            cAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCTGCCGGAATCGGGCGCAAATTGCCATCACGGCAGCCTGACGccc                        >  1:849867/1‑149 (MQ=255)
                                                                                                                             aGCCAGTAGTCTTTTTCGCCTGAAAAATAAAGTATTAGCGGTCTGCGTCAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCCGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCg                       <  2:150841/149‑1 (MQ=25)
                                                                                                                             aGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGCCGCCCg                       >  2:65094/1‑149 (MQ=255)
                                                                                                                               cAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGtt                     >  1:1075693/1‑149 (MQ=255)
                                                                                                                                   tAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCa                 >  1:360667/1‑149 (MQ=255)
                                                                                                                                   tAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCATCCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCa                 <  1:1273980/149‑1 (MQ=255)
                                                                                                                                    aaTCTTTTTTGCTTGAAAAATAAAGTATTAGCGGTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCAc                <  1:239849/149‑1 (MQ=255)
                                                                                                                                      tCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTAATGGGTGGCGTATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCCCGGCCGCCTGACGCCCGTTTTCACCt              >  2:622692/1‑149 (MQ=25)
                                                                                                                                        ttttttTGTTGAAAAATAAAGTATTGGCTGTCTGCGTTAAGATCTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGTATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTa            <  2:60983/149‑1 (MQ=18)
                                                                                                                                        ttttttGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTa            <  1:1165482/149‑1 (MQ=255)
                                                                                                                                        tttgttGCTTGAAAAATAAAGTATTAGCGTTATGCGGTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCGGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTa            <  2:732154/149‑1 (MQ=255)
                                                                                                                                         tttttGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTAc           >  1:1186947/1‑149 (MQ=255)
                                                                                                                                              gCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCg      <  1:189352/149‑1 (MQ=255)
                                                                                                                                               cTTGAATCATAACGTAGTTGCGTTCTGTCGTATGCGTCTTCTCATGGCTGCCGGATACAAAAAAGGCCGCAGGCTGGTACCCCGTCCGCCGGTTTCGCGCGCCTATTGCCATCACGTCAGCCTGACGCCCGTTTTCCCCTTACTTCCgg     <  2:1335417/149‑1 (MQ=255)
                                                                                                                                               cTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCgg     <  2:1204352/149‑1 (MQ=255)
                                                                                                                                               cTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCgg     <  1:627299/149‑1 (MQ=255)
                                                                                                                                                 tGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGtt   >  2:344832/1‑149 (MQ=255)
                                                                                                                                                  gAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTa  <  2:595343/149‑1 (MQ=255)
                                                                                                                                                  gAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTa  <  2:23499/149‑1 (MQ=255)
                                                                                                                                                  gAAAAATAAAGTATTAGCGTGCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTa  <  2:595566/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCCGGATCCGGAACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTA  >  NC_000913/1268886‑1269180

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