Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,759 A→G 16.0% intergenic (‑142/+22) ylaB ← / ← ylaC putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7590AG16.0% 278.8 / 42.0 124intergenic (‑142/+22)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (53/51);  new base G (11/9);  total (64/60)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 129 total aligned reads displayed.

TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478660‑478858
                                                                                                   |                                                                                                   
ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                                     >  1:1923910/1‑100 (MQ=255)
  gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                   >  1:75345/1‑100 (MQ=255)
  gtgGAAAAAGGCTGGCTGGAGATAAGGAGGGAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                   >  2:349464/1‑100 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTTGTATTTCATGACGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 <  1:546685/100‑1 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGcc                                                                                                 <  1:2899624/100‑1 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 >  2:51833/1‑100 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 >  2:2254809/1‑100 (MQ=255)
     gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                >  2:1768215/1‑100 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgn                                                                                              <  2:1773732/100‑2 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                              <  2:701472/100‑1 (MQ=255)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATTGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                              <  1:1211323/99‑1 (MQ=38)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                             <  1:1053779/100‑1 (MQ=39)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                             <  1:1497258/100‑1 (MQ=38)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                            <  1:1636080/100‑1 (MQ=255)
          aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                           <  2:2441421/100‑3 (MQ=38)
          aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                           <  2:1730178/100‑3 (MQ=38)
           ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                          >  2:1084177/1‑100 (MQ=255)
             ctggctggAGATAAGGAGGTAGAGCTTGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaa                                                                                         <  1:1675746/99‑5 (MQ=37)
              tggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGGGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATCTa                                                                                        <  1:2591394/99‑1 (MQ=255)
               ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                      <  2:1076404/100‑1 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTCCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCGTTACGCCTGGTATGTAAcc                                                                                     >  2:2153345/1‑100 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    <  1:694910/100‑1 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    <  1:51832/100‑1 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    >  1:30392/1‑100 (MQ=255)
                  tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCt                                                                                   >  2:991125/1‑100 (MQ=255)
                   ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                  <  2:2343984/100‑1 (MQ=255)
                   gcgcaTAAGGGCGTAGAGCTGGTATGTCATGAAGTGATGGGCGGAAACGATCATTGTCATTTGACGAGGCGGAGTGCTTTACTCTTGATCTGTAACCCt                                                                                   >  1:867963/5‑99 (MQ=255)
                    gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGCCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                  <  1:2562539/99‑1 (MQ=255)
                    gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                                 >  1:2581610/1‑100 (MQ=255)
                     agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                                >  2:2335902/1‑100 (MQ=255)
                      gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACAAGGCGCAGAGCTGTACGCCAGGTATGTAACCCTAAc                                                                                >  1:2335897/1‑99 (MQ=255)
                        tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                             <  1:2330074/100‑1 (MQ=255)
                          aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAg                                                                           >  2:1353700/1‑100 (MQ=255)
                           ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                          <  1:1471192/100‑1 (MQ=255)
                            gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGAt                                                                         <  1:73920/100‑1 (MQ=255)
                             aggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATg                                                                        <  1:1353698/100‑1 (MQ=255)
                                  gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   >  2:1221009/1‑100 (MQ=255)
                                    gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                 >  2:1342733/1‑100 (MQ=255)
                                     cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCg                                                                >  2:2750339/1‑100 (MQ=255)
                                      tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                               <  2:458889/100‑1 (MQ=255)
                                      tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                               <  2:2173182/100‑1 (MQ=255)
                                        gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                             <  2:2670037/100‑1 (MQ=255)
                                        gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                             <  2:1477070/100‑1 (MQ=255)
                                        gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                             <  2:1892793/100‑1 (MQ=255)
                                          aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                           >  2:701079/1‑100 (MQ=255)
                                          aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                           >  1:2774945/1‑100 (MQ=255)
                                            ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGTCCAGGCGCAGAGCTTTACGCCTAGTATGTGACCCTAACTTAAGATGTTGCTTCGgcg                                                             >  1:1794304/1‑96 (MQ=255)
                                              cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                       >  2:1706604/1‑100 (MQ=255)
                                               aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                      >  2:476173/1‑100 (MQ=255)
                                               aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCATAACTTAAGATGTTGATTCGGCgcgagcga                                                      >  1:2691961/1‑100 (MQ=255)
                                                tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                     >  1:384359/1‑100 (MQ=255)
                                                 gAAGTGATGGGTGGAAACGATCATTGTCATTTGGCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                    <  1:2675910/100‑1 (MQ=255)
                                                  aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                   <  1:1641302/100‑1 (MQ=255)
                                                   aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                  <  2:1247838/100‑1 (MQ=255)
                                                       aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGTTTCGGCGCGAGCGAGGGTaaaa                                              <  1:1686550/100‑1 (MQ=255)
                                                         gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                            <  2:1289298/100‑1 (MQ=255)
                                                          ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAACata                                           >  2:1734109/1‑100 (MQ=255)
                                                          ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                           >  1:2514994/1‑100 (MQ=255)
                                                            tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                         <  1:743664/100‑1 (MQ=255)
                                                             ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        <  2:2166815/100‑1 (MQ=255)
                                                             ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        <  2:1437139/100‑1 (MQ=255)
                                                               aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                      <  2:2686283/100‑1 (MQ=255)
                                                                  cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTCAAAATATCGTaaaa                                   >  1:1982646/1‑100 (MQ=255)
                                                                  cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                   >  2:1050823/1‑100 (MQ=255)
                                                                     tCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                                >  1:121535/1‑100 (MQ=255)
                                                                      cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                               <  2:1982654/100‑1 (MQ=255)
                                                                        ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             <  1:379422/100‑1 (MQ=255)
                                                                         tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                            >  2:2647691/1‑100 (MQ=255)
                                                                         tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                            <  1:1784469/100‑1 (MQ=255)
                                                                              tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                       <  1:2180950/100‑1 (MQ=255)
                                                                               ttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                      <  2:1691915/100‑1 (MQ=255)
                                                                                tGACCAGGCGCAGAGCTTTATGCCTGGTATGTACCACTCACATACGATGTTGACTAAGCGCGAGGGACGCTAAAAAGAACGTaaaca                                  >  1:2344364/1‑85 (MQ=255)
                                                                                tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:1205092/1‑100 (MQ=255)
                                                                                ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGGGCGAGGGTAAAAATATCGTAAAAAGACAGtgc                          >  1:583960/3‑93 (MQ=37)
                                                                                ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACATGTCACCtt                     >  1:533768/3‑100 (MQ=37)
                                                                                 gACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                    <  1:2153338/100‑1 (MQ=255)
                                                                                  aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCTGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  2:2787770/1‑100 (MQ=38)
                                                                                  aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                       <  2:2612555/96‑1 (MQ=255)
                                                                                  aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  1:497727/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  1:819154/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:2742741/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:1633094/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:2449899/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAAACCTAACTTAAGATGTTGATTAGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:313352/1‑97 (MQ=38)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAATACAGTTCACCTTTAc                  <  2:1186413/100‑1 (MQ=255)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  <  2:2502753/100‑1 (MQ=255)
                                                                                   ccACGCGTAAAGCTCTGCGCCTGGTATGTAACCATAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATAGCGTAAACCGACAGTCCACCttaag                  >  1:265378/1‑97 (MQ=255)
                                                                                    cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGACAGTTCACCTTTACg                 >  1:2357468/1‑100 (MQ=255)
                                                                                     aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg                <  2:1923918/100‑1 (MQ=255)
                                                                                     aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg                <  1:1871425/100‑1 (MQ=255)
                                                                                        cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg             >  2:1168158/1‑100 (MQ=255)
                                                                                         gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTATAAAGACAGTTCACCTTTACGGACg             >  2:1105471/1‑99 (MQ=255)
                                                                                          cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt           >  1:2281361/1‑100 (MQ=255)
                                                                                                 ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt    <  2:1489799/100‑1 (MQ=255)
                                                                                                 ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt    <  2:2193241/100‑1 (MQ=255)
                                                                                                 ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt    >  1:742806/1‑100 (MQ=255)
                                                                                                 ttACGCCTGGTATGTAACCCTAAATTAAGATGTTGATTCGGCGCGAGCGATGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt    <  1:896497/100‑1 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  >  1:712506/1‑100 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  1:386203/100‑1 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTAATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  1:2811089/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478660‑478858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: