Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,856,333 T→G 6.6% E199A (GAG→GCG) ‡ ydjH ← putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,856,3330TG 6.6% 257.3 / 10.6 76E199A (GAG→GCG) ‡ydjHputative kinase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base G (3/2);  total (38/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 = 1.93e-01

CTGCGACCGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTACC  >  NC_000913/1856184‑1856481
                                                                                                                                                     |                                                                                                                                                    
cTGCGACCGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCt                                                                                                                                                      >  3:540825/1‑149 (MQ=255)
   cGACCGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCg                                                                                                                                                    <  2:681450/149‑1 (MQ=255)
      ccgccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCa                                                                                                                                                 >  3:1071403/1‑149 (MQ=255)
         ccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                              <  4:405743/149‑1 (MQ=255)
               ccTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGa                                                                                                                                        <  2:584835/149‑1 (MQ=255)
                 ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAGAAAACAAGCCGCAATTTCATCCAGTATTTCTTTCCCGGTGAGTTATTTTTCCTCGGCAAAATTaaaacgg                                                                                                                                     >  1:905335/1‑143 (MQ=255)
                   cATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACa                                                                                                                                    <  4:1023888/149‑1 (MQ=255)
                   cATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACa                                                                                                                                    <  1:442377/149‑1 (MQ=255)
                       gTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATaa                                                                                                                               <  4:540825/149‑1 (MQ=255)
                         cATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATATTGCCTCGGCAAAATTAGGAAACAGATAAt                                                                                                                              >  1:745424/1‑149 (MQ=255)
                          aTGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATc                                                                                                                             <  2:280678/149‑1 (MQ=255)
                          aTGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGATGAGTAATTTTGCCTCGGCAATATTAGGAAACAGATAATc                                                                                                                             <  1:334745/149‑1 (MQ=255)
                            gTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           <  1:430728/149‑1 (MQ=255)
                               accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  4:1010644/1‑149 (MQ=255)
                               accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  2:1032942/1‑149 (MQ=255)
                               accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  2:577199/1‑149 (MQ=255)
                               accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  1:595023/1‑149 (MQ=255)
                                 cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           <  2:71308/144‑1 (MQ=255)
                                 cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           >  1:71308/1‑144 (MQ=255)
                                    gCTTGATAAAGCAGCCGTCTTTACCCGCTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                                                                   >  3:393782/1‑149 (MQ=255)
                                               cagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCa                                                                                                        >  2:812473/1‑149 (MQ=255)
                                                   cGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAata                                                                                                    <  3:274866/149‑1 (MQ=255)
                                                    gTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAatat                                                                                                   >  2:645934/1‑149 (MQ=255)
                                                     tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatc                                                                                                  <  3:570410/149‑1 (MQ=255)
                                                       tttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                                                                          >  3:252370/1‑106 (MQ=255)
                                                       tttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                                                                          <  4:252370/107‑2 (MQ=255)
                                                        ttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                                                                               <  2:819518/149‑1 (MQ=255)
                                                            ccGTTTTAATCACCACCGTTTTTACGCCGCACTCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGc                                                                                           <  4:816848/149‑1 (MQ=255)
                                                              gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGt                                                                                         >  4:181630/1‑149 (MQ=255)
                                                              gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGt                                                                                         >  4:248166/1‑149 (MQ=255)
                                                                 ttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttc                                                                                      <  3:653431/149‑1 (MQ=255)
                                                                           ccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTGGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                                                                            >  2:688128/1‑149 (MQ=255)
                                                                               ttttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGtt                                                                        <  4:385780/149‑1 (MQ=255)
                                                                                 ttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCgg                                                                                                                                                   >  3:29470/1‑72 (MQ=255)
                                                                                 ttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCgg                                                                                                                                                   <  4:29470/72‑1 (MQ=255)
                                                                                      ccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCat                                                                 <  1:812473/149‑1 (MQ=255)
                                                                                        gcacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCatat                                                               <  2:1149843/149‑1 (MQ=255)
                                                                                        gcacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCatat                                                               <  3:1010644/149‑1 (MQ=255)
                                                                                             caAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGca                                                          <  4:1071403/149‑1 (MQ=255)
                                                                                              aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGt                                                                                         <  4:403826/117‑1 (MQ=255)
                                                                                              aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGt                                                                                         >  3:403826/1‑117 (MQ=255)
                                                                                                gAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAg                                                       >  3:601032/1‑149 (MQ=255)
                                                                                                  aaGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt                                                     >  2:209743/1‑149 (MQ=255)
                                                                                                          aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCt                                             >  3:769773/1‑149 (MQ=255)
                                                                                                               tttCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCgg                                        <  3:248166/149‑1 (MQ=255)
                                                                                                                 tCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGc                                      <  1:577199/149‑1 (MQ=255)
                                                                                                                   aTCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCtt                                    >  1:797883/1‑149 (MQ=255)
                                                                                                                    tCCAGTGTCTCTTTCCCGGTGAGTCATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCACTGCTTCGCAAATATCCTCCCGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCttt                                   >  3:973287/1‑149 (MQ=255)
                                                                                                                     ccAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATCCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTa                                  <  1:645934/149‑1 (MQ=255)
                                                                                                                         tgtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCtt                              >  1:872168/1‑149 (MQ=255)
                                                                                                                           tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            >  1:1054065/1‑149 (MQ=255)
                                                                                                                           tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            <  2:595023/149‑1 (MQ=255)
                                                                                                                             tctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatc                                                 >  3:485455/1‑126 (MQ=255)
                                                                                                                             tctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatc                                                 <  4:485455/126‑1 (MQ=255)
                                                                                                                               tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                                                                    <  4:674109/105‑1 (MQ=255)
                                                                                                                               tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                                                                    >  3:674109/1‑105 (MQ=255)
                                                                                                                               tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATaa                                                   >  2:8361/1‑122 (MQ=255)
                                                                                                                               tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATaa                                                   <  1:8361/122‑1 (MQ=255)
                                                                                                                               tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTaa                        >  1:222510/1‑149 (MQ=255)
                                                                                                                                      gTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTc                 >  1:164251/1‑149 (MQ=255)
                                                                                                                                       tGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCtt                                    >  4:925168/1‑129 (MQ=255)
                                                                                                                                       tGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCtt                                    <  3:925168/129‑1 (MQ=255)
                                                                                                                                        gAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTg               >  4:408456/1‑149 (MQ=255)
                                                                                                                                         aGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                <  4:638768/51‑1 (MQ=255)
                                                                                                                                         aGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                >  3:638768/1‑51 (MQ=255)
                                                                                                                                         aGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              <  2:1054065/149‑1 (MQ=255)
                                                                                                                                         aGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              <  2:797883/149‑1 (MQ=255)
                                                                                                                                          gTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGtt             >  3:1066877/1‑149 (MQ=255)
                                                                                                                                          ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGttt                                                                       <  1:1048346/90‑1 (MQ=255)
                                                                                                                                          ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGttt                                                                       >  2:1048346/2‑91 (MQ=255)
                                                                                                                                           tAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATACAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa            >  3:876965/1‑149 (MQ=255)
                                                                                                                                           tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa            >  1:1013940/1‑149 (MQ=255)
                                                                                                                                             aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              <  3:700916/145‑1 (MQ=255)
                                                                                                                                             aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              >  4:700916/1‑145 (MQ=255)
                                                                                                                                                   ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTa    <  1:446325/149‑1 (MQ=255)
                                                                                                                                                    cTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt                <  4:362773/136‑1 (MQ=255)
                                                                                                                                                    cTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt                >  3:362773/1‑136 (MQ=255)
                                                                                                                                                     tCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAcc  <  2:222510/149‑1 (MQ=255)
                                                                                                                                                     tCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAcc  >  4:1134211/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCGACCGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTACC  >  NC_000913/1856184‑1856481

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