Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,701,850 C→T 6.0% intergenic (+16/+14) yhjV → / ← dppF putative transporter/dipeptide/heme ABC transporter ATPas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,701,8500CT 6.0% 698.6 / 21.4 215intergenic (+16/+14)yhjV/dppFputative transporter/dipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base C (85/116);  new base T (4/9);  total (89/126)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

Only 100 of 228 total aligned reads displayed.

TCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGGCGCAACCCGG  >  NC_000913/3701703‑3701998
                                                                                                                                                    |                                                                                                                                                    
tCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑Gccc                                                                                                                                                     >  4:1005929/1‑149 (MQ=255)
tCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑Gccc                                                                                                                                                     >  2:146741/1‑149 (MQ=255)
tCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑Gccc                                                                                                                                                     <  2:876882/149‑1 (MQ=255)
    gCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGt                                                                                                                                                 <  1:51428/149‑1 (MQ=255)
     ccTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGtt                                                                                                                                                <  3:566402/149‑1 (MQ=255)
      cTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCGCTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGttt                                                                                                                                               >  1:317187/1‑149 (MQ=255)
       tCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGtttt                                                                                                                                              >  4:878772/1‑149 (MQ=255)
        cATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTtg                                                                                                                                             <  3:789254/149‑1 (MQ=255)
         aTTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGC‑CATGCGCTCCGGtttac                                                                                                                                             >  3:1666734/1‑146 (MQ=255)
         aTTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTtgt                                                                                                                                            <  4:1593278/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTTTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                          <  2:1324684/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                          <  3:1005929/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                          <  2:1368097/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                          <  1:470015/149‑1 (MQ=255)
              gTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGtt                                                                                                                                       <  1:1254455/149‑1 (MQ=255)
                 tttCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAAc                                                                                                                                    <  2:1649715/149‑1 (MQ=255)
                   tCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGc                                                                                                                                  <  4:442125/149‑1 (MQ=255)
                      tGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGc                                                                                                                               >  2:493586/1‑149 (MQ=255)
                       gATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCg                                                                                                                              >  1:1499899/1‑149 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTACAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTACCCCGGAGCTCATG‑GCCCCGGTTTTGTGAGTTAGCGCTGCGGAtt                                                                                                                          <  3:249914/149‑1 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                          <  1:886525/133‑1 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAAc                                                                                                                                    >  1:939128/1‑139 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAtt                                                                                                                          >  3:1580593/1‑149 (MQ=255)
                            ataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAttt                                                                                                                         <  1:1110793/149‑1 (MQ=255)
                             taAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAtttt                                                                                                                        >  4:304418/1‑149 (MQ=255)
                               aaGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCa                                                                                                                      <  3:867646/149‑1 (MQ=255)
                               aaGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCa                                                                                                                      <  3:468843/149‑1 (MQ=255)
                                 gTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    >  1:289280/1‑149 (MQ=255)
                                  tCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTtgt                                                                                                                                            <  2:814828/124‑1 (MQ=255)
                                    gCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTg                                                                                                                 >  3:211260/1‑149 (MQ=255)
                                      acaACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                               >  2:1615648/1‑149 (MQ=255)
                                          cTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    >  2:734140/1‑140 (MQ=255)
                                           tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                          <  4:1276008/149‑1 (MQ=255)
                                           tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                          >  1:1674813/1‑149 (MQ=255)
                                           tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                          <  3:293579/149‑1 (MQ=255)
                                            ggAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATACTGATCAACAg                                                                                                         >  1:1749586/1‑149 (MQ=255)
                                             gAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGc                                                                                                        >  4:785619/1‑149 (MQ=255)
                                                aaaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaa                                                                                                     >  2:199217/1‑149 (MQ=255)
                                                 aaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaaa                                                                                                    <  3:737071/149‑1 (MQ=255)
                                                 aaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaaa                                                                                                    <  3:1537590/149‑1 (MQ=255)
                                                     tCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAg                                                                                                <  4:1132405/149‑1 (MQ=255)
                                                      cgcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                               <  2:1749586/149‑1 (MQ=255)
                                                      cgcaGATTTAAAGCCTGGCTGATTCCGGTGTACGGCATTTTGCTATGCTTGTCGCCATTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCTACGGTTTTGTGAGTTAACGCTGCGGATTTTCAGCCTGATCAACAGCAAACCaaac                                                                                               <  2:120213/149‑3 (MQ=255)
                                                       gcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                               >  4:610282/1‑148 (MQ=255)
                                                       gcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                               <  3:610282/148‑1 (MQ=255)
                                                       gcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCt                                                                                              >  1:276308/1‑149 (MQ=255)
                                                            tttAAAGCCTGGCTGATTCTGCTGTACAGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAg                                                                                         <  1:143322/149‑1 (MQ=255)
                                                              tAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                                       <  1:1585986/149‑1 (MQ=255)
                                                                    cTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGAccgc                                                                                 >  4:832252/1‑149 (MQ=255)
                                                                          gATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGt                                                                           <  2:1392501/149‑1 (MQ=255)
                                                                           aTTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTc                                                                          <  3:881271/149‑1 (MQ=255)
                                                                            ttCTGCTTTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                               >  2:936401/1‑111 (MQ=255)
                                                                            ttCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCt                                                                         >  1:548878/1‑149 (MQ=255)
                                                                              ctgctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCttt                                                                       <  3:969058/149‑1 (MQ=255)
                                                                               tgctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCtttt                                                                      <  1:68671/149‑1 (MQ=255)
                                                                               tgctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCtttt                                                                      <  4:1298077/149‑1 (MQ=255)
                                                                                gctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTa                                                                     <  1:1190752/149‑1 (MQ=255)
                                                                                    tACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTg                                                                 >  3:368013/1‑149 (MQ=255)
                                                                                     aCGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGc                                                                <  3:1233302/149‑1 (MQ=255)
                                                                                       ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAAc                                                                                                           <  1:813898/104‑1 (MQ=255)
                                                                                       ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCgg                                                              <  2:673458/149‑1 (MQ=255)
                                                                                         cATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCt                                                            <  3:672242/149‑1 (MQ=255)
                                                                                          aTTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTg                                                           >  2:1624521/1‑149 (MQ=255)
                                                                                           ttttGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGc                                                          >  4:348908/1‑149 (MQ=255)
                                                                                            tttGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAttt                                                                                                                         <  4:561811/85‑1 (MQ=255)
                                                                                              tGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAAc                                                       <  1:1624521/149‑1 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    <  4:1061241/149‑1 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    >  2:689366/1‑149 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    >  3:1247443/1‑149 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    >  3:1150891/1‑149 (MQ=255)
                                                                                                      ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    <  3:1405665/80‑1 (MQ=255)
                                                                                                      ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    >  4:1405665/1‑80 (MQ=255)
                                                                                                      ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCa                                               >  1:1331839/1‑149 (MQ=255)
                                                                                                        gTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAgg                                             <  2:276308/149‑1 (MQ=255)
                                                                                                              aCTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTAc                                                                                            >  4:1537589/1‑96 (MQ=255)
                                                                                                                  ttGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAATAAGCTACCAGTTGACCGTCGTAGTCTTTTAGCTGCGGCTGCAACTGGATGCAGGGGCCGAAgcg                                   >  4:1029455/1‑149 (MQ=255)
                                                                                                                  ttGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAgcg                                   >  4:1153630/1‑149 (MQ=255)
                                                                                                                  ttGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAgcg                                   >  2:431625/1‑149 (MQ=255)
                                                                                                                    gAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGcc                                 <  2:1116445/149‑1 (MQ=255)
                                                                                                                       gCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATGCTGATCAACAGCAAAACAAGCTAc                                                                                            <  2:59749/87‑1 (MQ=255)
                                                                                                                         tGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCcgacga                            >  2:95624/1‑149 (MQ=255)
                                                                                                                           aTTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                          >  1:820701/1‑149 (MQ=255)
                                                                                                                           aTTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                          <  2:289280/149‑1 (MQ=255)
                                                                                                                           aTTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                          <  4:1580593/149‑1 (MQ=255)
                                                                                                                             ttgatAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGc                        <  1:96906/145‑1 (MQ=255)
                                                                                                                             tcacaAAACCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGc                        <  2:375984‑M1/144‑1 (MQ=255)
                                                                                                                               acaAAACCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                                       >  3:25757‑M1/4‑84 (MQ=255)
                                                                                                                               acaAAACCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                                       <  4:25757‑M1/81‑1 (MQ=255)
                                                                                                                               acaAAACCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCgg                      <  1:681240‑M1/146‑1 (MQ=255)
                                                                                                                                gTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCggg                     >  3:823436/1‑149 (MQ=255)
                                                                                                                                 tAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGc                    >  1:1040749/1‑149 (MQ=255)
                                                                                                                                 aAAACCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGc                    >  1:375984‑M1/2‑149 (MQ=255)
                                                                                                                                    ACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGtt                 >  3:324404‑M1/1‑149 (MQ=255)
                                                                                                                                       GGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGttgga              >  2:1525801‑M1/1‑147 (MQ=255)
                                                                                                                                             cATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAggg                                            <  3:1297637/113‑1 (MQ=255)
                                                                                                                                             cATG‑GCCCCGGGTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCATTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCTACGACAGCGGGCGTTGAAGGCGCa        >  3:1514314/1‑149 (MQ=255)
                                                                                                                                                   ccccGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGGCGCAACCcg   >  1:858985/1‑149 (MQ=255)
                                                                                                                                                    cccGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATc                                                                                                              <  3:1173263/41‑1 (MQ=255)
                                                                                                                                                    cccGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGaa                                      <  4:855120/113‑1 (MQ=255)
                                                                                                                                                    cccGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGGCGCAACCcgg  >  4:1196224/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGGCGCAACCCGG  >  NC_000913/3701703‑3701998

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