Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,585 C→T 8.1% D108D (GAC→GAT) ‡ mcbR → colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,5850CT 8.1% 342.5 / 16.4 111D108D (GAC→GAT) ‡mcbRcolanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base C (53/49);  new base T (5/4);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 112 total aligned reads displayed.

TTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAG  >  NC_000913/1520437‑1520729
                                                                                                                                                     |                                                                                                                                                
tttCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACTGTTCCGGAAGTGGGGAAAAGTCAATTAGACGACAGCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGGTCTGGCTGTTGAAAACGTCACCCCGCAAGAc                                                                                                                                                  >  3:757193/1‑149 (MQ=255)
    ggTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCtt                                                                                                                                               >  2:364710/1‑148 (MQ=255)
     gTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGc                                                                                                                                             <  3:682771/149‑1 (MQ=255)
          tGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCggggt                                                                                                                                        <  3:272970/149‑4 (MQ=255)
           gcgcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAAc                                                                                                                                        >  2:421255/1‑148 (MQ=255)
            cgcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCggggtga                                                                                                                                      <  3:859961/149‑6 (MQ=255)
            cgcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACt                                                                                                                                       <  2:92615/148‑1 (MQ=255)
             gcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTcc                                                                                                                                     >  3:531480/1‑149 (MQ=255)
               tttCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCa                                                                                                                                    <  2:765271/148‑1 (MQ=255)
               tttCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAg                                                                                                                                   <  4:242257/149‑1 (MQ=255)
               tttCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAg                                                                                                                                   >  4:562153/1‑149 (MQ=255)
                        cACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTACTc                                                                                                                          >  3:132526/1‑149 (MQ=255)
                          ccTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCga                                                                                                                        <  4:339465/149‑1 (MQ=255)
                          ccTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCg                                                                                                                         <  1:421255/148‑1 (MQ=255)
                              cacaAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGaa                                                                                                                     >  1:733676/1‑148 (MQ=255)
                                caAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAt                                                                                                                   >  2:109264/1‑148 (MQ=255)
                                 aaGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAg                                                                                                                                   <  4:840866/131‑1 (MQ=255)
                                 aaGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAg                                                                                                                                   >  3:840866/1‑131 (MQ=255)
                                 aaGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtt                                                                                                                  <  1:109264/148‑1 (MQ=255)
                                   gCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                               <  3:16435/149‑1 (MQ=255)
                                   gCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                               <  4:549185/149‑1 (MQ=255)
                                     aTTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAcagc                                                                                                             <  4:531480/149‑1 (MQ=255)
                                        tACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAgg                                                                                                          <  3:163126/149‑1 (MQ=255)
                                         aCAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAgg                                                                                                          <  2:733676/148‑1 (MQ=255)
                                            gTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtt                                                                                                                  >  1:873167/1‑137 (MQ=255)
                                            gTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtt                                                                                                                  <  2:873167/137‑1 (MQ=255)
                                                  gAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCGGGCGCAAGaaaa                                                                                                >  3:351773/1‑149 (MQ=255)
                                                      tGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGaa                                                                                                                                         <  3:731525/104‑1 (MQ=255)
                                                       ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaa                                                                                                 <  1:807821/143‑1 (MQ=255)
                                                       ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaa                                                                                                 >  2:807821/1‑143 (MQ=255)
                                                         ggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                          >  1:539612/1‑148 (MQ=255)
                                                                tCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAc                                                                                  >  4:743222/1‑149 (MQ=255)
                                                                   aTTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTc                                                                                                                          >  3:571103/1‑106 (MQ=255)
                                                                   aTTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTc                                                                                                                          <  4:571103/106‑1 (MQ=255)
                                                                     taagcGAAATCCATCCGATCCCCCAAGGACTGGAATTAATGGGAAGTGCTTTGGGTGTTGAAAACCCCCCCCCGGCAGACCCTGGCGAAACTCAGGAATTGGCTGAAAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAAt                                                                              <  1:530129/144‑1 (MQ=255)
                                                                     tAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAAt                                                                              <  1:414784/148‑1 (MQ=255)
                                                                      aGACGAAATCACTCGGATCCGCTACGAACTCGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACATTGCGGAACTCCAGGAATTGCTCGAGATATTACAGCAGGCGCAAGAAAAGGTCGACATGGAACAAATCa                                                                            <  4:757193/149‑1 (MQ=255)
                                                                       gACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCa                                                                            <  2:375675/148‑1 (MQ=255)
                                                                          gAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        <  4:622524/149‑1 (MQ=255)
                                                                           aaaTCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        >  2:726657/1‑148 (MQ=255)
                                                                           aaaTCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        >  2:349878/1‑148 (MQ=255)
                                                                            aatcaatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATg                                                                      >  3:400676/1‑149 (MQ=255)
                                                                               caatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTaa                                                                   >  4:137206/1‑149 (MQ=255)
                                                                                 atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTaaa                                                                  <  2:792329/148‑1 (MQ=255)
                                                                                 atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTaaa                                                                  >  1:449872/1‑148 (MQ=255)
                                                                                 atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAAc                                                                 >  3:303264/1‑149 (MQ=255)
                                                                                  tcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCggggt                                                                                                                                        <  3:793326/77‑4 (MQ=255)
                                                                                  tcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCggggt                                                                                                                                        >  4:793326/1‑74 (MQ=255)
                                                                                   cGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGCCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACa                                                                >  2:427727/1‑148 (MQ=255)
                                                                                     gATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGc                                                                                                                            >  4:150235/1‑86 (MQ=255)
                                                                                        ccGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTAt                                                          >  4:341464/1‑149 (MQ=255)
                                                                                           cTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTc                                                       <  3:562153/149‑1 (MQ=255)
                                                                                             aCGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCg                                                      >  2:799045/1‑148 (MQ=255)
                                                                                             aCGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCACGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCg                                                      <  1:587950/148‑1 (MQ=255)
                                                                                               gAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCtt                                                   >  4:245652/1‑149 (MQ=255)
                                                                                                 aCTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAg                                                 >  4:835030/1‑149 (MQ=255)
                                                                                                   tGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCTCCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGgcgc                                                                                                       <  4:623597/93‑1 (MQ=255)
                                                                                                   tGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGgcgc                                                                                                       >  3:623597/1‑93 (MQ=255)
                                                                                                     gAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAt                                             >  3:24666/1‑149 (MQ=255)
                                                                                                     gAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAt                                             >  4:507623/1‑149 (MQ=255)
                                                                                                      aaTTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTaa                                                                   <  2:15064/126‑1 (MQ=255)
                                                                                                      aaTTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTaa                                                                   >  1:15064/1‑126 (MQ=255)
                                                                                                      aaTTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAt                                             <  1:349878/148‑1 (MQ=255)
                                                                                                       aTTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAtt                                            >  1:631313/1‑148 (MQ=255)
                                                                                                        ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGTCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGa                                              >  3:855314/1‑145 (MQ=255)
                                                                                                        ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGAc                                                                                         <  1:875332/102‑1 (MQ=255)
                                                                                                        ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGAc                                                                                         >  2:875332/1‑102 (MQ=255)
                                                                                                        ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGa                                              <  4:855314/145‑1 (MQ=255)
                                                                                                        ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAttt                                           >  2:278261/1‑148 (MQ=255)
                                                                                                        ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTa                                          >  3:439147/1‑149 (MQ=255)
                                                                                                          aaTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAttt                                           >  3:258781/1‑146 (MQ=255)
                                                                                                          aaTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAttt                                           <  4:258781/146‑1 (MQ=255)
                                                                                                            tGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatca                                       >  2:310067/1‑148 (MQ=255)
                                                                                                            tGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatca                                       >  1:240982/1‑148 (MQ=255)
                                                                                                               cAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGc                                   <  4:439147/149‑1 (MQ=255)
                                                                                                                         ggCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGc                         <  4:100363/149‑1 (MQ=255)
                                                                                                                         ggCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGATCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATg                          >  2:297413/1‑148 (MQ=255)
                                                                                                                          gCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGcc                        >  3:881332/1‑149 (MQ=255)
                                                                                                                             gTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCa                      <  1:115043/148‑1 (MQ=255)
                                                                                                                             gTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCa                      >  2:453920/1‑148 (MQ=255)
                                                                                                                              ttGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATAAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCAt                     <  2:539612/148‑1 (MQ=255)
                                                                                                                               tGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        >  1:34333/1‑96 (MQ=255)
                                                                                                                               tGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        <  2:34333/96‑1 (MQ=255)
                                                                                                                                gAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCt                  >  4:870714/1‑149 (MQ=255)
                                                                                                                                 aaaaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCt                  >  2:169854/1‑148 (MQ=255)
                                                                                                                                      cTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtgtg             <  1:726657/148‑1 (MQ=255)
                                                                                                                                      cTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga            <  3:811156/149‑1 (MQ=255)
                                                                                                                                      cTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga            >  4:235234/1‑149 (MQ=255)
                                                                                                                                       tCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgag           >  3:607865/1‑149 (MQ=255)
                                                                                                                                        cACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgaga          <  3:137206/149‑1 (MQ=255)
                                                                                                                                        cACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgag           <  1:310067/148‑1 (MQ=255)
                                                                                                                                          ttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACa                                                                                        >  4:59896/3‑69 (MQ=255)
                                                                                                                                          ttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACa                                                                                        <  3:59896/67‑1 (MQ=255)
                                                                                                                                          ttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgat         >  1:667715/3‑148 (MQ=255)
                                                                                                                                           tccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatg        >  1:299731/2‑148 (MQ=255)
                                                                                                                                            ccGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatgat      <  4:79939/149‑1 (MQ=255)
                                                                                                                                             cGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatgat      <  1:397315/148‑1 (MQ=255)
                                                                                                                                                 aGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga            >  2:742457/1‑138 (MQ=255)
                                                                                                                                                 aGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga            <  1:742457/138‑1 (MQ=255)
                                                                                                                                                 aGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAg  <  1:867442/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAG  >  NC_000913/1520437‑1520729

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: