Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,759 A→G 16.4% 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.4% 270.2 / 43.7 122intergenic (‑142/+22)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (60/42);  new base G (16/4);  total (76/46)
Fisher's exact test for biased strand distribution p-value = 8.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 123 total aligned reads displayed.

TGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478661‑478858
                                                                                                  |                                                                                                   
   ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 >  1:1270245/1‑100 (MQ=255)
   ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 >  2:220747/1‑100 (MQ=255)
    gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAGGTGATGAGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                >  2:1736276/1‑100 (MQ=255)
    gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGAGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                >  2:1737387/1‑100 (MQ=255)
     aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                               <  2:297674/100‑1 (MQ=255)
     aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                               <  1:1650514/100‑1 (MQ=255)
     aaaaaGGCTGGCTGGAGATAAGGAGGGAGAGGTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                               <  2:2524159/100‑1 (MQ=255)
      aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                              >  1:720103/1‑100 (MQ=255)
       aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                             <  1:1922902/100‑1 (MQ=38)
       aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                             <  2:7909/100‑1 (MQ=38)
        aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                            <  1:2239194/100‑1 (MQ=255)
        aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                            <  1:1405542/100‑1 (MQ=255)
        aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGGATTTCAGGACGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  2:1797038/100‑2 (MQ=37)
         aGGGTGGCTGGAGATAAGGATGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATGTGACAAGGCGCAGAGCTTTACGCCTGGTAt                                                                                           <  2:1154596/100‑1 (MQ=255)
           gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGt                                                                                         >  2:919840/1‑100 (MQ=255)
              ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                      <  2:1222957/100‑1 (MQ=255)
               gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                     >  1:1970033/1‑100 (MQ=255)
                ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    >  2:357280/1‑100 (MQ=255)
                ccggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGTAACGATCATTGTCATTTGACCAGGCGCAGAGCTTGACGCCTGGTATGTATccc                                                                                    >  2:1503775/3‑100 (MQ=255)
                   gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                                 >  1:296546/1‑100 (MQ=255)
                        aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                            <  1:145410/100‑1 (MQ=255)
                           gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGTAAACGATCATTGTCATTTGCCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGAt                                                                         <  2:494786/100‑1 (MQ=255)
                            aggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACATAAGATg                                                                        <  1:1348366/100‑1 (MQ=255)
                            aggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCACTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATg                                                                        <  2:1635943/100‑1 (MQ=255)
                             ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCACTGTCATTTGACCAGACGCAGAAATTTACGCCTGGCATGTAACGCTAAc                                                                                >  2:2199427/1‑91 (MQ=255)
                                agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                    >  2:2103368/1‑100 (MQ=255)
                                 gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   <  2:1117769/100‑1 (MQ=255)
                                 gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   >  1:2503329/1‑100 (MQ=255)
                                   gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                 >  2:439919/1‑100 (MQ=255)
                                      ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                              >  1:2263016/1‑100 (MQ=255)
                                      ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                              >  1:1864580/1‑100 (MQ=255)
                                      ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGAGGTTGATTCGgc                                                              >  2:1250322/1‑100 (MQ=255)
                                      ggTATTGCATGAAGGGATGGGTGGGAATGATCATTGTCATTTTACCTGGCGCTGTGCTTTACGCGTGTTATTTAATCGTAACTTAAGATGTTGAtttg                                                                >  2:2172079/1‑96 (MQ=255)
                                       gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                             <  2:1480901/100‑1 (MQ=255)
                                         aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                           >  2:830265/1‑100 (MQ=255)
                                           ttCATGAAGTGATGGGTGGAAACGATCATTATCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                         >  1:408236/1‑100 (MQ=255)
                                            tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                        >  2:1034438/1‑100 (MQ=255)
                                             cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagtg                                                       <  1:2239655/100‑3 (MQ=255)
                                               tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                     <  1:2053317/100‑1 (MQ=255)
                                                gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                    >  2:1521740/1‑100 (MQ=255)
                                                gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                    >  1:638664/1‑100 (MQ=255)
                                                gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                    <  1:490515/100‑1 (MQ=255)
                                                 aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                   >  2:358169/1‑100 (MQ=255)
                                                 aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                   <  1:1657853/100‑1 (MQ=255)
                                                  aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                  <  2:296542/100‑1 (MQ=255)
                                                  aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                  <  2:298939/100‑1 (MQ=255)
                                                  aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                  >  2:71734/1‑100 (MQ=255)
                                                   gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                 >  2:438967/1‑100 (MQ=255)
                                                   gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                 >  2:2275236/1‑100 (MQ=255)
                                                   gTGATGGGTGGAAACGACCATTGCCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTTACTCAAAATGTTGATTCGGTGCGATAGAggg                                                   >  2:101552/1‑98 (MQ=255)
                                                     gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                               >  1:475757/1‑100 (MQ=255)
                                                     gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                               >  2:611439/1‑100 (MQ=255)
                                                     gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaa                                                >  2:202395/1‑99 (MQ=255)
                                                          gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                          >  2:1883310/1‑100 (MQ=255)
                                                            ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        <  1:2001307/100‑1 (MQ=255)
                                                            ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        <  2:1696648/100‑1 (MQ=255)
                                                              aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        <  2:467774/98‑1 (MQ=255)
                                                               aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                     <  1:1356175/100‑1 (MQ=255)
                                                                  gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                  >  2:52310/1‑100 (MQ=255)
                                                                  gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                  >  2:186907/1‑100 (MQ=255)
                                                                   ataaTTGTCATTTAAACAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAAAAAAATAAAAAg                                 <  2:5298/97‑1 (MQ=255)
                                                                     cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                               >  1:989683/1‑100 (MQ=255)
                                                                       ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             >  1:571799/1‑100 (MQ=255)
                                                                        tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                            <  1:1052050/100‑1 (MQ=255)
                                                                        tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGGATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                            <  2:512203/100‑1 (MQ=255)
                                                                         ntCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                           >  1:189370/2‑100 (MQ=255)
                                                                            aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                        >  1:865931/1‑100 (MQ=255)
                                                                             tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                       <  1:1176896/100‑1 (MQ=255)
                                                                             tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                       <  1:1176981/100‑1 (MQ=255)
                                                                              cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAgatcggaa                      >  2:1301532/4‑93 (MQ=37)
                                                                               tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:1379137/1‑100 (MQ=255)
                                                                               tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:684925/1‑100 (MQ=255)
                                                                               ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  2:1749626/3‑100 (MQ=38)
                                                                                taCCAGGCGTAAAGCTCTGCGCCTGGTATGTGACCCTAACTTAAGATGTTGATTCGGCGCGAGCGGGGGTAAAAATATCGTAAGAAGACAg                             >  2:2401139/2‑91 (MQ=37)
                                                                                taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                    >  2:474485/2‑100 (MQ=38)
                                                                                taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                    >  1:267800/2‑100 (MQ=38)
                                                                                 aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTTAAAAGACAGTTCACCTTTa                   >  2:1843256/1‑100 (MQ=38)
                                                                                 aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  1:1706733/1‑100 (MQ=38)
                                                                                 aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  1:1809691/1‑100 (MQ=39)
                                                                                 aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  1:424732/1‑100 (MQ=38)
                                                                                  ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:1934172/1‑100 (MQ=255)
                                                                                  ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:958005/1‑100 (MQ=255)
                                                                                  ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:269819/1‑100 (MQ=255)
                                                                                  ccAGGCGTAAAGCTATGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGACAGTTCACCttt                    >  2:1812778/1‑98 (MQ=255)
                                                                                  ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  <  1:1857877/100‑1 (MQ=255)
                                                                                  ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:2400757/1‑100 (MQ=255)
                                                                                   cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                 >  1:277945/1‑100 (MQ=255)
                                                                                   cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                 <  2:2476845/100‑1 (MQ=255)
                                                                                     ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa               >  2:2056502/1‑100 (MQ=255)
                                                                                     ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa               >  2:667049/1‑100 (MQ=255)
                                                                                     ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa               <  1:2389127/100‑1 (MQ=255)
                                                                                      gcgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc              >  1:1491734/1‑100 (MQ=255)
                                                                                       cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGagg             >  1:2014235/1‑98 (MQ=255)
                                                                                       cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg             <  1:907704/100‑1 (MQ=255)
                                                                                       cgaagagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTACAAAGACAGTTCACCTTTACGGACg             >  1:1559388/4‑100 (MQ=255)
                                                                                             gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt       >  1:2031504/1‑100 (MQ=255)
                                                                                               tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGCTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGGTCAgtt      <  2:2174892/99‑1 (MQ=255)
                                                                                                ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt    >  1:1293603/1‑100 (MQ=255)
                                                                                                 tACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt   >  1:1236249/1‑100 (MQ=255)
                                                                                                  aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  1:1428914/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478661‑478858

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