Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,759 A→G 18.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,7590AG18.4% 422.6 / 83.8 194intergenic (‑142/+22)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (86/72);  new base G (21/15);  total (107/87)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTG  >  NC_000913/478660‑478858
                                                                                                   |                                                                                                   
 tgtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGTCGCAGAGCTTTa                                                                                                     >  1:2294313/1‑99 (MQ=255)
 tgtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTAc                                                                                                    <  2:1354075/100‑1 (MQ=255)
  gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                   >  1:2540658/1‑100 (MQ=255)
   tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                                  >  2:3719254/1‑100 (MQ=255)
   tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                                  >  1:192145/1‑100 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                 >  1:2750786/1‑100 (MQ=255)
     gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                >  2:1015779/1‑100 (MQ=255)
      naaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                               >  2:3326239/2‑100 (MQ=255)
      aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTg                                                                                               <  2:3010376/100‑1 (MQ=255)
      aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                               <  2:2460/100‑1 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgt                                                                                              <  2:870210/100‑2 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                              >  2:640920/1‑100 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGACCGGGTATTTCATGAAGTGATGGGTGGAAGCGATCATTGTCATTTGACCATGCGTAAAGCTCTGCGCCtag                                                                                              <  2:2364517/100‑3 (MQ=255)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                             <  1:2973928/100‑1 (MQ=255)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  2:2990367/100‑2 (MQ=38)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  1:1438585/100‑2 (MQ=38)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  1:2509389/100‑2 (MQ=38)
         aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                            <  2:113051/100‑1 (MQ=255)
            gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGt                                                                                         >  2:1209736/1‑100 (MQ=255)
            gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGCAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                            <  2:3950107/97‑2 (MQ=38)
             ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                        >  1:397107/1‑100 (MQ=255)
             ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                        <  1:2243428/100‑1 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                     >  2:3925994/1‑100 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                     <  1:2960264/100‑1 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                     >  1:237419/1‑100 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                     >  2:1300672/1‑100 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                    >  1:1462355/1‑100 (MQ=255)
                   ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                  >  2:2645067/1‑100 (MQ=255)
                     agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                                >  2:3845216/1‑100 (MQ=255)
                      gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACt                                                                               <  2:2243147/100‑1 (MQ=255)
                        tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                             >  1:3759908/1‑100 (MQ=255)
                          aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACTAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAg                                                                           >  2:1432985/1‑100 (MQ=255)
                          aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAg                                                                           >  2:1618651/1‑100 (MQ=255)
                          aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCGTTACGCCTGGTATGTAACGCTAACTTAAg                                                                           >  2:2587651/1‑100 (MQ=255)
                           ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                          <  2:159288/100‑1 (MQ=255)
                           ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                          >  2:2623634/1‑100 (MQ=255)
                             aggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATg                                                                        >  1:2872590/1‑100 (MQ=255)
                              ggTAGAGCTGGTATTTCATGAAGTGATTGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGt                                                                       >  2:790185/1‑100 (MQ=255)
                               gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                      <  1:1084653/100‑1 (MQ=255)
                                tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                     <  2:2689520/100‑1 (MQ=255)
                                  gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   <  1:1578536/100‑1 (MQ=255)
                                  gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   <  2:1462346/100‑1 (MQ=255)
                                  gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                   >  1:3116835/1‑100 (MQ=255)
                                   agCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                  >  2:2209534/1‑100 (MQ=255)
                                    gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                 >  1:1141610/1‑100 (MQ=255)
                                      tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                               >  1:2045388/1‑100 (MQ=255)
                                       ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                              <  2:225302/100‑1 (MQ=255)
                                        gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcg                                                             <  2:3723955/100‑1 (MQ=255)
                                          aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                           <  1:3450868/100‑1 (MQ=255)
                                             tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                        <  2:854456/100‑1 (MQ=255)
                                              cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                       >  2:2973899/1‑100 (MQ=255)
                                               aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                      >  1:3226156/1‑100 (MQ=255)
                                               aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                      >  2:3796010/1‑100 (MQ=255)
                                                 gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                    >  2:2333120/1‑100 (MQ=255)
                                                  aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                   <  1:2577025/100‑1 (MQ=255)
                                                    gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                 <  1:2886240/100‑1 (MQ=255)
                                                        tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                             >  1:1216805/1‑100 (MQ=255)
                                                         gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                            >  1:1248209/1‑100 (MQ=255)
                                                          ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                           >  1:2923201/1‑100 (MQ=255)
                                                          ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                           <  2:2798529/100‑1 (MQ=255)
                                                          ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                           <  1:582128/100‑1 (MQ=255)
                                                           gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                          >  2:1949323/1‑100 (MQ=255)
                                                           gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                          >  2:3225197/1‑100 (MQ=255)
                                                            tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                         <  2:3361884/100‑1 (MQ=255)
                                                             ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                        <  1:3325283/100‑1 (MQ=255)
                                                               ttccGATCTTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                      <  2:3192862/97‑1 (MQ=255)
                                                                       aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                              >  1:2501817/1‑100 (MQ=255)
                                                                        ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                             <  2:3870823/100‑1 (MQ=255)
                                                                            cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                         <  2:55179/100‑1 (MQ=255)
                                                                             aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                        >  2:1524246/1‑100 (MQ=255)
                                                                               ttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAATTCACCt                      >  1:2926763/1‑100 (MQ=255)
                                                                               cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                      >  1:3066542/4‑100 (MQ=38)
                                                                                ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTATCCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  2:2607960/3‑100 (MQ=37)
                                                                                ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  1:3785089/3‑100 (MQ=38)
                                                                                ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                     >  2:2154340/3‑100 (MQ=38)
                                                                                 taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAGCTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                    >  2:1470165/2‑100 (MQ=38)
                                                                                 gACCGGGCGGATGGCTTTGCGCCTGGTATTTAACTCTTACTTATTATGTTGATTCGGCGCGGGCGTGGGTAAATGTATCGTAAAAAGACAGTTCACCttt                    >  2:3144049/1‑100 (MQ=255)
                                                                                 gACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                    >  2:433392/1‑100 (MQ=255)
                                                                                 gACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                    >  1:1759778/1‑100 (MQ=255)
                                                                                  aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGTCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  1:1009739/1‑100 (MQ=38)
                                                                                  aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                   >  1:99139/1‑100 (MQ=38)
                                                                                  aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACACTACCTTAAGATGTTGCTTCGGCGCAAGCGAGGGTACAAATATCGTACAAAGACAGATCAcccaca                   >  2:2710506/1‑96 (MQ=19)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  1:3626592/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:3475484/1‑100 (MQ=255)
                                                                                   ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                  >  2:240043/1‑100 (MQ=255)
                                                                                    cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                 >  1:2611791/1‑100 (MQ=255)
                                                                                       gcgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc              <  2:1088419/100‑1 (MQ=255)
                                                                                       gcgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc              >  1:521636/1‑100 (MQ=255)
                                                                                       gcgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc              >  1:456482/1‑100 (MQ=255)
                                                                                          cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt           >  1:3703142/1‑100 (MQ=255)
                                                                                          cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt           >  1:2804795/1‑100 (MQ=255)
                                                                                              gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt       >  1:2810138/1‑100 (MQ=255)
                                                                                              gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCGTTACGGACGATCAtt       >  1:3887279/1‑100 (MQ=255)
                                                                                                tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt     <  1:777551/100‑1 (MQ=255)
                                                                                                  tACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAGAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt   >  2:2551013/1‑100 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  1:1818837/100‑1 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  1:3172430/100‑1 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  >  2:708619/1‑100 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  >  1:3303623/1‑100 (MQ=255)
                                                                                                   aCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg  <  1:640923/100‑1 (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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: