Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 144,977 (A)7→8 coding (401/441 nt) yadI → putative PTS enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,9701.A99.3% 537.6 / ‑3.7 148coding (394/441 nt)yadIputative PTS enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (64/83);  total (65/83)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 185 total aligned reads displayed.

TGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCGTT  >  NC_000913/144878‑145057
                                                                                              |                                                                                       
tgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaac                                                                                          <  1:3528168/93‑1 (MQ=255)
 ggCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaacaaca                                                                                        >  2:2474089/1‑93 (MQ=255)
 ggCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaacaaca                                                                                        <  2:381114/93‑1 (MQ=255)
  gCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaacaacaa                                                                                       <  2:3434178/93‑1 (MQ=255)
   cTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaa                                                                                      <  1:572515/93‑1 (MQ=255)
    ttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaa                                                                                     >  2:4030085/1‑93 (MQ=255)
    ttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaa                                                                                     <  1:1446844/93‑1 (MQ=255)
    ttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaa                                                                                     >  2:4159797/1‑93 (MQ=255)
      gCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaa                                                                                   <  2:24725/93‑1 (MQ=255)
       ccGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaa                                                                                  <  1:4211750/93‑1 (MQ=255)
        cGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaa                                                                                 <  1:2817314/93‑2 (MQ=255)
         gTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaat                                                                                >  1:2261021/1‑91 (MQ=255)
         gTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaat                                                                                >  2:52030/1‑91 (MQ=255)
          tGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaatc                                                                               >  2:1034955/1‑90 (MQ=255)
           gAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaatcc                                                                              <  2:4613424/93‑5 (MQ=255)
             aaCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                                                             >  2:794891/1‑93 (MQ=255)
             aaCCATGACCAGGTCAGAGTTGCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                                                             <  2:3389581/93‑1 (MQ=255)
              aCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGcc                                                                            >  2:2191208/1‑93 (MQ=255)
               ccATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCt                                                                           >  2:2424344/1‑93 (MQ=255)
                cATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCtt                                                                          >  1:1678219/1‑93 (MQ=255)
                cATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCtt                                                                          >  2:2938394/1‑93 (MQ=255)
                  tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                                                        <  2:1313467/93‑1 (MQ=255)
                  tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                                                        >  2:2287865/1‑93 (MQ=255)
                  tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                                                        >  2:2648645/1‑93 (MQ=255)
                   gACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCg                                                                       <  2:3464464/93‑1 (MQ=255)
                   gACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCg                                                                       >  1:3331462/1‑93 (MQ=255)
                      cAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTcc                                                                    <  2:73550/93‑1 (MQ=255)
                      cAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTcc                                                                    >  1:1695498/1‑93 (MQ=255)
                         ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                                                 >  2:2535801/1‑93 (MQ=255)
                          tCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaa                                                                >  1:376626/1‑93 (MQ=255)
                          tCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAAGCCGCCTTTCGTCCTCaa                                                                <  1:3241789/93‑1 (MQ=255)
                           cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                                               <  1:2651659/93‑1 (MQ=255)
                           cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                                               <  2:325787/93‑1 (MQ=255)
                           cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                                               >  1:4171046/1‑93 (MQ=255)
                             gagTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACa                                                             <  2:4347613/93‑1 (MQ=255)
                             gagTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACa                                                             >  2:1673337/1‑93 (MQ=255)
                              agTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACAt                                                            >  1:4574591/1‑93 (MQ=255)
                               gTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATa                                                           <  1:2233226/93‑1 (MQ=255)
                                tttCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATaa                                                          <  1:2119122/93‑1 (MQ=255)
                                tttCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATaa                                                          <  2:1575061/93‑1 (MQ=255)
                                  tCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAtt                                                        <  1:2975539/93‑1 (MQ=255)
                                  tCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAtt                                                        <  2:3687969/93‑1 (MQ=255)
                                  tCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAtt                                                        <  1:4315427/93‑1 (MQ=255)
                                  tCGCCAGCGTAATGGTGAACCGGGTTCCCCCGAAGTTGATAAGTCTTTGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAtt                                                        <  1:189601/93‑1 (MQ=255)
                                     cgAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGt                                                     >  2:1180169/1‑93 (MQ=255)
                                       aGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTAt                                                   >  1:3214889/1‑93 (MQ=255)
                                        gCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                                                  <  1:2944223/93‑1 (MQ=255)
                                        gCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                                                  <  1:111816/93‑1 (MQ=255)
                                        gCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                                                  >  1:1553078/1‑93 (MQ=255)
                                         cGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGa                                                 <  2:3620398/93‑1 (MQ=255)
                                          gTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAg                                                <  2:182764/93‑1 (MQ=255)
                                           tATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGt                                               >  2:3142710/1‑93 (MQ=255)
                                            aTTGTCGCACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTa                                              <  2:4662691/93‑1 (MQ=255)
                                            aTTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTa                                              <  1:133245/93‑1 (MQ=255)
                                            aTTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTa                                              <  2:67923/93‑1 (MQ=255)
                                             ttGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAt                                             >  1:757546/1‑93 (MQ=255)
                                             ttGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAt                                             >  2:1353350/1‑93 (MQ=255)
                                              tGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAtt                                            <  1:4393407/93‑1 (MQ=255)
                                               gTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTa                                           >  2:2238420/1‑93 (MQ=255)
                                                tCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTaa                                          <  1:4215097/93‑1 (MQ=255)
                                                tCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTaa                                          <  1:2070698/93‑1 (MQ=255)
                                                   aaCTGGTTCCACGGGAGGTGAGTAGTCTTGGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAAccc                                       <  2:3033931/93‑1 (MQ=255)
                                                   aaCTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAAccc                                       >  1:884425/1‑93 (MQ=255)
                                                     cTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATGTGTATGAGTATTAACCcgc                                     <  1:1306307/93‑1 (MQ=255)
                                                      tGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCcgcg                                    >  2:446167/1‑93 (MQ=255)
                                                      tGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCcgcg                                    <  2:3880330/93‑1 (MQ=255)
                                                       gggTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGa                                   >  1:4104950/1‑93 (MQ=255)
                                                        ggTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGAt                                  <  1:4081438/93‑1 (MQ=255)
                                                           tGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTc                                <  2:4246440/93‑1 (MQ=255)
                                                           tGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTc                                <  1:1982018/93‑1 (MQ=255)
                                                              accGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGa                             <  2:1939651/92‑1 (MQ=255)
                                                              accGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGa                             <  2:1518064/92‑1 (MQ=255)
                                                               ccGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGAt                            >  2:285701/1‑93 (MQ=255)
                                                               ccGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGAt                            <  1:3818240/93‑1 (MQ=255)
                                                                    ggTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGgcgc                       <  1:3074083/93‑1 (MQ=255)
                                                                     gtaagtggtCTTGGGCCCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtt                     <  2:947573/90‑1 (MQ=255)
                                                                     gTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCt                      <  2:1003702/93‑1 (MQ=255)
                                                                     gTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCt                      >  1:1193642/1‑93 (MQ=255)
                                                                     gTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCt                      <  1:374793/93‑1 (MQ=255)
                                                                      tGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtt                     >  1:2366842/1‑93 (MQ=255)
                                                                      tGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtt                     <  1:3046685/93‑1 (MQ=255)
                                                                       gAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCttt                    <  1:146722/93‑1 (MQ=255)
                                                                        agtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtttt                   >  2:2720932/1‑93 (MQ=255)
                                                                        agtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtttt                   <  1:2500085/93‑1 (MQ=255)
                                                                        agtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtttt                   >  1:3123449/1‑93 (MQ=255)
                                                                         gtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTg                  >  2:3427414/1‑93 (MQ=255)
                                                                          tagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGc                 <  2:218822/93‑1 (MQ=255)
                                                                             tCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTAc              >  2:4013384/1‑93 (MQ=255)
                                                                             tCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTAc              >  2:663058/1‑93 (MQ=255)
                                                                             tCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTAc              >  1:2125648/1‑93 (MQ=255)
                                                                             tCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAACATTAACACGCGATTCTGATGGCGGTTTTGCttc              >  1:108106/1‑91 (MQ=255)
                                                                              cTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACa             <  1:108411/93‑1 (MQ=255)
                                                                                  ggCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATaa         >  2:848881/1‑93 (MQ=255)
                                                                                     aCCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAg      <  1:2636691/93‑1 (MQ=255)
                                                                                     aCCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAg      >  2:922233/1‑93 (MQ=255)
                                                                                      ccAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGc     >  1:3949976/1‑93 (MQ=255)
                                                                                      ccAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGc     >  1:3433475/1‑93 (MQ=255)
                                                                                       caacaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCg    <  1:4677141/93‑1 (MQ=255)
                                                                                        aacaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgt   <  1:654473/93‑1 (MQ=255)
                                                                                         acaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgtg  >  2:2208196/1‑92 (MQ=255)
                                                                                              |                                                                                       
TGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCGTT  >  NC_000913/144878‑145057

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