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

ATCCGGATCCGGAACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTA  >  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)
                                                 ggCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTCTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑Ga                                                                                                   <  1:1055781/149‑1 (MQ=255)
                                                   cTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATa                                                                                                 <  1:409860/149‑1 (MQ=255)
                                                    tAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATAc                                                                                                >  1:1222425/1‑149 (MQ=255)
                                                    tAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATAc                                                                                                >  1:433699/1‑149 (MQ=255)
                                                    tAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATGATCGTTCTGCGTTAAGACTTTTTTCATGGGTGACG‑GATAc                                                                                                >  1:1223184/1‑149 (MQ=255)
                                                          tgAACCGGTCCTCAAGCCGATGAAAGCGATTGATGATCTGGTGAAAGGGTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACaaaaaa                                                                                          <  2:194032/148‑1 (MQ=255)
                                                               cGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCg                                                                                     <  1:784862/149‑1 (MQ=255)
                                                                     tCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCt                                                                               >  2:496600/1‑149 (MQ=255)
                                                                           agagGAAAGCGATTGATGACCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTAcc                                                                         <  1:416581/149‑1 (MQ=255)
                                                                           aGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAGCTGGATACCAGCAAGGATTCTTTTGTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTAcc                                                                         <  1:1118864/149‑1 (MQ=255)
                                                                           aGATGAAAGCGATTGATGATCTGGTGAAAGGTATCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTAcc                                                                         >  2:627299/1‑149 (MQ=255)
                                                                                 aaGCGATGGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCg                                                                   <  2:525256/149‑1 (MQ=255)
                                                                                      aTTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAACTAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAATGGCCGCAGACTGTTACGCCTGCGGCCgg                                                              >  2:1346279/1‑149 (MQ=255)
                                                                                      aTTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCgg                                                              >  2:1346261/1‑149 (MQ=255)
                                                                                          atgatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATTGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTc                                                          >  2:462254/1‑149 (MQ=255)
                                                                                            gatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCgg                                                        >  2:965091/1‑149 (MQ=255)
                                                                                            gatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCAACGGAAAGACTTTTTT‑ATGAGTCAAGAAACTCTAAAAAGGGAGCAGTCACTTACTCCTGCGGCCGGTTTCgg                                                        >  1:642458/1‑149 (MQ=255)
                                                                                             gacTGGGGAAAGTTTTACAAACACCGGCTACCAGCCAGTAATCTCTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                                                       <  2:229358/147‑1 (MQ=37)
                                                                                             atCTGGTGAAAGGTTTCGAAGAACCGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                                                       >  1:1004466/1‑149 (MQ=255)
                                                                                             atCTAGTGAAAGTTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                                                       <  2:140601/149‑1 (MQ=255)
                                                                                                 ggTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCa                                                   >  2:189352/1‑149 (MQ=255)
                                                                                                   tGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCGCCGGCCGGTTGCGGGCGCata                                                 >  1:368748/1‑149 (MQ=255)
                                                                                                       aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTTTTTTTATTTTAAAGTATTAGCGGTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGc                                             >  2:454254/1‑149 (MQ=37)
                                                                                                       aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTAGCCCTGCGGCCGGTTTCGGGCGCATATTGc                                             <  1:836528/149‑1 (MQ=255)
                                                                                                       aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGc                                             <  2:97614/149‑1 (MQ=255)
                                                                                                         gTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTTTTTTTAAAATAAAGTATTAGCGTTCCGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGGGGCCGGTTTCGGGCGCATATTGCCa                                           >  1:732202/1‑149 (MQ=38)
                                                                                                           ttCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGCATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATc                                         <  2:388279/149‑1 (MQ=255)
                                                                                                               aagaaCTGGATACCAGCAAGTAATCTTTTTTTCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGACCGGACTCGAGCGCATATTGCCGTCGCgg                                     >  2:715853/1‑149 (MQ=38)
                                                                                                               aagaaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACgg                                     <  1:32016/149‑1 (MQ=255)
                                                                                                                    cTGGATACCAGCAAGTAATCGTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCC‑GATACAAACAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGcc                                <  1:289577/149‑1 (MQ=255)
                                                                                                                     tGGATACCAGCAAGTAATCTTTTTTTTTTTTTTTTTAAAGTATTAGCGGTCTGCGTTAAGACTTTTTTCATGGGTGGCG‑GATACAAAAAAGGCCTCCGGCTTTTACCGCTGCGGCCGGTTTCGGGCGCTTTTTTCCATTACCGCAGctt                               >  2:276577/1‑147 (MQ=19)
                                                                                                                            cAGCAAGTAATTTTTTTTGCTTTAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGccc                        <  1:496600/149‑1 (MQ=255)
                                                                                                                            cAGCAAGTAATCTTTTTTTTTTTTTTTTTAAAGTATTAACGTCCTGCGTTAAGACTTTTTTCATGGGGGCCG‑AATACAGGAAAGGCCGCAGGCTGTTACCCCTGAGGCCGTTTTCGTTCGCTTATTGCCATAACGGCCGCCTGACGccc                        >  2:674642/1‑149 (MQ=16)
                                                                                                                            cAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCTGCCGGAATCGGGCGCAAATTGCCATCACGGCAGCCTGACGccc                        >  1:849867/1‑149 (MQ=255)
                                                                                                                             aGCCAGTAGTCTTTTTCGCCTGAAAAATAAAGTATTAGCGGTCTGCGTCAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCCGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCg                       <  2:150841/149‑1 (MQ=25)
                                                                                                                             aGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGCCGCCCg                       >  2:65094/1‑149 (MQ=255)
                                                                                                                               cAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGtt                     >  1:1075693/1‑149 (MQ=255)
                                                                                                                                   tAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCa                 >  1:360667/1‑149 (MQ=255)
                                                                                                                                   tAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCATCCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCa                 <  1:1273980/149‑1 (MQ=255)
                                                                                                                                    aaTCTTTTTTGCTTGAAAAATAAAGTATTAGCGGTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCAc                <  1:239849/149‑1 (MQ=255)
                                                                                                                                      tCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTAATGGGTGGCG‑TATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCCCGGCCGCCTGACGCCCGTTTTCACCt              >  2:622692/1‑149 (MQ=25)
                                                                                                                                        ttttttTGTTGAAAAATAAAGTATTGGCTGTCTGCGTTAAGATCTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGTATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTa            <  2:60983/149‑1 (MQ=18)
                                                                                                                                        ttttttGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTa            <  1:1165482/149‑1 (MQ=255)
                                                                                                                                        tttgttGCTTGAAAAATAAAGTATTAGCGTTATGCGGTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCGGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTa            <  2:732154/149‑1 (MQ=255)
                                                                                                                                         tttttGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTAc           >  1:1186947/1‑149 (MQ=255)
                                                                                                                                              gCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCg      <  1:189352/149‑1 (MQ=255)
                                                                                                                                               cTTGAATCATAACGTAGTTGCGTTCTGTCGTATGCGTCTTCTCATGGCTGCCG‑GATACAAAAAAGGCCGCAGGCTGGTACCCCGTCCGCCGGTTTCGCGCGCCTATTGCCATCACGTCAGCCTGACGCCCGTTTTCCCCTTACTTCCgg     <  2:1335417/149‑1 (MQ=255)
                                                                                                                                               cTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCgg     <  2:1204352/149‑1 (MQ=255)
                                                                                                                                               cTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCgg     <  1:627299/149‑1 (MQ=255)
                                                                                                                                                 tGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGtt   >  2:344832/1‑149 (MQ=255)
                                                                                                                                                  gAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTa  <  2:595343/149‑1 (MQ=255)
                                                                                                                                                  gAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTa  <  2:23499/149‑1 (MQ=255)
                                                                                                                                                  gAAAAATAAAGTATTAGCGTGCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTa  <  2:595566/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCCGGATCCGGAACATGCGAAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCG‑GATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTA  >  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: ‑