Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,750 C→T 16.0% intergenic (‑133/+31) 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,7500CT16.0% 297.5 / 38.2 124intergenic (‑133/+31)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (51/52);  new base T (11/9);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 130 total aligned reads displayed.

TATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAT  >  NC_000913/478651‑478849
                                                                                                   |                                                                                                   
tataTCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGgcgc                                                                                                     <  1:459550/100‑1 (MQ=255)
 atatCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCa                                                                                                    >  2:1477623/1‑100 (MQ=255)
   atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCaga                                                                                                  <  2:1349349/100‑1 (MQ=255)
   atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCaga                                                                                                  <  1:2228063/100‑1 (MQ=255)
    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCagag                                                                                                 >  2:2404240/1‑100 (MQ=255)
     cTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                                <  1:304759/100‑1 (MQ=255)
     cTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                                <  1:2755204/100‑1 (MQ=255)
     cTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTAGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                                <  1:782340/100‑1 (MQ=255)
        aTTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCttt                                                                                             >  1:1106403/1‑100 (MQ=255)
         ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                            <  1:2564764/100‑1 (MQ=255)
         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:2254809/1‑100 (MQ=255)
             ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        >  2:51833/1‑100 (MQ=255)
             ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        <  1:1511501/100‑1 (MQ=255)
             ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTAGACCAGGCGCAGAGCTTTACGcc                                                                                        >  2:1131814/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)
                  aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                   <  2:221578/100‑2 (MQ=38)
                  aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                   <  1:1636080/100‑1 (MQ=255)
                   aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                  <  2:2441421/100‑3 (MQ=38)
                      ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                               <  1:47689/100‑1 (MQ=255)
                        ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             <  2:1076404/100‑1 (MQ=255)
                         gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             >  1:1728310/1‑99 (MQ=255)
                         gctggAGATAAGGAGGTAGAGCTGGTATTCCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCGTTACGCCTGGTATGTAAcc                                                                            >  2:2153345/1‑100 (MQ=255)
                          ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                           <  1:694910/100‑1 (MQ=255)
                          ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                           <  1:51832/100‑1 (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:336594/1‑100 (MQ=255)
                             gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                        >  1:2581610/1‑100 (MQ=255)
                              agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                       <  2:1816420/100‑1 (MQ=255)
                               gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACAAGGCGCAGAGCTGTACGCCAGGTATGTAACCCTAAc                                                                       >  1:2335897/1‑99 (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)
                                          agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCCTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                           >  2:1423743/1‑100 (MQ=255)
                                           gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                          >  2:1221009/1‑100 (MQ=255)
                                             gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                        >  2:1342733/1‑100 (MQ=255)
                                             gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                        >  1:1281219/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:1477070/100‑1 (MQ=255)
                                                 gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                    <  2:2670037/100‑1 (MQ=255)
                                                  tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                   >  2:1223158/1‑100 (MQ=255)
                                                   aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                  >  2:701079/1‑100 (MQ=255)
                                                   aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                  >  1:2774945/1‑100 (MQ=255)
                                                   aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                  >  2:23603/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)
                                                         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:2342114/1‑100 (MQ=255)
                                                                    gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                 <  1:913015/100‑1 (MQ=255)
                                                                        aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                               <  1:2890091/98‑1 (MQ=255)
                                                                        aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                             <  2:2686283/100‑1 (MQ=255)
                                                                           cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTCAAAATATCGTaaaa                          >  1:1982646/1‑100 (MQ=255)
                                                                           cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAACAATATCGTaaaa                          >  1:2665896/1‑100 (MQ=255)
                                                                           cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                          >  2:1050823/1‑100 (MQ=255)
                                                                               cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                      <  2:1982654/100‑1 (MQ=255)
                                                                                aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCGAACTTAAGATGTTGTTTTGGCCCGACCGAGGGTAAAAATATCAtaataaaa                       >  1:792430/1‑96 (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)
                                                                                         ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACATGTCACCtt            >  1:533768/3‑100 (MQ=37)
                                                                                          taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  2:1723006/2‑100 (MQ=38)
                                                                                           aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCTGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          >  2:2787770/1‑100 (MQ=38)
                                                                                           aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc              <  2:2612555/96‑1 (MQ=255)
                                                                                            ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  2:2449899/1‑100 (MQ=255)
                                                                                            ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  1:717042/1‑100 (MQ=255)
                                                                                            ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  2:2742741/1‑100 (MQ=255)
                                                                                            ccAGGCGTAAAGCTCTGCGCCTGGTATGTAAACCTAACTTAAGATGTTGATTAGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt            >  1:313352/1‑97 (MQ=38)
                                                                                            ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         <  2:2502753/100‑1 (MQ=255)
                                                                                            ccACGCGTAAAGCTCTGCGCCTGGTATGTAACCATAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATAGCGTAAACCGACAGTCCACCttaag         >  1:265378/1‑97 (MQ=255)
                                                                                             cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGACAGTTCACCTTTACg        >  1:2357468/1‑100 (MQ=255)
                                                                                             cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTAACCTTTACg        >  1:566712/1‑100 (MQ=255)
                                                                                              aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg       <  2:1923918/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)
                                                                                                   cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt  <  2:153194/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAT  >  NC_000913/478651‑478849

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: