Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,757 T→C 15.4% 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,7570TC15.4% 260.8 / 37.9 118intergenic (‑140/+24)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (58/41);  new base C (14/4);  total (73/45)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-02

Only 100 of 122 total aligned reads displayed.

ATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478659‑478856
                                                                                                  |                                                                                                   
aTTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCttt                                                                                                    <  1:2174808/100‑1 (MQ=255)
     ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                               >  1:1270245/1‑100 (MQ=255)
      gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAGGTGATGAGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                              >  2:1736276/1‑100 (MQ=255)
      gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGAGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                              >  2:1737387/1‑100 (MQ=255)
      gAAAAAGGCTGGCCGGAGATAAGGAGGTAGAGCTGGTATTTCATAAAGTGATGGGTGGAAACGATGATTGTCATTTGACAAGGCGCAGAGCTTTACg                                                                                                 >  2:2167303/1‑97 (MQ=255)
       aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             <  2:297674/100‑1 (MQ=255)
       aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             <  1:1650514/100‑1 (MQ=255)
        aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                            >  1:720103/1‑100 (MQ=255)
         aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                           <  2:7909/100‑1 (MQ=38)
         aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                           <  1:1922902/100‑1 (MQ=38)
          aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                          <  1:1405542/100‑1 (MQ=255)
          aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                          <  1:2239194/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)
                 gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                   >  1:1970033/1‑100 (MQ=255)
                 gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                   >  2:1774701/1‑100 (MQ=255)
                  ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                  >  2:357280/1‑100 (MQ=255)
                  ccggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGTAACGATCATTGTCATTTGACCAGGCGCAGAGCTTGACGCCTGGTATGTATccc                                                                                  >  2:1503775/3‑100 (MQ=255)
                    ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGGGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCt                                                                                 <  2:2228009/99‑1 (MQ=255)
                     gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                               >  1:298943/1‑100 (MQ=255)
                     gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                               >  1:296546/1‑100 (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)
                                gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                    <  2:526793/100‑1 (MQ=255)
                                 tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                   >  2:1904503/1‑100 (MQ=255)
                                  agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                  <  1:1841272/100‑1 (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)
                                             ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                       <  1:140172/100‑1 (MQ=255)
                                             ttCATGAAGTGATGGGTGGAAACGATCATTATCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                       >  1:408236/1‑100 (MQ=255)
                                              tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                      >  2:1034438/1‑100 (MQ=255)
                                                 tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                   <  1:2053317/100‑1 (MQ=255)
                                                  gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                  <  1:490515/100‑1 (MQ=255)
                                                  gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                  >  1:638664/1‑100 (MQ=255)
                                                   aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                 >  2:358169/1‑100 (MQ=255)
                                                    aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                <  2:296542/100‑1 (MQ=255)
                                                    aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                >  2:71734/1‑100 (MQ=255)
                                                     gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                               >  2:2275236/1‑100 (MQ=255)
                                                     gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                               >  2:438967/1‑100 (MQ=255)
                                                       gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                             >  2:611439/1‑100 (MQ=255)
                                                       gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaa                                              >  2:202395/1‑99 (MQ=255)
                                                           ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                         >  2:1400001/1‑100 (MQ=255)
                                                            gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        >  2:1883310/1‑100 (MQ=255)
                                                            gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        >  2:35601/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)
                                                                aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                    <  1:1645368/100‑1 (MQ=255)
                                                                aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                    <  1:877857/100‑1 (MQ=255)
                                                                 aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                   <  1:1356175/100‑1 (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)
                                                                           ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                         >  1:1945375/9‑100 (MQ=25)
                                                                              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                   <  2:266393/100‑1 (MQ=255)
                                                                                  taCCAGGCGTAAAGCTCTGCGCCTGGTATGTGACCCTAACTTAAGATGTTGATTCGGCGCGAGCGGGGGTAAAAATATCGTAAGAAGACAg                           >  2:2401139/2‑91 (MQ=37)
                                                                                  taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  2:474485/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)
                                                                                   aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACATTTa                 >  1:2394655/1‑100 (MQ=255)
                                                                                    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                >  2:2400757/1‑100 (MQ=255)
                                                                                    ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                <  1:2206086/100‑1 (MQ=255)
                                                                                    ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                <  1:1857877/100‑1 (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)
                                                                                         cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTACGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           >  2:1076126/1‑100 (MQ=255)
                                                                                         cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           >  2:1662535/1‑100 (MQ=255)
                                                                                         cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           <  1:907704/100‑1 (MQ=255)
                                                                                         cgaagagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTACAAAGACAGTTCACCTTTACGGACg           >  1:1559388/4‑100 (MQ=255)
                                                                                               gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt     >  1:2031504/1‑100 (MQ=255)
                                                                                               gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt     >  1:575577/1‑100 (MQ=255)
                                                                                                 tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGCTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGGTCAgtt    <  2:2174892/99‑1 (MQ=255)
                                                                                                  ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  >  1:1293306/1‑100 (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: ‑