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.

TCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGGCGCAACCCGG  >  NC_000913/3701703‑3701998
                                                                                                                                                    |                                                                                                                                                    
tCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGccc                                                                                                                                                     >  4:1005929/1‑149 (MQ=255)
tCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGccc                                                                                                                                                     >  2:146741/1‑149 (MQ=255)
tCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGccc                                                                                                                                                     <  2:876882/149‑1 (MQ=255)
    gCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGt                                                                                                                                                 <  1:51428/149‑1 (MQ=255)
     ccTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGtt                                                                                                                                                <  3:566402/149‑1 (MQ=255)
      cTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCGCTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGttt                                                                                                                                               >  1:317187/1‑149 (MQ=255)
       tCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGtttt                                                                                                                                              >  4:878772/1‑149 (MQ=255)
        cATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTtg                                                                                                                                             <  3:789254/149‑1 (MQ=255)
         aTTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGCCATGCGCTCCGGtttac                                                                                                                                             >  3:1666734/1‑146 (MQ=255)
         aTTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTtgt                                                                                                                                            <  4:1593278/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTTTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGa                                                                                                                                          <  2:1324684/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGa                                                                                                                                          <  3:1005929/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGa                                                                                                                                          <  2:1368097/149‑1 (MQ=255)
           tCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGa                                                                                                                                          <  1:470015/149‑1 (MQ=255)
              gTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGtt                                                                                                                                       <  1:1254455/149‑1 (MQ=255)
                 tttCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAAc                                                                                                                                    <  2:1649715/149‑1 (MQ=255)
                   tCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGc                                                                                                                                  <  4:442125/149‑1 (MQ=255)
                      tGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGc                                                                                                                               >  2:493586/1‑149 (MQ=255)
                       gATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCg                                                                                                                              >  1:1499899/1‑149 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTACAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTACCCCGGAGCTCATGGCCCCGGTTTTGTGAGTTAGCGCTGCGGAtt                                                                                                                          <  3:249914/149‑1 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGa                                                                                                                                          <  1:886525/133‑1 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAAc                                                                                                                                    >  1:939128/1‑139 (MQ=255)
                           tataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGAtt                                                                                                                          >  3:1580593/1‑149 (MQ=255)
                            ataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGAttt                                                                                                                         <  1:1110793/149‑1 (MQ=255)
                             taAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGAtttt                                                                                                                        >  4:304418/1‑149 (MQ=255)
                               aaGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCa                                                                                                                      <  3:867646/149‑1 (MQ=255)
                               aaGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCa                                                                                                                      <  3:468843/149‑1 (MQ=255)
                                 gTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    >  1:289280/1‑149 (MQ=255)
                                  tCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTtgt                                                                                                                                            <  2:814828/124‑1 (MQ=255)
                                    gCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTg                                                                                                                 >  3:211260/1‑149 (MQ=255)
                                      acaACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                               >  2:1615648/1‑149 (MQ=255)
                                          cTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    >  2:734140/1‑140 (MQ=255)
                                           tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                          <  4:1276008/149‑1 (MQ=255)
                                           tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                          >  1:1674813/1‑149 (MQ=255)
                                           tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                          <  3:293579/149‑1 (MQ=255)
                                            ggAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATACTGATCAACAg                                                                                                         >  1:1749586/1‑149 (MQ=255)
                                             gAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGc                                                                                                        >  4:785619/1‑149 (MQ=255)
                                                aaaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaa                                                                                                     >  2:199217/1‑149 (MQ=255)
                                                 aaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaaa                                                                                                    <  3:737071/149‑1 (MQ=255)
                                                 aaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaaa                                                                                                    <  3:1537590/149‑1 (MQ=255)
                                                     tCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAg                                                                                                <  4:1132405/149‑1 (MQ=255)
                                                      cgcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                               <  2:1749586/149‑1 (MQ=255)
                                                      cgcaGATTTAAAGCCTGGCTGATTCCGGTGTACGGCATTTTGCTATGCTTGTCGCCATTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCTACGGTTTTGTGAGTTAACGCTGCGGATTTTCAGCCTGATCAACAGCAAACCaaac                                                                                               <  2:120213/149‑3 (MQ=255)
                                                       gcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                               >  4:610282/1‑148 (MQ=255)
                                                       gcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                               <  3:610282/148‑1 (MQ=255)
                                                       gcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCt                                                                                              >  1:276308/1‑149 (MQ=255)
                                                            tttAAAGCCTGGCTGATTCTGCTGTACAGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAg                                                                                         <  1:143322/149‑1 (MQ=255)
                                                              tAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                                       <  1:1585986/149‑1 (MQ=255)
                                                                    cTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGAccgc                                                                                 >  4:832252/1‑149 (MQ=255)
                                                                          gATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGt                                                                           <  2:1392501/149‑1 (MQ=255)
                                                                           aTTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTc                                                                          <  3:881271/149‑1 (MQ=255)
                                                                            ttCTGCTTTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                               >  2:936401/1‑111 (MQ=255)
                                                                            ttCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCt                                                                         >  1:548878/1‑149 (MQ=255)
                                                                              ctgctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCttt                                                                       <  3:969058/149‑1 (MQ=255)
                                                                               tgctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCtttt                                                                      <  1:68671/149‑1 (MQ=255)
                                                                               tgctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCtttt                                                                      <  4:1298077/149‑1 (MQ=255)
                                                                                gctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTa                                                                     <  1:1190752/149‑1 (MQ=255)
                                                                                    tACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTg                                                                 >  3:368013/1‑149 (MQ=255)
                                                                                     aCGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGc                                                                <  3:1233302/149‑1 (MQ=255)
                                                                                       ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAAc                                                                                                           <  1:813898/104‑1 (MQ=255)
                                                                                       ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCgg                                                              <  2:673458/149‑1 (MQ=255)
                                                                                         cATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCt                                                            <  3:672242/149‑1 (MQ=255)
                                                                                          aTTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTg                                                           >  2:1624521/1‑149 (MQ=255)
                                                                                           ttttGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGc                                                          >  4:348908/1‑149 (MQ=255)
                                                                                            tttGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGAttt                                                                                                                         <  4:561811/85‑1 (MQ=255)
                                                                                              tGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAAc                                                       <  1:1624521/149‑1 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    <  4:1061241/149‑1 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    >  2:689366/1‑149 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    >  3:1247443/1‑149 (MQ=255)
                                                                                                 tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                                    >  3:1150891/1‑149 (MQ=255)
                                                                                                      ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    <  3:1405665/80‑1 (MQ=255)
                                                                                                      ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                                    >  4:1405665/1‑80 (MQ=255)
                                                                                                      ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCa                                               >  1:1331839/1‑149 (MQ=255)
                                                                                                        gTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAgg                                             <  2:276308/149‑1 (MQ=255)
                                                                                                              aCTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTAc                                                                                            >  4:1537589/1‑96 (MQ=255)
                                                                                                                  ttGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAATAAGCTACCAGTTGACCGTCGTAGTCTTTTAGCTGCGGCTGCAACTGGATGCAGGGGCCGAAgcg                                   >  4:1029455/1‑149 (MQ=255)
                                                                                                                  ttGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAgcg                                   >  4:1153630/1‑149 (MQ=255)
                                                                                                                  ttGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAgcg                                   >  2:431625/1‑149 (MQ=255)
                                                                                                                    gAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGcc                                 <  2:1116445/149‑1 (MQ=255)
                                                                                                                       gCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATGCTGATCAACAGCAAAACAAGCTAc                                                                                            <  2:59749/87‑1 (MQ=255)
                                                                                                                         tGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCcgacga                            >  2:95624/1‑149 (MQ=255)
                                                                                                                           aTTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                          >  1:820701/1‑149 (MQ=255)
                                                                                                                           aTTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                          <  2:289280/149‑1 (MQ=255)
                                                                                                                           aTTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                          <  4:1580593/149‑1 (MQ=255)
                                                                                                                             ttgatAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGc                        <  1:96906/145‑1 (MQ=255)
                                                                                                                             tcacaAAACCGGGGCCATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGc                        <  2:375984‑M1/144‑1 (MQ=255)
                                                                                                                               acaAAACCGGGGCCATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                                       >  3:25757‑M1/4‑84 (MQ=255)
                                                                                                                               acaAAACCGGGGCCATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                                       <  4:25757‑M1/81‑1 (MQ=255)
                                                                                                                               acaAAACCGGGGCCATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCgg                      <  1:681240‑M1/146‑1 (MQ=255)
                                                                                                                                gTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCggg                     >  3:823436/1‑149 (MQ=255)
                                                                                                                                 tAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGc                    >  1:1040749/1‑149 (MQ=255)
                                                                                                                                 aAAACCGGGGCCATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGc                    >  1:375984‑M1/2‑149 (MQ=255)
                                                                                                                                    ACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGtt                 >  3:324404‑M1/1‑149 (MQ=255)
                                                                                                                                       GGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGttgga              >  2:1525801‑M1/1‑147 (MQ=255)
                                                                                                                                             cATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAggg                                            <  3:1297637/113‑1 (MQ=255)
                                                                                                                                             cATGGCCCCGGGTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCATTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCTACGACAGCGGGCGTTGAAGGCGCa        >  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)
                                                                                                                                                    |                                                                                                                                                    
TCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATGGCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGGCGCAACCCGG  >  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: