Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,750 C→T 18.5% 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,7500CT18.5% 318.1 / 55.2 142intergenic (‑133/+31)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (62/52);  new base T (16/10);  total (79/63)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 151 total aligned reads displayed.

ATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAT  >  NC_000913/478654‑478849
                                                                                                |                                                                                                   
atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCaga                                                                                                  >  1:760882/1‑100 (MQ=255)
    gATTGTGGAAAAAGGCTGGCTGGAGATAAGGGGGTAGGGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGGAGAGCtt                                                                                              >  1:3615330/1‑100 (MQ=255)
        gtgGAAACAGCCTGGATGGAAATAAGTAGGTAGATCTGGTATTTCATGTAGTGATGGGTGCTAACGATCATTGTCATTTTACTAGTCGCAGAGCTTTACg                                                                                          <  1:1315206/100‑1 (MQ=255)
        gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                          >  2:2090600/1‑100 (MQ=255)
        gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                          >  1:3577715/1‑100 (MQ=255)
        gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                          >  2:1556207/1‑100 (MQ=255)
        gtgGAAAAAGGCTGGATGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                          <  2:3759282/100‑1 (MQ=255)
          ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        >  1:1765776/1‑100 (MQ=255)
          ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGATGGTATTTCATGAAGTGATGGGTGGAAACGATCATTCTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                        >  1:2299506/1‑100 (MQ=255)
           gAAAAAGACTGGCTGGAGATAAGGAGGTCGAGCTGGTATTTAATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                          >  1:1309714/1‑97 (MQ=255)
            aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                      <  2:924103/100‑1 (MQ=255)
            aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                      <  2:2965986/100‑1 (MQ=255)
            aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                      >  2:1666150/1‑100 (MQ=255)
             aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                     <  2:57222/100‑1 (MQ=255)
              aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                    <  2:1361773/100‑1 (MQ=38)
              aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                    <  1:2830133/100‑1 (MQ=255)
               aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                   <  2:618487/100‑2 (MQ=38)
               aaGGCTGGCTAGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGAGCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                   <  1:2531970/100‑2 (MQ=37)
                aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                  <  2:1980857/100‑3 (MQ=38)
                aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                  <  1:2096561/100‑3 (MQ=38)
                aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                  <  2:1633150/100‑3 (MQ=38)
                aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTAt                                                                                  >  1:3727317/1‑100 (MQ=255)
                 ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                 >  1:364234/1‑100 (MQ=255)
                 ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                 <  2:441631/100‑1 (MQ=255)
                  gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGt                                                                                <  2:2549124/100‑1 (MQ=255)
                   ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                               >  1:2771282/1‑100 (MQ=255)
                     ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             <  2:1356729/100‑1 (MQ=255)
                     ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             <  2:2970398/100‑1 (MQ=255)
                      gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                             <  2:3064257/99‑1 (MQ=255)
                       cttccgatctggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                           <  1:1189073/90‑1 (MQ=255)
                       ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                           <  2:2047040/100‑1 (MQ=255)
                       ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                           >  1:176773/1‑100 (MQ=255)
                          gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                        >  1:1001193/1‑100 (MQ=255)
                          gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                        >  2:2479682/1‑100 (MQ=255)
                            gaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACt                                                                      >  1:435176/1‑100 (MQ=255)
                              tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                    <  1:891717/100‑1 (MQ=255)
                               aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                   <  2:3352817/100‑1 (MQ=255)
                                 ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcccagatcggaag                                                                 >  2:1189080/1‑90 (MQ=255)
                                 ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                 <  1:2198880/100‑1 (MQ=255)
                                 ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                 <  1:2375453/100‑1 (MQ=255)
                                    ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGt                                                              >  2:727980/1‑100 (MQ=255)
                                    ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGt                                                              >  2:2224737/1‑100 (MQ=255)
                                        gtgCTGGTGTTTCATGAAGTGATGGGTGGAATCGATCATTGTCATTTGTCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                          <  1:980304/98‑1 (MQ=255)
                                        gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                          >  1:1672955/1‑100 (MQ=255)
                                         agCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                         >  2:3100314/1‑100 (MQ=255)
                                           cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCg                                                       >  2:1105072/1‑100 (MQ=255)
                                            tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                      <  2:3495799/100‑1 (MQ=255)
                                             ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                     >  1:359890/1‑100 (MQ=255)
                                             ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                     >  1:3192929/1‑100 (MQ=255)
                                                aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                  >  1:666263/1‑100 (MQ=255)
                                                aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                  <  1:784490/100‑1 (MQ=255)
                                                       gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                           >  1:2943510/1‑100 (MQ=255)
                                                        aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                          <  1:1502932/100‑1 (MQ=255)
                                                         aGTGATGGGTGGAAACTATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                         <  2:315268/100‑1 (MQ=255)
                                                          gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                        <  2:364239/100‑1 (MQ=255)
                                                           tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaa                                       >  2:1988870/1‑100 (MQ=255)
                                                            gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                      >  2:3047631/1‑100 (MQ=255)
                                                              tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGATGGTaaaaa                                    >  1:285712/1‑100 (MQ=255)
                                                                ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTCAAGACGTTGATTCGGCGCGAGCGAGGATCAAAata                                  >  1:1947284/1‑100 (MQ=255)
                                                                 gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                 >  2:2888230/1‑100 (MQ=255)
                                                                 gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                 >  1:1916377/1‑100 (MQ=255)
                                                                  tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                >  1:1131666/1‑100 (MQ=255)
                                                                  tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                <  2:1926987/100‑1 (MQ=255)
                                                                   ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                               <  2:3616916/100‑1 (MQ=255)
                                                                   ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                               <  2:1443995/100‑1 (MQ=255)
                                                                        cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                          >  2:43564/1‑100 (MQ=255)
                                                                          aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                        >  1:1831051/1‑100 (MQ=255)
                                                                          aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                        <  1:944710/100‑1 (MQ=255)
                                                                           tCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                       <  2:273389/100‑1 (MQ=255)
                                                                              ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTGACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                    >  1:2120683/1‑100 (MQ=255)
                                                                                ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGAGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                  >  1:2197932/9‑100 (MQ=25)
                                                                                gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                  >  2:2417061/1‑100 (MQ=255)
                                                                                gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                  >  1:1273398/1‑100 (MQ=255)
                                                                                    tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc              <  2:409148/100‑1 (MQ=255)
                                                                                    acataCCAGGCGTAAAGCTCTGCGCCTGGTCTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc              >  2:1095448/5‑100 (MQ=37)
                                                                                     ttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt             <  2:1248513/100‑1 (MQ=255)
                                                                                      tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt            <  1:1744055/100‑1 (MQ=255)
                                                                                       taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  1:1843399/2‑100 (MQ=38)
                                                                                       taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  2:1734169/2‑100 (MQ=38)
                                                                                       taCCAGGCGTAAAGCTATGCGCTTGGAATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAGAAGACAGtcaatcttt           >  1:3700881/2‑92 (MQ=255)
                                                                                       gACCAGGAGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGGGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  1:2005808/1‑100 (MQ=255)
                                                                                        acgaGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTGActgtt           >  1:2664585/4‑95 (MQ=255)
                                                                                        aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt           >  1:2840669/1‑99 (MQ=38)
                                                                                        aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa          >  1:1546215/1‑100 (MQ=38)
                                                                                         ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  2:2056780/1‑100 (MQ=255)
                                                                                         ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         <  2:1928356/100‑1 (MQ=255)
                                                                                         ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         >  2:3355594/1‑100 (MQ=255)
                                                                                         ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc         <  1:3798760/100‑1 (MQ=255)
                                                                                         ccAGACGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAGTATCGTAAAACGACAGTTCACCttt           >  1:2880789/1‑98 (MQ=255)
                                                                                          cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATTTTGATTCGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGACAGTTCACCTTTACg        >  2:1254592/1‑100 (MQ=255)
                                                                                          cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg        >  2:3205055/1‑100 (MQ=255)
                                                                                          cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg        <  1:2791363/100‑1 (MQ=255)
                                                                                            ggCGTAGAGCTTTACGCCTGGTATGGAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg        <  1:1667080/98‑1 (MQ=255)
                                                                                            ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa      >  2:832959/1‑100 (MQ=255)
                                                                                               gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa   <  2:3466975/100‑1 (MQ=255)
                                                                                               gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa   >  2:3126269/1‑100 (MQ=255)
                                                                                                cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt  >  1:2968741/1‑100 (MQ=255)
                                                                                                cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt  <  1:2985643/100‑1 (MQ=255)
                                                                                                cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt  <  2:2299517/100‑1 (MQ=255)
                                                                                                cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt  <  2:1908097/100‑1 (MQ=255)
                                                                                                |                                                                                                   
ATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAT  >  NC_000913/478654‑478849

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: