Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,701,840 A→G 7.1% intergenic (+6/+24) 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,8400AG 7.1% 395.1 / 11.4 126intergenic (+6/+24)yhjV/dppFputative transporter/dipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base A (63/54);  new base G (6/3);  total (69/57)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 132 total aligned reads displayed.

ATGGTATTGTGTCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGG  >  NC_000913/3701692‑3701988
                                                                                                                                                    |                                                                                                                                                     
atGGTATTGTGTCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGa                                                                                                                                                       <  1:110718/149‑1 (MQ=255)
atGGTATTGTGTCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGa                                                                                                                                                       <  2:475790/149‑1 (MQ=255)
atGGTATTGTGTCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGa                                                                                                                                                       <  3:185460/149‑1 (MQ=255)
      ttGTGTCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCAt                                                                                                                                                 >  2:780646/1‑149 (MQ=255)
        gtgtCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACGGTTGAAGCTGATTGAGTAAACCGGAGCGCATg‑g                                                                                                                                              <  1:657580/149‑1 (MQ=255)
         tgtCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑Gc                                                                                                                                             >  4:428396/1‑149 (MQ=255)
             gTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCg                                                                                                                                         >  2:580795/1‑149 (MQ=255)
             gTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTACTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGTAGCGCATG‑GCCCCg                                                                                                                                         >  2:497954/1‑149 (MQ=255)
                ccTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGtt                                                                                                                                      <  2:531396/149‑1 (MQ=255)
                 cTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGttt                                                                                                                                     >  3:213807/1‑149 (MQ=255)
                 cTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGttt                                                                                                                                     >  1:523208/1‑149 (MQ=255)
                  tCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGC‑CATGCGCTCCGGttta                                                                                                                                    <  1:475175/149‑2 (MQ=255)
                  tCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGtttt                                                                                                                                    <  2:383797/149‑1 (MQ=255)
                  tCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGtttt                                                                                                                                    >  3:533657/1‑149 (MQ=255)
                   cATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTtg                                                                                                                                   <  4:242564/149‑1 (MQ=255)
                    aTTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTtgt                                                                                                                                  <  4:239053/149‑1 (MQ=255)
                         gTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGtt                                                                                                                             >  3:477608/1‑149 (MQ=255)
                                 tGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGc                                                                                                                     >  4:841177/1‑149 (MQ=255)
                                 tGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGc                                                                                                                     >  3:283305/1‑149 (MQ=255)
                                 tGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGc                                                                                                                     >  4:392650/1‑149 (MQ=255)
                                  gATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCg                                                                                                                    <  2:929467/149‑1 (MQ=255)
                                  gATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCg                                                                                                                    >  3:830952/1‑149 (MQ=255)
                                       ataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGa                                                                                                                                                       <  3:619049/110‑1 (MQ=255)
                                       ataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAttt                                                                                                               <  4:283305/149‑1 (MQ=255)
                                       ataAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAttt                                                                                                               >  2:849993/1‑149 (MQ=255)
                                        taAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAtttt                                                                                                              <  3:369765/149‑1 (MQ=255)
                                        taAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGAtttt                                                                                                              <  3:706109/149‑1 (MQ=255)
                                         aaaGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTc                                                                                                             >  4:17091/1‑149 (MQ=255)
                                           aGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCAt                                                                                                           <  4:368906/149‑1 (MQ=255)
                                              cGCACAACTGGAAAAACTTTGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTTATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCt                                                                                                        <  2:270569/149‑1 (MQ=255)
                                                cacaACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGTATTTTCATCCTGa                                                                                                      >  2:626229/1‑149 (MQ=255)
                                                 acaACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGC‑CATGCGCTCCGGttta                                                                                                                                    >  4:6689/1‑117 (MQ=255)
                                                 acaACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGC‑CATGCGCTCCGGttta                                                                                                                                    <  3:6689/118‑2 (MQ=255)
                                                 acaACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                     <  3:692839/149‑1 (MQ=255)
                                                   aaCTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCa                                                                                                   >  2:338188/1‑149 (MQ=255)
                                                   aaCTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCa                                                                                                   >  1:585790/1‑149 (MQ=255)
                                                    aCTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGC‑CATGCGCTCCGGttta                                                                                                                                    <  1:630434/115‑2 (MQ=255)
                                                    aCTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGGGC‑CATGCGCTCCGGttta                                                                                                                                    >  2:630434/1‑114 (MQ=255)
                                                     cTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAAc                                                                                                 <  2:582651/149‑1 (MQ=255)
                                                      tGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACa                                                                                                <  3:17091/149‑1 (MQ=255)
                                                       ggAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAg                                                                                               >  2:156872/1‑149 (MQ=255)
                                                        gAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                     >  1:892610/1‑142 (MQ=255)
                                                        gAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                     <  2:892610/142‑1 (MQ=255)
                                                        gAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGc                                                                                              >  3:219244/1‑149 (MQ=255)
                                                         aaaaaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCa                                                                                             >  3:900788/1‑149 (MQ=255)
                                                            aaCTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAACCCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCaaaa                                                                                          <  4:533657/149‑1 (MQ=255)
                                                               ttCGTGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACaa                                                                                       <  1:351079/149‑1 (MQ=255)
                                                                tCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAg                                                                                      <  4:744711/149‑1 (MQ=255)
                                                                tCGCGGATTTAAAGCCTGGCTGATTCTACTGTACGGCATTTTGCTATGCTTGTCGCCACGGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCCGCGGATTTTCATCCTGATCAACAGCAAAACAAg                                                                                      <  1:497954/149‑1 (MQ=255)
                                                                 cgcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                     <  4:466099/149‑1 (MQ=255)
                                                                 cgcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                     <  1:125639/149‑1 (MQ=255)
                                                                 cgcgGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACCCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGc                                                                                     <  4:876195/149‑1 (MQ=255)
                                                                    ggATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTAc                                                                                  <  1:780646/149‑1 (MQ=255)
                                                                      aTTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCa                                                                                >  1:17025/1‑149 (MQ=255)
                                                                      aTTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCa                                                                                >  3:174617/1‑149 (MQ=255)
                                                                       tttAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAg                                                                               >  2:200601/1‑149 (MQ=255)
                                                                           aaGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGa                                                                           <  3:822664/149‑1 (MQ=255)
                                                                               cTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGAccgc                                                                       <  3:151242/149‑1 (MQ=255)
                                                                                       ttCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTaa                                                                                                                           <  2:101887/89‑1 (MQ=255)
                                                                                       ttCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTaa                                                                                                                           >  1:101887/1‑89 (MQ=255)
                                                                                        tCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCtt                                                              >  3:433768/1‑149 (MQ=255)
                                                                                        tCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCtt                                                              <  2:835018/149‑1 (MQ=255)
                                                                                           gctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTa                                                           <  1:750328/149‑1 (MQ=255)
                                                                                           gctgTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTa                                                           <  2:904852/149‑1 (MQ=255)
                                                                                              gTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAAAGCTGCGGATTTTCATCCTTATCAACAGCAAAACAAGCTACCAGTTGACCGCTGTAGTCTTTTAGCt                                                        <  4:213807/149‑1 (MQ=255)
                                                                                               tACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTg                                                       >  1:14023/1‑149 (MQ=255)
                                                                                                  ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGTAGCGCATG‑GCCCCGGCTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAACTACCAGTTGACCGCCGTAGTCTTTTAGCTGCgg                                                    >  4:743158/1‑149 (MQ=255)
                                                                                                  ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                          <  4:733575/95‑1 (MQ=255)
                                                                                                  ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATc                                                                                                          >  3:733575/1‑95 (MQ=255)
                                                                                                  ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCgg                                                    >  4:231713/1‑149 (MQ=255)
                                                                                                  ggCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCgg                                                    >  2:406877/1‑149 (MQ=255)
                                                                                                    cATGTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCATTTAGCTGCGGCt                                                  >  4:787813/1‑149 (MQ=255)
                                                                                                       tttGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCa                                               >  3:311271/1‑149 (MQ=255)
                                                                                                            tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACGGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGa                                                                                                      >  1:767165/1‑89 (MQ=255)
                                                                                                            tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGa                                                                                                      <  2:767165/89‑1 (MQ=255)
                                                                                                            tATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTgg                                          <  2:621954/149‑1 (MQ=255)
                                                                                                              tGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGt                                        <  3:787813/149‑1 (MQ=255)
                                                                                                                 ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                <  1:868386/58‑1 (MQ=255)
                                                                                                                 ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGa                                                                                                                                >  2:868386/1‑58 (MQ=255)
                                                                                                                 ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                     >  4:553264/1‑85 (MQ=255)
                                                                                                                 ttGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGAt                                                                                                     <  3:553264/85‑1 (MQ=255)
                                                                                                                        cACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                             >  1:343540/1‑102 (MQ=255)
                                                                                                                        cACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGtt                                                                             <  2:343540/102‑1 (MQ=255)
                                                                                                                           tGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAg                           >  1:683268/1‑149 (MQ=255)
                                                                                                                            gTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAgc                          <  4:183725/149‑1 (MQ=255)
                                                                                                                             ttGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAgcg                         >  1:400684/1‑149 (MQ=255)
                                                                                                                                     gATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCgacgac                 >  1:462980/1‑149 (MQ=255)
                                                                                                                                      aTTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACa                >  2:168473/1‑149 (MQ=255)
                                                                                                                                       ttGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAg               >  4:423108/1‑149 (MQ=255)
                                                                                                                                        tcacaaaaCCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGAATTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGAGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGc              >  2:351309/6‑149 (MQ=255)
                                                                                                                                          acaaaaCCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCgg            >  4:883412/4‑149 (MQ=255)
                                                                                                                                          acaaaaCCGGGGC‑CATGCGCTCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCgg            >  2:515787/4‑149 (MQ=255)
                                                                                                                                           gTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCt                                                               >  3:418252/1‑97 (MQ=255)
                                                                                                                                            tAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGc          >  4:676711/1‑149 (MQ=255)
                                                                                                                                               aCCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGtt       >  3:667792/1‑149 (MQ=255)
                                                                                                                                               aCCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGtt       >  2:742792/1‑149 (MQ=255)
                                                                                                                                                 cGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGa     >  4:180570/1‑149 (MQ=255)
                                                                                                                                                   gAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAg   <  4:504962/149‑1 (MQ=255)
                                                                                                                                                   gAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAg   <  2:400684/149‑1 (MQ=255)
                                                                                                                                                    aGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAgg  <  3:15069/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGGTATTGTGTCGTGCCTCATTCCGTTTTTCCTGATCTATAAAGTCGCACAACTGGAAAAACTTCGCGGATTTAAAGCCTGGCTGATTCTGCTGTACGGCATTTTGCTATGCTTGTCGCCACTGTTGAAGCTGATTGAGTAAACCGGAGCGCATG‑GCCCCGGTTTTGTGAGTTAACGCTGCGGATTTTCATCCTGATCAACAGCAAAACAAGCTACCAGTTGACCGCCGTAGTCTTTTAGCTGCGGCTGCAACTGGGTGCAGGGGCCGAAGCGCCGACGACAGCGGGCGTTGAAGG  >  NC_000913/3701692‑3701988

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