Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,750 C→T 15.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,7500CT15.0% 427.4 / 51.5 172intergenic (‑133/+31)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (72/73);  new base T (18/8);  total (91/81)
Fisher's exact test for biased strand distribution p-value = 8.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 176 total aligned reads displayed.

ATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAT  >  NC_000913/478652‑478849
                                                                                                  |                                                                                                   
atatCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCa                                                                                                    >  1:1820927/1‑100 (MQ=255)
   tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCagag                                                                                                 >  2:3216949/1‑100 (MQ=255)
   tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCagag                                                                                                 >  1:2625511/1‑100 (MQ=255)
    cTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                                >  2:2793171/1‑100 (MQ=255)
      gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                                >  1:1339457/1‑98 (MQ=255)
      gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCtt                                                                                              >  1:1475951/1‑100 (MQ=255)
          gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                          >  2:1004346/1‑100 (MQ=255)
            ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAATGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        <  2:1942709/100‑1 (MQ=255)
            ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        <  1:3580013/100‑1 (MQ=255)
            ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        <  2:915925/100‑1 (MQ=255)
             gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                       >  2:18805/1‑100 (MQ=255)
              aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                      <  2:965241/100‑1 (MQ=255)
              aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                      <  2:2405751/100‑1 (MQ=255)
              aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                      <  1:3479700/100‑1 (MQ=255)
               aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                     <  1:2760696/100‑1 (MQ=255)
               aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                     >  2:360416/1‑100 (MQ=255)
                aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                     <  1:1984090/99‑1 (MQ=39)
                aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                    <  1:3334364/100‑1 (MQ=38)
                aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                    <  1:2804808/100‑1 (MQ=38)
                 aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                   <  1:2103221/100‑2 (MQ=38)
                 aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                   <  1:1169301/100‑2 (MQ=38)
                 aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                   <  2:1451889/100‑1 (MQ=255)
                  aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                  <  2:2690056/100‑3 (MQ=38)
                   ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                 >  2:2314564/1‑100 (MQ=255)
                    gctggctgGAGATGAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGt                                                                                >  1:1108657/1‑100 (MQ=255)
                    gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGt                                                                                >  2:2777570/1‑100 (MQ=255)
                       ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             >  2:3308877/1‑100 (MQ=255)
                       ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             <  1:2270104/100‑1 (MQ=255)
                        gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaatgac                                                                            <  2:2611313/100‑9 (MQ=25)
                        gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTCAcc                                                                            <  2:1764988/100‑1 (MQ=255)
                        gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                            >  2:2346180/1‑100 (MQ=255)
                          tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCt                                                                          <  2:1475943/100‑1 (MQ=255)
                             agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                       >  1:1582648/1‑100 (MQ=255)
                             agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                       <  1:1951279/100‑1 (MQ=255)
                              gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACt                                                                      >  1:2422206/1‑100 (MQ=255)
                                tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                    <  1:1757514/100‑1 (MQ=255)
                                tAAGGACGTACAGCTGTTATTTCATGAAGTGATGGGTGGAAACGATCATTTTCTTTTGACCATGCGCATAGCTTTACGCCTTGTATGTAACCCTAACt                                                                      <  2:424317/98‑1 (MQ=255)
                                 aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                   <  1:702502/100‑1 (MQ=255)
                                 aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                   <  1:2777588/100‑1 (MQ=255)
                                  aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGATCAGGCGCAGAGCTTTACGCCTGGTATGTCACCCTAACTTatg                                                                  <  1:2467384/100‑3 (MQ=255)
                                  aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAg                                                                  >  1:73914/1‑100 (MQ=255)
                                   ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAACCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                 <  2:174954/100‑1 (MQ=255)
                                   ggaggTAGAGCTGGTATTGCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGAGCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                 >  2:2174247/1‑100 (MQ=255)
                                    gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGAt                                                                >  2:2853803/1‑100 (MQ=255)
                                        ttgagCTTGTAATTCATGAAGTGACGGATCGAAACGATCATTGTCAGGTGGCCAGGAGTAGACCTATACGCCTGGTAGGTAACGCTTACTGTAGATGTTg                                                            >  2:1498675/3‑100 (MQ=255)
                                              tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                      >  2:43158/1‑100 (MQ=255)
                                              tGGTATTTCATGAAGTGACGAGTGGGAAGGATCATTGTCATTTGACCAGACAAAGCGCTTTCCGCCTGGTATGTAAGCCTAACTTAAGACGTTGATTc                                                        >  2:2340201/1‑98 (MQ=255)
                                               ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGACATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                     <  1:2563369/100‑1 (MQ=255)
                                                gTATTTCATGAAGTGATGGGTGGAAATGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                    <  2:2033205/100‑1 (MQ=255)
                                                gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                    <  1:603681/100‑1 (MQ=255)
                                                  aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                  >  2:3053975/1‑100 (MQ=255)
                                                     tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                               <  1:3117007/100‑1 (MQ=255)
                                                     tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                               >  1:349286/1‑100 (MQ=255)
                                                        tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                            <  1:2326383/100‑1 (MQ=255)
                                                         gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                           <  1:3044222/100‑1 (MQ=255)
                                                               aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                     >  2:1205705/1‑100 (MQ=255)
                                                                tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                    <  1:3075956/100‑1 (MQ=255)
                                                                  ggTGGAAACGATCATTGTCATTTGACTAGGCGCAGAGCTTTACGCTTGGTATGTAACTCTAACTTAAGATTTTGATTCGGCTCGAGCGTGGGTAataata                                  >  2:3270525/1‑100 (MQ=255)
                                                                   gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTCACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                 >  2:1603271/1‑100 (MQ=255)
                                                                   gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                 <  2:1135800/100‑1 (MQ=255)
                                                                    tagAAACGCTCATTGTCATTTGACCAGGAGCAGAGCTTTACGCCTGGTATGTAACCCTCACTTAAGATGTTGATTCGGCGCGAGAGAGGGTAAAAATATc                                <  2:2148608/98‑1 (MQ=255)
                                                                    tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                <  2:1616343/100‑1 (MQ=255)
                                                                      gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                              <  1:3290890/100‑1 (MQ=255)
                                                                      gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                              >  1:3260734/1‑100 (MQ=255)
                                                                      gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                              <  1:393888/100‑1 (MQ=255)
                                                                       aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                             >  1:1025712/1‑100 (MQ=255)
                                                                       aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                             >  1:309423/1‑100 (MQ=255)
                                                                         aCGATAATTGTCATTTGACCAAGCGCAGAGTTTTACGCCTGGGATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                               <  2:3127854/96‑1 (MQ=255)
                                                                           gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGGTGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                         <  2:1562261/100‑1 (MQ=255)
                                                                           gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                         <  1:3325887/100‑1 (MQ=255)
                                                                               aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                     >  1:2193198/1‑100 (MQ=255)
                                                                               aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                     >  1:2697995/1‑100 (MQ=255)
                                                                                ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                    <  1:367707/100‑1 (MQ=255)
                                                                                ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                    <  2:2422190/100‑1 (MQ=255)
                                                                                  ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                    >  2:1716853/9‑98 (MQ=25)
                                                                                  ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                  >  1:2077089/9‑100 (MQ=25)
                                                                                  gTCATTTGACCAGGCGCAGAGCTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                 <  2:3089771/100‑1 (MQ=255)
                                                                                    cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                >  2:3300364/1‑100 (MQ=255)
                                                                                    cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                >  1:1988287/1‑100 (MQ=255)
                                                                                    cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTAATTCGGAGCGAGCGAGGGTAAAAATATCGTAAAACGACAGTTCa                <  2:1259525/100‑1 (MQ=255)
                                                                                     tacataCCAGGCGTAAAGCTCTGCGCCTGGTATATAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAAGGTAAAAATATCGTAAAATGACAGTTccc               >  2:2185510/6‑98 (MQ=25)
                                                                                      tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc              >  1:1044480/1‑100 (MQ=255)
                                                                                       cataCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt             >  1:1840294/4‑100 (MQ=255)
                                                                                        tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt            >  1:90477/1‑100 (MQ=255)
                                                                                        ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt            >  1:761863/3‑100 (MQ=38)
                                                                                        ataCCAGGCGTAAAGCTCTGCGCATGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt            >  2:732771/3‑100 (MQ=255)
                                                                                         taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAAGCCTAACGTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  2:301960/2‑100 (MQ=37)
                                                                                         gACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  2:2590721/1‑100 (MQ=255)
                                                                                          aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          <  1:3290819/100‑1 (MQ=255)
                                                                                          aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          <  1:1379453/100‑1 (MQ=255)
                                                                                          aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          >  2:3583120/1‑100 (MQ=255)
                                                                                          aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          >  2:209301/1‑100 (MQ=255)
                                                                                          aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          <  2:2411404/100‑1 (MQ=255)
                                                                                           ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  1:404729/1‑100 (MQ=255)
                                                                                           ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  2:16499/1‑100 (MQ=255)
                                                                                           ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCTAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         <  2:1988276/100‑1 (MQ=255)
                                                                                           ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  2:2513331/1‑100 (MQ=255)
                                                                                            cAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg        >  1:2159836/1‑100 (MQ=255)
                                                                                               gcacAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc     >  2:18667/4‑100 (MQ=255)
                                                                                                  cAGAGCTTTACGCCTGGTATGTAACCCGAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt  <  1:800056/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAT  >  NC_000913/478652‑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: