Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,757 T→C 16.7% intergenic (‑140/+24) 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,7570TC16.7% 267.6 / 42.8 120intergenic (‑140/+24)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (51/48);  new base C (11/9);  total (63/57)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 125 total aligned reads displayed.

ATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478659‑478856
                                                                                                  |                                                                                                   
aTTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCttt                                                                                                    >  1:1106403/1‑100 (MQ=255)
 ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                                   <  1:2564764/100‑1 (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)
     ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTAGACCAGGCGCAGAGCTTTACGcc                                                                                               >  2:1131814/1‑100 (MQ=255)
      gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                              >  2:1768215/1‑100 (MQ=255)
        aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                            <  2:701472/100‑1 (MQ=255)
         aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATTGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                            <  1:1211323/99‑1 (MQ=38)
         aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                           <  1:1497258/100‑1 (MQ=38)
         aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                           <  1:1053779/100‑1 (MQ=39)
          aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                          <  2:221578/100‑2 (MQ=38)
          aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                          <  1:1636080/100‑1 (MQ=255)
           aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                         <  2:1730178/100‑3 (MQ=38)
           aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                         <  2:2441421/100‑3 (MQ=38)
            ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                        >  2:1084177/1‑100 (MQ=255)
              ctggctggAGATAAGGAGGTAGAGCTTGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaa                                                                                       <  1:1675746/99‑5 (MQ=37)
              ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                      <  1:47689/100‑1 (MQ=255)
               tggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGGGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATCTa                                                                                      <  1:2591394/99‑1 (MQ=255)
                ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                    <  2:1076404/100‑1 (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)
                     gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                               >  1:336594/1‑100 (MQ=255)
                      agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                              >  2:2335902/1‑100 (MQ=255)
                      agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                              <  2:1816420/100‑1 (MQ=255)
                       gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACAAGGCGCAGAGCTGTACGCCAGGTATGTAACCCTAAc                                                                              >  1:2335897/1‑99 (MQ=255)
                          aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGACCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                          <  2:443978/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)
                               ggTAGAGCTGTTATTTCATGAAGTGATGGGTCGAAACGATCATTGTCATTTTACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATg                                                                      <  1:1891195/99‑1 (MQ=255)
                                 tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                    <  2:2623056/99‑1 (MQ=255)
                                  agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCCTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                  >  2:1423743/1‑100 (MQ=255)
                                   gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                 >  2:1221009/1‑100 (MQ=255)
                                     gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                               >  1:1281219/1‑100 (MQ=255)
                                     gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                               >  2:1342733/1‑100 (MQ=255)
                                      cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCg                                                              >  2:2750339/1‑100 (MQ=255)
                                       tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                             <  2:2173182/100‑1 (MQ=255)
                                       tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                             <  2:458889/100‑1 (MQ=255)
                                         gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                           <  2:1892793/100‑1 (MQ=255)
                                          tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                          >  2:1223158/1‑100 (MQ=255)
                                           aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         >  1:2774945/1‑100 (MQ=255)
                                           aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         >  2:23603/1‑100 (MQ=255)
                                           aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         >  2:701079/1‑100 (MQ=255)
                                             ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGTCCAGGCGCAGAGCTTTACGCCTAGTATGTGACCCTAACTTAAGATGTTGCTTCGgcg                                                           >  1:1794304/1‑96 (MQ=255)
                                              tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                      >  2:617653/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)
                                                           ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAACata                                         >  2:1734109/1‑100 (MQ=255)
                                                            gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        <  1:913015/100‑1 (MQ=255)
                                                            gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        >  1:2342114/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)
                                                                   cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAACAATATCGTaaaa                                 >  1:2665896/1‑100 (MQ=255)
                                                                       cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                             <  2:1982654/100‑1 (MQ=255)
                                                                          tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                          <  1:1784469/100‑1 (MQ=255)
                                                                          tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                          >  2:2647691/1‑100 (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)
                                                                                  taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  2:1723006/2‑100 (MQ=38)
                                                                                  gACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  <  1:2153338/100‑1 (MQ=255)
                                                                                   aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                     <  2:2612555/96‑1 (MQ=255)
                                                                                   aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 >  1:497727/1‑100 (MQ=255)
                                                                                    ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  2:1633094/1‑100 (MQ=255)
                                                                                    ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  1:819154/1‑100 (MQ=255)
                                                                                    ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  2:2742741/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)
                                                                                     cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGACAGTTCACCTTTACg               >  1:2357468/1‑100 (MQ=255)
                                                                                     cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTAACCTTTACg               >  1:566712/1‑100 (MQ=255)
                                                                                      aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg              <  1:1871425/100‑1 (MQ=255)
                                                                                      aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg              <  2:1923918/100‑1 (MQ=255)
                                                                                         cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           >  2:2697492/1‑100 (MQ=255)
                                                                                          gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTATAAAGACAGTTCACCTTTACGGACg           >  2:1105471/1‑99 (MQ=255)
                                                                                           cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt         >  1:2281361/1‑100 (MQ=255)
                                                                                           cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt         <  2:153194/100‑1 (MQ=255)
                                                                                                  ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  >  1:742806/1‑100 (MQ=255)
                                                                                                  ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  <  2:1489799/100‑1 (MQ=255)
                                                                                                  ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  <  2:2193241/100‑1 (MQ=255)
                                                                                                  ttACGCCTGGTATGTAACCCTAAATTAAGATGTTGATTCGGCGCGAGCGATGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  <  1:896497/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478659‑478856

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: