Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 874,628:1 +G 5.5% coding (1650/2349 nt) yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913874,6281.G 5.5% 494.2 / 22.4 109coding (1650/2349 nt)yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base . (52/51);  new base G (3/3);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 112 total aligned reads displayed.

CGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAA  >  NC_000913/874480‑874776
                                                                                                                                                     |                                                                                                                                                    
cGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATc                                                                                                                                                       <  2:1114790/149‑1 (MQ=255)
   tttCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGTGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGtt                                                                                                                                                   <  1:1223219/149‑1 (MQ=255)
        aGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGatat                                                                                                                                              <  1:970327/149‑1 (MQ=255)
         gTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATAtt                                                                                                                                             >  1:945190/1‑149 (MQ=255)
         gTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACGATCGGTTGat                                                                                                                                                <  3:450218/146‑1 (MQ=255)
         gTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACGATCGGTTGat                                                                                                                                                >  4:450218/1‑146 (MQ=255)
            tCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAac                                                                                                                                          <  2:395716/149‑1 (MQ=255)
               gTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGGAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGatat                                                                                                                                              <  4:1492670/142‑1 (MQ=255)
               gTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGGAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGatat                                                                                                                                              >  3:1492670/1‑142 (MQ=255)
               gTGACCATACGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACAcc                                                                                                                                       >  4:1800409/1‑149 (MQ=255)
                     atatGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGcc                                                                                                                                 >  3:1359290/1‑149 (MQ=255)
                      tatGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCa                                                                                                                                >  3:1824733/1‑149 (MQ=255)
                          ttCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCt                                                                                                                            >  4:896155/1‑149 (MQ=255)
                            ctcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTgg                                                                                                                          <  2:945190/149‑1 (MQ=255)
                             tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTggg                                                                                                                         <  1:1081618/149‑1 (MQ=255)
                             tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTggg                                                                                                                         >  2:1203823/1‑149 (MQ=255)
                             tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTggg                                                                                                                         >  4:787269/1‑149 (MQ=255)
                               agagTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGgcg                                                                                                                       >  2:1354448/1‑149 (MQ=255)
                                   tttATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGa                                                                                                                   >  3:724663/1‑149 (MQ=255)
                                    ttATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAg                                                                                                                  <  3:787269/149‑1 (MQ=255)
                                         ttCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGtt                                                                                                             <  3:896155/149‑1 (MQ=255)
                                             gtgATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGc                                                                                                         <  1:1203823/149‑1 (MQ=255)
                                                 tgaCGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGtt                                                                                                     >  2:62050/1‑149 (MQ=255)
                                                 tgaCGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGtt                                                                                                     <  1:1354448/149‑1 (MQ=255)
                                                   aCGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTGATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTgg                                                                                                   <  4:1244752/149‑1 (MQ=255)
                                                    cgcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTggg                                                                                                  <  1:281329/149‑1 (MQ=255)
                                                     gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTaa                                                                                                                                           <  1:381801/107‑1 (MQ=255)
                                                     gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTaa                                                                                                                                           >  2:381801/1‑107 (MQ=255)
                                                     gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCgcgc                                                                                                                                  <  2:1542799/116‑1 (MQ=255)
                                                     gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCgcgc                                                                                                                                  >  1:1542799/1‑116 (MQ=255)
                                                       gAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAac                                                                                                                                          <  4:809836/106‑1 (MQ=255)
                                                        aaaCTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTc                                                                                              <  2:1634039/149‑1 (MQ=255)
                                                            tCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCg                                                                                          >  1:1804288/1‑149 (MQ=255)
                                                              aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTaa                                                                                                                                           >  3:1118470/1‑98 (MQ=255)
                                                              aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTaa                                                                                                                                           <  4:1118470/98‑1 (MQ=255)
                                                              aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAacac                                                                                                                                        >  4:152205/1‑101 (MQ=255)
                                                              aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAacac                                                                                                                                        <  3:152205/101‑1 (MQ=255)
                                                              aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGa                                                                                        >  4:1232986/1‑149 (MQ=255)
                                                              aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGa                                                                                        >  2:351216/1‑149 (MQ=255)
                                                                 gAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCg                                                                                     >  1:1111653/1‑149 (MQ=255)
                                                                       gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCa                                                                                                                                >  2:343873/1‑100 (MQ=255)
                                                                       gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCa                                                                                                                                <  1:343873/100‑1 (MQ=255)
                                                                       gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTgg                                                                                                   >  4:646567/1‑129 (MQ=255)
                                                                       gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTgg                                                                                                   <  3:646567/129‑1 (MQ=255)
                                                                          gagaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCAt                                                                            <  4:1824733/149‑1 (MQ=255)
                                                                            gaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGata                                                                                                                                               <  3:1519725/80‑1 (MQ=255)
                                                                            gaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGata                                                                                                                                               >  4:1519725/1‑80 (MQ=255)
                                                                            gaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCAttt                                                                          >  2:1026263/1‑149 (MQ=255)
                                                                             aaTATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTc                                                                         >  3:1088863/1‑149 (MQ=255)
                                                                             aaTATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTc                                                                         >  4:247049/1‑149 (MQ=255)
                                                                                atCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGc                                                                                                         >  2:983577/1‑114 (MQ=255)
                                                                                atCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGc                                                                                                         <  1:983577/114‑1 (MQ=255)
                                                                                 tCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTg                                                                                                    >  2:912682/1‑118 (MQ=255)
                                                                                  cTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGc                                                                    <  3:247049/149‑1 (MQ=255)
                                                                                      aTTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGaa                                                                >  1:538494/1‑149 (MQ=255)
                                                                                       ttGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAg                                                               >  4:844467/1‑149 (MQ=255)
                                                                                            cgTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCAt                                                          <  2:1804288/149‑1 (MQ=255)
                                                                                             gTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTtgt                                                                                                           >  1:1353844/1‑99 (MQ=255)
                                                                                             gTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTtgt                                                                                                           <  2:1353844/99‑1 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTg                                                                                                                                                  >  2:1744981/1‑58 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTg                                                                                                                                                  <  1:1744981/58‑1 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTaa                                                                                                                                           <  3:1397981/65‑1 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTaa                                                                                                                                           >  4:1397981/1‑65 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgggg                                                                                  <  3:1524483/119‑1 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgggg                                                                                  >  4:1524483/1‑119 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                                    >  1:1830911/1‑117 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                                    <  2:1830911/117‑1 (MQ=255)
                                                                                                   ttCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCAtt                                                   >  1:401699/1‑149 (MQ=255)
                                                                                                     cGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCAttac                                                 <  2:538494/149‑3 (MQ=255)
                                                                                                        cATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                              >  1:177670/1‑149 (MQ=255)
                                                                                                         aTGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                             <  4:1088863/149‑1 (MQ=255)
                                                                                                         aTGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                             <  3:1232986/149‑1 (MQ=255)
                                                                                                            cAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                                   >  2:375452/1‑108 (MQ=255)
                                                                                                            cAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                                   <  1:375452/108‑1 (MQ=255)
                                                                                                              aTGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGatat                                        <  2:1111653/149‑1 (MQ=255)
                                                                                                                ggAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCg                                      <  1:62050/149‑1 (MQ=255)
                                                                                                                     ttGGCGAACTATTCCTCGTATATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGttt                                 >  3:1117287/1‑149 (MQ=255)
                                                                                                                        gCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTa                                >  1:423133/1‑147 (MQ=255)
                                                                                                                        gCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTa                                <  2:423133/147‑1 (MQ=255)
                                                                                                                          gAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGatat                                        <  2:318154/137‑1 (MQ=255)
                                                                                                                          gAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGatat                                        >  1:318154/1‑137 (MQ=255)
                                                                                                                           aaCTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAAt                           <  3:844467/149‑1 (MQ=255)
                                                                                                                              tATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                                   <  1:1457546/90‑1 (MQ=255)
                                                                                                                              tATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                                   >  2:1457546/1‑90 (MQ=255)
                                                                                                                              tATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCAcc                                                                      >  3:1450108/1‑103 (MQ=255)
                                                                                                                                  ccTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTgg                    <  1:1382604/149‑1 (MQ=255)
                                                                                                                                   cTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGa                         >  4:799898/1‑143 (MQ=255)
                                                                                                                                   cTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGa                         <  3:799898/143‑1 (MQ=255)
                                                                                                                                       tCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCgg                                     >  1:581847/1‑127 (MQ=255)
                                                                                                                                        cTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTcg                                                                        >  3:1790534/1‑90 (MQ=255)
                                                                                                                                        cTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTATGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCa                                                                        <  4:1790534/91‑1 (MQ=255)
                                                                                                                                        atATCAACGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                             >  3:258904/2‑118 (MQ=255)
                                                                                                                                         tATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCa             <  2:177670/149‑1 (MQ=255)
                                                                                                                                         tATCAACGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                              >  1:188547/1‑116 (MQ=255)
                                                                                                                                         tATCAACGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                              <  2:188547/116‑1 (MQ=255)
                                                                                                                                              aCCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGat                                          <  1:584746/115‑1 (MQ=255)
                                                                                                                                              aCCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGat                                          >  2:584746/1‑115 (MQ=255)
                                                                                                                                               ccGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCaa                            <  3:971650/128‑1 (MQ=255)
                                                                                                                                                    cGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTaaaa  <  4:1117287/149‑1 (MQ=255)
                                                                                                                                                    cGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTaaaa  <  4:724663/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATCGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAA  >  NC_000913/874480‑874776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: