Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,757 T→C 15.7% 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.7% 379.9 / 61.8 169intergenic (‑140/+24)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (65/75);  new base C (18/8);  total (86/83)
Fisher's exact test for biased strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 176 total aligned reads displayed.

GATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478658‑478856
                                                                                                   |                                                                                                   
gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCtt                                                                                                     >  1:1475951/1‑100 (MQ=255)
gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCtt                                                                                                     <  2:2787519/100‑1 (MQ=255)
  ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                                   >  1:960491/1‑100 (MQ=255)
      ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAATGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                               <  2:1942709/100‑1 (MQ=255)
      ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                               <  2:915925/100‑1 (MQ=255)
      ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                               >  1:197782/1‑100 (MQ=255)
       gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                              >  2:18805/1‑100 (MQ=255)
        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             <  2:965241/100‑1 (MQ=255)
        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             <  1:3479700/100‑1 (MQ=255)
        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             <  2:2405751/100‑1 (MQ=255)
         aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                            >  2:360416/1‑100 (MQ=255)
         aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                            <  1:2760696/100‑1 (MQ=255)
          aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                           <  1:3334364/100‑1 (MQ=38)
           aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                          <  2:1451889/100‑1 (MQ=255)
           aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                          >  1:3537840/1‑100 (MQ=255)
             ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                        >  2:2314564/1‑100 (MQ=255)
                 ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                    <  1:669065/100‑1 (MQ=255)
                 ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                    >  2:3308877/1‑100 (MQ=255)
                  gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaatgac                                                                                   <  2:2611313/100‑9 (MQ=25)
                  gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTCAcc                                                                                   <  2:1764988/100‑1 (MQ=255)
                  gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                   >  2:2346180/1‑100 (MQ=255)
                  gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCACAGAGCTATACGCCTGGTATGTAAcc                                                                                   >  1:790787/1‑100 (MQ=255)
                    tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCt                                                                                 >  1:1828103/1‑100 (MQ=255)
                    tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCt                                                                                 <  2:1475943/100‑1 (MQ=255)
                       agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                              >  1:1582648/1‑100 (MQ=255)
                        gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACt                                                                             >  1:2422206/1‑100 (MQ=255)
                          tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                           <  1:1757514/100‑1 (MQ=255)
                            aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAg                                                                         <  1:736451/100‑1 (MQ=255)
                             ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAACCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                        <  2:174954/100‑1 (MQ=255)
                             ggaggTAGAGCTGGTATTGCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGAGCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                        >  2:2174247/1‑100 (MQ=255)
                                  ttgagCTTGTAATTCATGAAGTGACGGATCGAAACGATCATTGTCAGGTGGCCAGGAGTAGACCTATACGCCTGGTAGGTAACGCTTACTGTAGATGTTg                                                                   >  2:1498675/3‑100 (MQ=255)
                                  tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                   <  1:2721901/100‑1 (MQ=255)
                                        tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                             <  1:185762/100‑1 (MQ=255)
                                        tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                             >  2:3278140/1‑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)
                                           tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTCACATGTTGATTCGgcgc                                                          >  2:2890560/1‑100 (MQ=255)
                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         >  2:3053975/1‑100 (MQ=255)
                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         <  1:2229409/100‑1 (MQ=255)
                                             tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcga                                                        >  1:2102563/1‑100 (MQ=255)
                                                cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                     >  1:2145758/1‑100 (MQ=255)
                                                   gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                  <  1:2065677/100‑1 (MQ=255)
                                                   gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                  <  1:3044222/100‑1 (MQ=255)
                                                    aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                 <  1:3097774/100‑1 (MQ=255)
                                                         aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                            >  2:1371199/1‑100 (MQ=255)
                                                         aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                            >  2:1205705/1‑100 (MQ=255)
                                                           gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                          <  1:2399054/100‑1 (MQ=255)
                                                             gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTCACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        >  2:1603271/1‑100 (MQ=255)
                                                             gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        <  1:2279442/100‑1 (MQ=255)
                                                              tagAAACGCTCATTGTCATTTGACCAGGAGCAGAGCTTTACGCCTGGTATGTAACCCTCACTTAAGATGTTGATTCGGCGCGAGAGAGGGTAAAAATATc                                       <  2:2148608/98‑1 (MQ=255)
                                                              tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                       <  2:1106194/100‑1 (MQ=255)
                                                              tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                       <  1:1334705/100‑1 (MQ=255)
                                                              tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                       <  2:1616343/100‑1 (MQ=255)
                                                              tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                       >  1:1078834/1‑100 (MQ=255)
                                                                gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                     <  1:393888/100‑1 (MQ=255)
                                                                gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                     >  1:1157924/1‑100 (MQ=255)
                                                                gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                     <  1:3290890/100‑1 (MQ=255)
                                                                gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                     >  1:3260734/1‑100 (MQ=255)
                                                                 aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                    >  1:309423/1‑100 (MQ=255)
                                                                 aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                    >  1:1025712/1‑100 (MQ=255)
                                                                   aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                                  <  2:2625493/100‑1 (MQ=255)
                                                                    cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                 >  2:2741214/1‑100 (MQ=255)
                                                                     gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGGTGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                <  2:1562261/100‑1 (MQ=255)
                                                                     gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                <  1:3325887/100‑1 (MQ=255)
                                                                         aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                            <  2:124458/100‑1 (MQ=255)
                                                                         aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                            >  1:2193198/1‑100 (MQ=255)
                                                                          ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                           <  2:2422190/100‑1 (MQ=255)
                                                                          ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                           <  1:367707/100‑1 (MQ=255)
                                                                            ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                           >  2:1716853/9‑98 (MQ=25)
                                                                            gTCATTTGACCAGGCGCAGAGCTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                        <  2:3089771/100‑1 (MQ=255)
                                                                              cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                       >  1:1988287/1‑100 (MQ=255)
                                                                               tacataCCAGGCGTAAAGCTCTGCGCCTGGTATATAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAAGGTAAAAATATCGTAAAATGACAGTTccc                      >  2:2185510/6‑98 (MQ=25)
                                                                                tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                     >  1:1044480/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)
                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  2:2583227/2‑100 (MQ=38)
                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  2:3232235/2‑100 (MQ=38)
                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACATAAGATGTTGATTCGGCGCGACCGAGGGTAAAAATATCGTAACAAGACAGCTCACCttt                  >  1:882749/2‑100 (MQ=37)
                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 <  2:2411404/100‑1 (MQ=255)
                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 >  2:3583120/1‑100 (MQ=255)
                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 <  1:1379453/100‑1 (MQ=255)
                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 <  1:3290819/100‑1 (MQ=255)
                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 >  2:209301/1‑100 (MQ=255)
                                                                                     ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  2:16499/1‑100 (MQ=255)
                                                                                     ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  1:404729/1‑100 (MQ=255)
                                                                                     ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  1:2633811/1‑100 (MQ=255)
                                                                                     ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCTAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                <  2:1988276/100‑1 (MQ=255)
                                                                                        ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa             <  1:2793190/100‑1 (MQ=255)
                                                                                         gcacAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc            >  2:18667/4‑100 (MQ=255)
                                                                                          cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           >  1:3158595/1‑100 (MQ=255)
                                                                                              gagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa       <  1:2494050/100‑1 (MQ=255)
                                                                                                gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt     >  1:3476358/1‑100 (MQ=255)
                                                                                                gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt     <  2:675211/100‑1 (MQ=255)
                                                                                                 cTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttt    <  1:2314580/100‑1 (MQ=255)
                                                                                                  tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt   >  2:3540719/1‑100 (MQ=255)
                                                                                                   ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  <  2:3223427/100‑1 (MQ=255)
                                                                                                   ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  >  2:391324/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478658‑478856

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: