Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 564,046 T→G S145S (TCA→TCC sfmZ ← response regulator family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913564,0460TG98.5% 696.0 / 0.3 195S145S (TCA→TCCsfmZresponse regulator family protein
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base G (101/91);  total (101/94)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 195 total aligned reads displayed.

TGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTTGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGAAAA  >  NC_000913/563949‑564142
                                                                                                 |                                                                                                
tGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGCTAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGa                                                                                                <  1:5245891/100‑1 (MQ=255)
tGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGa                                                                                                >  1:5397656/1‑100 (MQ=255)
 gCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGaa                                                                                               >  2:5085555/1‑100 (MQ=255)
 gCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGaa                                                                                               >  1:1525199/1‑100 (MQ=255)
  cTAAGTAATAACTTATCGGCAATTTCTTTATTAGATATTCCGCTAATCAGATAACGTAATATGGTTACTTCACGATTAGATAGCACAGTGACCGTGGAAc                                                                                              <  1:5179064/100‑1 (MQ=255)
  cTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAAc                                                                                              >  1:740779/1‑100 (MQ=255)
   tAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACt                                                                                             >  2:3288873/1‑100 (MQ=255)
   tAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACt                                                                                             >  1:3083543/1‑100 (MQ=255)
      gTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTAtt                                                                                          <  1:270544/100‑1 (MQ=255)
       taataaCTTATCGTCAATTTCTTTATTAGATAATCCGCTAACCAGACAACGTAATATGGTCACTTCACCATTAGATAGCACAGTGACCGTGGAACTATTc                                                                                         <  2:1190791/100‑1 (MQ=255)
         ataaCTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGt                                                                                       <  1:1086103/100‑1 (MQ=255)
         ataaCTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGt                                                                                       >  2:5220646/1‑100 (MQ=255)
            aCTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGtact                                                                                    <  2:2109625/100‑1 (MQ=255)
            aCTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGtact                                                                                    >  1:4048313/1‑100 (MQ=255)
             cTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGtacta                                                                                   >  2:308197/1‑100 (MQ=255)
             cTTATCGGCAATTTCTTTATGAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGtacta                                                                                   >  1:4341191/1‑100 (MQ=255)
              ttATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGtactac                                                                                  >  1:3689596/1‑100 (MQ=255)
                   ggCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTa                                                                             >  2:4046252/1‑100 (MQ=255)
                     cAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTAtt                                                                           <  2:1553229/100‑1 (MQ=255)
                        tttCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCt                                                                        >  1:2909899/1‑100 (MQ=255)
                        tttCTTTATTAGATAATCCGCTAACCAGATAACCTAATATGGTCACTTCACGATTAGACAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCt                                                                        <  2:660220/100‑1 (MQ=255)
                         ttCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCtt                                                                       <  1:3435354/100‑1 (MQ=255)
                         ttCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCtt                                                                       >  2:3319640/1‑100 (MQ=255)
                         ttCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACAGTGGAACTATTCGTACTACATTTATTGCtt                                                                       >  1:5561782/1‑100 (MQ=255)
                          tCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCttt                                                                      <  1:2340972/100‑1 (MQ=255)
                           cTTTATTAGATAATCCGCTAACCAGATAACGTGATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCtttt                                                                     >  1:381969/1‑100 (MQ=255)
                           cTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCtttt                                                                     <  1:2437345/100‑1 (MQ=255)
                           cTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCtttt                                                                     >  1:428338/1‑100 (MQ=255)
                            tttATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCttttt                                                                    <  1:1695072/100‑1 (MQ=255)
                             ttattaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTGGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTta                                                                   <  1:4959775/100‑1 (MQ=255)
                              tattaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtat                                                                  <  2:3207550/100‑1 (MQ=255)
                              tattaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtat                                                                  <  1:3559166/100‑1 (MQ=255)
                              tattaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtat                                                                  <  1:2510828/100‑1 (MQ=255)
                              tattaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtat                                                                  >  2:2696015/1‑100 (MQ=255)
                              tattaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtat                                                                  <  2:4249985/100‑1 (MQ=255)
                               attaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtata                                                                 >  2:4722783/1‑100 (MQ=255)
                               attaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtata                                                                 >  1:1315486/1‑100 (MQ=255)
                               attaGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtata                                                                 >  1:3322492/1‑100 (MQ=255)
                                 taGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTtatata                                                               >  2:3194993/1‑100 (MQ=255)
                                   gATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGt                                                             <  2:1992354/100‑1 (MQ=255)
                                     tAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTa                                                           >  2:3829737/1‑100 (MQ=255)
                                      aaTCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTaa                                                          >  1:3100825/1‑100 (MQ=255)
                                       aTCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTATATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAg                                                         <  1:64354/100‑1 (MQ=255)
                                       aTCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAg                                                         >  1:954738/1‑100 (MQ=255)
                                        tCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGc                                                        <  1:3496503/100‑1 (MQ=255)
                                         cccccaaCCAGACAACATAACAGTTTCACTGCACGAGTAGAGAGCACAGGGACCGTGGAACTATTCGGACTACATTTATTGCTTTGGATATAGTTAATCg                                                       <  2:5487875/95‑1 (MQ=255)
                                             tAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTc                                                   <  1:3865108/100‑1 (MQ=255)
                                             tAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTc                                                   <  1:5026435/100‑1 (MQ=255)
                                               aCCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTTGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGc                                                 <  2:1441912/100‑1 (MQ=255)
                                               aCCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGc                                                 <  1:1214463/100‑1 (MQ=255)
                                               aCCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGc                                                 >  2:1728227/1‑100 (MQ=255)
                                                ccAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCt                                                <  2:2167225/100‑1 (MQ=255)
                                                  aGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTCAATTTATTGCTTTTTATATATCTAACCGTCTCGCttt                                              >  2:4592373/1‑98 (MQ=255)
                                                  aGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTgg                                              <  2:3617635/100‑1 (MQ=255)
                                                   gATAACGTAATATGGTCACTTCACTATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTggg                                             <  2:219218/100‑1 (MQ=255)
                                                    aTAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGa                                            >  1:3956700/1‑100 (MQ=255)
                                                     tAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGaa                                           <  2:3438280/100‑1 (MQ=255)
                                                       aCGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGaaaa                                         >  2:4294423/1‑100 (MQ=255)
                                                       aCGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGaaaa                                         <  1:1884604/100‑1 (MQ=255)
                                                           aaTATGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACg                                     <  1:1064534/100‑1 (MQ=255)
                                                            atatGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACgt                                    <  2:4630432/100‑1 (MQ=255)
                                                             tatGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACgtg                                   <  2:894563/100‑1 (MQ=255)
                                                             tatGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACgtg                                   <  1:2740453/100‑1 (MQ=255)
                                                             tatGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACgtg                                   <  1:3044746/100‑1 (MQ=255)
                                                              atGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACgtgt                                  >  2:2400960/1‑100 (MQ=255)
                                                              atGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACgtgt                                  <  1:1727029/100‑1 (MQ=255)
                                                               tGGTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGAAAAAAACGTGTa                                 <  1:4789495/100‑1 (MQ=255)
                                                                ggTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACGTTTATTGCTTTTTGTATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTAt                                >  2:2502588/1‑100 (MQ=255)
                                                                ggTCACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTAt                                >  1:3002929/1‑100 (MQ=255)
                                                                   cACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCg                             >  2:3088694/1‑100 (MQ=255)
                                                                   cACTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCg                             >  2:317146/1‑100 (MQ=255)
                                                                    aCTTCCCGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCgg                            <  1:258209/100‑1 (MQ=255)
                                                                     cTTCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGga                           >  2:404669/1‑100 (MQ=255)
                                                                      ttCACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGgag                          <  2:4118613/100‑1 (MQ=255)
                                                                       tCACCATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTGATATAGTTAAGCGTTTCGCTGGGAAATAACGTGTAGCCGgaga                         <  2:5330034/100‑1 (MQ=255)
                                                                        cACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGCAAAAACGTGTATCCGGagag                        >  2:3899834/1‑100 (MQ=255)
                                                                        cACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGagag                        <  2:853586/100‑1 (MQ=255)
                                                                        cACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGagag                        >  1:1909528/1‑100 (MQ=255)
                                                                        cACGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGagag                        >  2:1405427/1‑100 (MQ=255)
                                                                         aCGATTAGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCGTTTTATAGAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAgg                       <  1:5292549/100‑1 (MQ=255)
                                                                               aGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGatcatc                 >  1:2943198/1‑100 (MQ=255)
                                                                               aGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGatcatc                 <  2:2882430/100‑1 (MQ=255)
                                                                               aGATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGatcatc                 <  2:4332606/100‑1 (MQ=255)
                                                                                gATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCt                >  1:2062155/1‑100 (MQ=255)
                                                                                gATAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCt                >  2:49923/1‑100 (MQ=255)
                                                                                 aTAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTg               >  2:2307406/1‑100 (MQ=255)
                                                                                 aTAGCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTg               <  1:1748567/100‑1 (MQ=255)
                                                                                    gCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAAc            <  1:5065790/100‑1 (MQ=255)
                                                                                    gCACAGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAAc            >  1:1584992/1‑100 (MQ=255)
                                                                                      acaGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACgg          >  1:4574512/1‑100 (MQ=255)
                                                                                       caGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGc         >  2:4875362/1‑100 (MQ=255)
                                                                                       caGTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGc         >  1:5479697/1‑100 (MQ=255)
                                                                                         gTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCAt       >  2:3144985/1‑100 (MQ=255)
                                                                                         gTGACCGTGGAACTATTCGTACTACATTTATTGCTTTTTACATTGTTAATCGTTTCGCTGGGAAAAAACGTGTATCCGGAGCGGATCATCTGAACGGCAt       >  2:1850130/1‑100 (MQ=255)
                                                                                           gACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGa     <  1:3594516/100‑1 (MQ=255)
                                                                                           gACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGa     >  2:5141995/1‑100 (MQ=255)
                                                                                           gACCGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGa     <  1:2227620/100‑1 (MQ=255)
                                                                                             ccGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGaaa   <  2:2824833/100‑1 (MQ=255)
                                                                                             ccGTGGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGaaa   >  1:1803864/1‑100 (MQ=255)
                                                                                              cGTGGATCTATTCGTACTACATTTATTCCTTTTTGTATAGTTAAGCGTTTCGCTGTGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGaaaa  <  2:5252515/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATATGGTCACTTCACGATTAGATAGCACAGTGACCGTTGAACTATTCGTACTACATTTATTGCTTTTTATATAGTTAAGCGTTTCGCTGGGAAAAAACGTGTATCCGGAGAGGATCATCTGAACGGCATGAAAA  >  NC_000913/563949‑564142

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: