Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,759 A→G 17.7% 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,7590AG17.7% 325.8 / 61.6 147intergenic (‑142/+22)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (64/57);  new base G (18/8);  total (82/65)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 155 total aligned reads displayed.

TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478660‑478858
                                                                                                   |                                                                                                   
ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                                     >  1:2765365/1‑100 (MQ=255)
 tgtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTAc                                                                                                    >  1:3371566/1‑100 (MQ=255)
  gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                   >  1:2295643/1‑100 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCGTTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 >  2:3666225/1‑100 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 <  1:682814/100‑1 (MQ=255)
      aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                               <  1:1355025/100‑1 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                              >  2:847384/1‑100 (MQ=255)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                             <  2:3445662/100‑1 (MQ=39)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                             <  1:1343030/100‑1 (MQ=255)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  2:1080424/100‑2 (MQ=38)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  1:2953355/100‑2 (MQ=38)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                            <  1:48510/100‑1 (MQ=255)
             ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                        <  1:1131979/100‑1 (MQ=255)
                 ttggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    <  2:1022979/99‑1 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    >  1:1697859/1‑100 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    <  2:1001785/100‑1 (MQ=255)
                  tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTGATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccca                                                                                   >  1:1022988/1‑99 (MQ=255)
                  tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCt                                                                                   >  1:2823276/1‑100 (MQ=255)
                   ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                  >  2:1075570/1‑100 (MQ=255)
                       aTAAGGAGGTAGAGCTGGTATTTCATGCAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACtt                                                                              <  1:116704/100‑1 (MQ=255)
                        tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                             <  2:3680604/100‑1 (MQ=255)
                            gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGAt                                                                         <  1:3548125/100‑1 (MQ=255)
                              ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGt                                                                       >  2:2954932/1‑100 (MQ=255)
                                 agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                    <  2:1406904/100‑1 (MQ=255)
                                  gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   <  1:2398770/100‑1 (MQ=255)
                                  gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   <  2:487189/100‑1 (MQ=255)
                                    gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                 <  1:581022/100‑1 (MQ=255)
                                    gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                 <  1:3073223/100‑1 (MQ=255)
                                        gcaTTTAATGAGGGGAGTGGTGGACGCGAGCATGGTCATTTGACCAGGCGCAGTGTTTAAAGAGTGGGATGTAACCCTATCTCAAggagatga                                                                    <  2:3115204/91‑8 (MQ=255)
                                        gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                             <  2:3416788/100‑1 (MQ=255)
                                         tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                            <  2:3225904/100‑1 (MQ=255)
                                          aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                           <  2:1963479/100‑1 (MQ=255)
                                          aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                           <  1:2135659/100‑1 (MQ=255)
                                              cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGCTTCGGCgcgagcg                                                       <  2:2360991/100‑1 (MQ=255)
                                               aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                      >  1:1898339/1‑100 (MQ=255)
                                               aTGAAGTGATGGGTGGAAACGATCATGGTCATATGACCAGGCGCAGAGCTTTACACCTGGCATGTAACCCTAACTCAAGATGTTGATTCGgcgcgcgc                                                        >  2:2905893/1‑98 (MQ=255)
                                                tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                     >  2:3179484/1‑100 (MQ=255)
                                                tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                     <  1:1419334/100‑1 (MQ=255)
                                                 gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                    >  1:895046/1‑100 (MQ=255)
                                                    gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                  >  2:2629248/1‑99 (MQ=255)
                                                      gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                               >  1:3279635/1‑100 (MQ=255)
                                                       aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                              <  2:214967/100‑1 (MQ=255)
                                                           gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                          <  2:1209471/100‑1 (MQ=255)
                                                            tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                         <  1:2898383/100‑1 (MQ=255)
                                                             gaaaaCGATCATTATCATTTGACCAGGCGCAAAGCTTTACGCGTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        >  1:2578257/3‑100 (MQ=255)
                                                              gAAACGATCATTGTCATGTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                       <  1:500873/100‑1 (MQ=255)
                                                                  cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGCTGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                                    >  1:2410440/1‑99 (MQ=255)
                                                                  cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                   <  1:253892/100‑1 (MQ=255)
                                                                   gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGTCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                  <  2:16370/100‑1 (MQ=255)
                                                                   gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                  >  1:1555666/1‑100 (MQ=255)
                                                                       aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                              <  1:400455/100‑1 (MQ=255)
                                                                        ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             >  2:1391266/1‑100 (MQ=255)
                                                                        ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             <  1:3666252/100‑1 (MQ=255)
                                                                        ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             <  2:1220686/100‑1 (MQ=255)
                                                                        ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             >  1:255865/1‑100 (MQ=255)
                                                                          gtgaTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                           <  2:3073490/97‑1 (MQ=255)
                                                                          ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                           >  2:603615/9‑100 (MQ=25)
                                                                            ttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                         >  2:916959/7‑100 (MQ=37)
                                                                            cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                         <  2:2343611/100‑1 (MQ=255)
                                                                             aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                        >  1:3704659/1‑100 (MQ=255)
                                                                               ttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                      >  2:2249952/1‑100 (MQ=255)
                                                                               cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAACATATCGTAAAAAGACAGTTCAcc                       >  2:1435905/4‑99 (MQ=37)
                                                                               cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                      >  1:2698897/4‑100 (MQ=38)
                                                                               cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                      >  1:1553858/4‑100 (MQ=38)
                                                                               cataCCAGGCGTAAAGCTCTGCGCCTGGTATGGAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                      >  1:2523635/4‑100 (MQ=38)
                                                                                tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  2:1921989/1‑100 (MQ=255)
                                                                                tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:2045791/1‑100 (MQ=255)
                                                                                tGACCAGGAGCAGAGCTTTACGACTGGTATGTAAACCTAAGTTAAGATGTTGATTCGGAGAGAGCGAGGGTAAAAATATCGTAAAAAGAGAGATCa                         >  2:2924225/1‑96 (MQ=255)
                                                                                ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:321819/3‑100 (MQ=38)
                                                                                  aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAGCTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  2:2042480/1‑100 (MQ=38)
                                                                                  aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  2:881253/1‑100 (MQ=38)
                                                                                  aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACATTTCACCttt                    >  2:809985/1‑99 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCGTGGTATGTAACCCTAACTTAAGCTGTTGATTGGGGGAGAGCGTGGATAAAAATATCGTAAACAGATAGTTAATCtt                     >  2:3548051/1‑97 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  1:1910323/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  1:162292/1‑100 (MQ=255)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   <  2:2138375/99‑1 (MQ=255)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  1:3505206/1‑100 (MQ=255)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  <  2:3748505/100‑1 (MQ=255)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:3572642/1‑100 (MQ=255)
                                                                                   ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  1:2968100/1‑100 (MQ=255)
                                                                                    cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                 >  2:1035524/1‑100 (MQ=255)
                                                                                     aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATGGTAAAAAGACAGTTCACCttt                    >  2:652076/1‑96 (MQ=255)
                                                                                     aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg                >  1:1254801/1‑100 (MQ=255)
                                                                                      ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg                >  2:2436434/1‑99 (MQ=255)
                                                                                      ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa               >  1:3394829/1‑100 (MQ=255)
                                                                                       gcgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc              <  2:2019091/100‑1 (MQ=255)
                                                                                          cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa            >  1:1757505/1‑99 (MQ=255)
                                                                                          cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt           >  2:1550197/1‑100 (MQ=255)
                                                                                            gagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa         <  1:2018115/100‑1 (MQ=255)
                                                                                              gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCCCCTTTACGGACGATCAtt       >  1:317243/1‑100 (MQ=255)
                                                                                              gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt       >  2:2592691/1‑100 (MQ=255)
                                                                                              gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt       >  1:216252/1‑100 (MQ=255)
                                                                                               cTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttt      >  1:2887708/1‑100 (MQ=255)
                                                                                               cTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttt      <  2:898751/100‑1 (MQ=255)
                                                                                                tttACGCCGGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAGCAt        <  2:1031476/97‑1 (MQ=255)
                                                                                                 ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt    <  2:253052/100‑1 (MQ=255)
                                                                                                  tACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt   >  2:1951124/1‑100 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  >  2:535348/1‑100 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  2:1919103/100‑1 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  >  1:1440865/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478660‑478858

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: