Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,181,824 A→G 75.0% A243A (GCA→GCG yegV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,8240AG75.0% 224.3 / 75.8 124A243A (GCA→GCGyegVputative kinase
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base G (46/47);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 127 total aligned reads displayed.

CGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTA  >  NC_000913/2181725‑2181920
                                                                                                   |                                                                                                
cgAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCa                                                                                                  >  2:309040/1‑100 (MQ=255)
  aaCGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTTGTCTCGATAAAGAAGGCGCGTg                                                                                                <  2:2421656/100‑1 (MQ=255)
   aCGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTgg                                                                                               <  2:3569468/100‑1 (MQ=255)
    cGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTTCCGCGCCGTTCATCGTTCGTCTCGATAACGAAGGCGCGTgtt                                                                                              <  1:2814962/100‑3 (MQ=255)
    cGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTgg                                                                                               <  2:2789595/99‑1 (MQ=255)
     gTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTa                                                                                             <  1:3971450/100‑1 (MQ=255)
     gTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTa                                                                                             <  2:2447368/100‑1 (MQ=255)
      ttttGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTAt                                                                                            <  2:20504/100‑1 (MQ=255)
       tttGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTAtt                                                                                           <  1:24033/100‑1 (MQ=255)
         tGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTc                                                                                         <  2:499644/100‑1 (MQ=255)
         tGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTc                                                                                         >  1:943393/1‑100 (MQ=255)
          gCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCa                                                                                        >  1:2864661/1‑100 (MQ=255)
          gCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCa                                                                                        <  1:83913/100‑1 (MQ=255)
           cTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAg                                                                                       <  1:1742046/100‑1 (MQ=255)
           cTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAg                                                                                       >  2:2118067/1‑100 (MQ=255)
           cTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAg                                                                                       <  1:3980195/100‑1 (MQ=255)
           cTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAg                                                                                       >  2:551684/1‑100 (MQ=255)
            tttATCCGCAGAGATAACAACACTTGGCAAGCAATGTCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGc                                                                                      >  2:3545347/1‑100 (MQ=255)
             ttATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCa                                                                                     <  2:3769782/100‑1 (MQ=255)
                tCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAacg                                                                                  >  1:2571773/1‑100 (MQ=255)
                 ccGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGCAATTTGCCGCGCCGTTGATCGTTCGACTCGCTAAAGAAGGCGCGTGGTATTTCCGCAacga                                                                                 >  2:1011190/1‑100 (MQ=255)
                 ccGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAacga                                                                                 >  2:57274/1‑100 (MQ=255)
                 ccGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAacga                                                                                 >  1:1766060/1‑100 (MQ=255)
                  cGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAacgac                                                                                >  1:492936/1‑100 (MQ=255)
                   gCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAacgacg                                                                               >  1:2149236/1‑100 (MQ=255)
                   gCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAacgacg                                                                               >  1:1919612/1‑100 (MQ=255)
                     agagATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCATTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCt                                                                             >  1:2443512/1‑100 (MQ=255)
                       agaTAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTc                                                                           <  2:943389/100‑1 (MQ=255)
                       agaTAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTc                                                                           >  1:2189849/1‑100 (MQ=255)
                         aTAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTg                                                                         <  2:1938314/100‑1 (MQ=255)
                          tAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTgg                                                                        >  1:248598/1‑100 (MQ=255)
                          tAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTgg                                                                        >  1:3203360/1‑100 (MQ=255)
                          tAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTgg                                                                        >  1:1614583/1‑100 (MQ=255)
                          tAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCCCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTgg                                                                        >  2:2300412/1‑100 (MQ=255)
                           aacaacACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGc                                                                       <  1:2765668/100‑1 (MQ=255)
                             caacaCTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTg                                                                     <  1:2898510/100‑1 (MQ=255)
                             caacaCTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTg                                                                     <  2:3731908/100‑1 (MQ=255)
                             caacaCTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTg                                                                     >  1:1758348/1‑100 (MQ=255)
                             caacaCTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTg                                                                     <  2:3778108/100‑1 (MQ=255)
                                cacTTTGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCa                                                                   <  2:3092064/99‑1 (MQ=255)
                                cacTTGGCAAGGCATGGCTGACGAAATTTGACGCGCCGTTGATCGTTCGTCTCCATGAAGTGGGCACGTGGAATATCAGCAACGTGGTTTCTCGCTGc                                                                    <  2:510021/98‑1 (MQ=255)
                                cacTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCAt                                                                  >  1:10626/1‑100 (MQ=255)
                                  cTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTc                                                                >  2:3451285/1‑100 (MQ=255)
                                  cTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTc                                                                >  1:2141428/1‑100 (MQ=255)
                                    tGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCg                                                              <  1:2496440/100‑1 (MQ=255)
                                     ggCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCgg                                                             <  1:2701474/100‑1 (MQ=255)
                                     ggCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCgg                                                             >  2:3654750/1‑100 (MQ=255)
                                     ggCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCgg                                                             <  1:3810589/100‑1 (MQ=255)
                                         agcaaTGGCAGGAGAATTTTGCCGCGCCATTGATCTGTCGTTTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGTCGCTTCTGGCTGCATTCCgg                                                             <  1:3905659/96‑1 (MQ=255)
                                         agcaaTGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCAtt                                                         >  2:762391/1‑100 (MQ=255)
                                           caaTGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTc                                                       <  1:719620/100‑1 (MQ=255)
                                             aTGGCAGGAGAAATATGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTc                                                       <  2:2438320/98‑1 (MQ=255)
                                              tGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCaa                                                    <  1:273489/100‑1 (MQ=255)
                                               ggCAGGAGCAATTTGCCGCGCCGTTGATCGGGCGTCTCGATAAAGAAGGCGCGTGGGATTTAAGCAACGCCGCTTCTGGCTGCCTTCCGGCATTTCGAAc                                                   >  2:632808/1‑100 (MQ=255)
                                               ggCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAAc                                                   <  1:1325270/100‑1 (MQ=255)
                                                gCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACg                                                  >  1:4129546/1‑100 (MQ=255)
                                                 cAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGc                                                 >  2:4084841/1‑100 (MQ=255)
                                                   ggAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCa                                                >  2:1967780/1‑99 (MQ=255)
                                                    gagaAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAg                                              <  1:2958905/100‑1 (MQ=255)
                                                      gaAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGtt                                            >  2:3452867/1‑100 (MQ=255)
                                                       aaaTTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGtt                                            <  2:2903802/99‑1 (MQ=255)
                                                        aaTTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGt                                          <  2:711332/100‑1 (MQ=255)
                                                         aTTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTc                                         <  2:2769361/100‑1 (MQ=255)
                                                            tGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGac                                      <  1:2325028/100‑1 (MQ=255)
                                                            tGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGac                                      <  1:3021025/100‑1 (MQ=255)
                                                             gCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGac                                      <  2:1649989/99‑1 (MQ=255)
                                                              ccGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGacac                                    >  2:2858929/1‑100 (MQ=255)
                                                              ccGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGacac                                    >  1:3034664/1‑100 (MQ=255)
                                                               cgcgCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACAcc                                   >  1:11264/1‑100 (MQ=255)
                                                               cgcgCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACAcc                                   >  1:895368/1‑100 (MQ=255)
                                                                gcgcCGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCa                                  >  2:4082486/1‑100 (MQ=255)
                                                                 cgcCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCAt                                 <  1:1842305/100‑1 (MQ=255)
                                                                 cgcCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCAt                                 <  1:1842327/100‑1 (MQ=255)
                                                                    cGTTGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTgg                              >  1:1462615/1‑100 (MQ=255)
                                                                      ttGATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTgggg                            <  1:160636/100‑1 (MQ=255)
                                                                        gATCGTTCGTCTCGTTAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCAGTTCCAACGGAAGTTGTCGACACCATTGgggcg                          <  2:778647/100‑1 (MQ=255)
                                                                        gATCGTTCGTCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcg                          >  1:3922525/1‑100 (MQ=255)
                                                                          tcgttcgtCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcggg                        >  1:1238566/1‑100 (MQ=255)
                                                                          tcgttcgtCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcggg                        >  2:307205/1‑100 (MQ=255)
                                                                           cgttcgtCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcgggc                       <  2:1104524/100‑1 (MQ=255)
                                                                           cgttcgtCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcgggc                       <  1:933000/100‑1 (MQ=255)
                                                                           cgttcgtCTCGATAAAGAAGGCGCATGGTACTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcgggc                       <  1:846490/100‑1 (MQ=255)
                                                                            gttcgtCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGgggcgggc                       <  2:1268305/99‑1 (MQ=255)
                                                                            gttcgtCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGggcgggcg                      <  1:2474985/100‑1 (MQ=255)
                                                                            gttcgtCTCGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGggcgggcg                      >  1:78851/1‑100 (MQ=255)
                                                                               cgtCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACa                   <  1:1676124/100‑1 (MQ=255)
                                                                                 tctcGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGt                 <  1:3662674/100‑1 (MQ=255)
                                                                                  ctcGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTc                <  1:2231729/100‑1 (MQ=255)
                                                                                   tcGATAAAGAAGGCGCGTGGTATTTCAGCAACGGCGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCa               <  2:1864480/100‑1 (MQ=255)
                                                                                   tcGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCCGGCGACAGTCa               >  2:1773830/1‑100 (MQ=255)
                                                                                   tcGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCa               <  1:3566366/100‑1 (MQ=255)
                                                                                   tcGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCa               >  1:1245486/1‑100 (MQ=255)
                                                                                    cGATAAAGAAGGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCAt              >  1:3023371/1‑100 (MQ=255)
                                                                                         aagaagGCGCGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCgg         >  2:852348/1‑100 (MQ=255)
                                                                                         aagaagGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCgg         <  1:702013/100‑1 (MQ=255)
                                                                                            aagGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTgg      <  1:496574/100‑1 (MQ=255)
                                                                                               gcgcGTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGt   <  2:912029/100‑1 (MQ=255)
                                                                                                cgcgTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTa  >  1:433169/1‑100 (MQ=255)
                                                                                                cgcgTGGTATTTCAGCAACGACGCGTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTa  >  2:3911314/1‑100 (MQ=255)
                                                                                                cgcATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTa  >  1:2448961/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTA  >  NC_000913/2181725‑2181920

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: