Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,332,871 A→C 13.3% S134R (AGT→AGG)  basS ← sensory histidine kinase in two‑component regulatory system with BasR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,8710AC13.3% 151.6 / 5.5 72S134R (AGT→AGG) basSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base C (4/7);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGAT  >  NC_000913/4332724‑4333019
                                                                                                                                                   |                                                                                                                                                    
aaGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAACCAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                                                                                                     >  2:703327/1‑149 (MQ=255)
  gCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGCCTg                                                                                                                                                   >  2:812225/1‑149 (MQ=255)
  gCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                                                                                                   >  2:123129/1‑149 (MQ=255)
    cgggCCACAAACGGTGCTACAACAAGGTGTGGCGGTAGCGCCAGTAGGTCCAGATGCACACGCACCCCCGCCAGTGGCGTTCGCAGGTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTCTCCAGCGTGCTGGTCAGGCGACTGAc                                                                                                                                                 <  2:965714/146‑1 (MQ=255)
     gTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                                                                                                >  2:611947/1‑149 (MQ=255)
      tGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGCGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCa                                                                                                                                               >  1:102440/1‑149 (MQ=255)
      gcccACTAACGGTGCGACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTGCGCAGTTCGTGCGCGACGGCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCa                                                                                                                                               <  1:806040/147‑1 (MQ=255)
       gCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                                                                                                              >  1:229833/1‑149 (MQ=255)
       gCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                                                                                                              >  1:233630/1‑149 (MQ=255)
          aCTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                                                                                                           <  1:289323/149‑1 (MQ=255)
            tAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                                                                                                         <  2:404710/149‑1 (MQ=255)
             aaCGGTGCTACATCAATGTGCTGCGTTTTCGCCTGCAGTTCCAGATGGAATCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcacgct                                                                                                                                        <  1:1017061/149‑7 (MQ=255)
              actgTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCTCCCCCGCCAGTGGCTTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTGAGGCGACTGACCAGATCGttt                                                                                                                                       <  2:739489/146‑1 (MQ=255)
              aCGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcacgctg                                                                                                                                       <  1:501279/149‑8 (MQ=255)
                   gCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTGCGCAGGTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                                                                                                                  <  1:602608/149‑1 (MQ=255)
                       cATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTCAACGCCGaa                                                                                                                              >  2:680958/1‑149 (MQ=255)
                          cAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                                                                                                           >  1:2564/1‑149 (MQ=255)
                              gtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCa                                                                                                                                               >  1:1006025/1‑125 (MQ=255)
                              gtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCa                                                                                                                                               <  2:1006025/125‑1 (MQ=255)
                                gATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCt                                                                                                                     >  2:192234/1‑149 (MQ=255)
                                gATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                                                                                                                     >  2:605789/1‑149 (MQ=255)
                                 acgCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                                                                                                    <  1:234681/147‑1 (MQ=255)
                                 aTGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                                                                                                    >  1:774745/1‑149 (MQ=255)
                                 aTGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                                                                                                    <  2:42881/149‑1 (MQ=255)
                                 aTGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCATGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCtt                                                                                                                    >  1:906635/1‑149 (MQ=255)
                                  tGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCTGTGGCGTTCGCAGTTCGTGCGCGACGTCAGGGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCTCCGCTTc                                                                                                                   <  1:1028815/149‑1 (MQ=255)
                                      ttttCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACCACCTTTCGTAATCCAGCCTTCTCGTCAGGAGACTGACCAGATCGTTTAAAGCAGAAACCACCGCTTCGAtt                                                                                                               >  2:181423/1‑149 (MQ=255)
                                            ccAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                                                                                                         <  2:233630/149‑1 (MQ=255)
                                            ccAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                                                                                                         <  2:229833/149‑1 (MQ=255)
                                            ccAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGGACACATCCTGTCGTTATCCAGCTTGCTGGGAAGGCGACTGCCCAGAGCGTTTACCGCCGAAACCTCGGCTTCGATTTCGCgg                                                                                                         >  1:335134/1‑149 (MQ=255)
                                                     ccAGATGCAAACGCACCCCCGCCAGTGGTGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGGtg                                                                                                >  1:253036/1‑149 (MQ=255)
                                                     ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                                                                                                >  2:371087/1‑149 (MQ=255)
                                                      cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt                                                                                               >  2:68339/1‑149 (MQ=255)
                                                      cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt                                                                                               >  2:333173/1‑149 (MQ=255)
                                                      cAGATGCAAACGCACCCACGCGAGTGGCGTGGGCAGTTCGTGCGCGACGGCCGCGGTAAACAACCTTTCGTTATCCAGCGGGCTGGTCAGGCGACTGTCCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGAGCtgt                                                                                               <  2:587578/149‑1 (MQ=255)
                                                           gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCTTCGATTTCGAGGGTGGCGCTGTGAATg                                                                                          >  2:133588/1‑149 (MQ=255)
                                                           gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGCCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCAGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg                                                                                          >  2:737714/1‑149 (MQ=255)
                                                             aaaCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                                                                        <  2:931560/149‑1 (MQ=255)
                                                             aaaCGCACCCCCGCAAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTAGTTATCCCTAGTGCTTGTCAGGAGACTGCACAGATCGTTTAAAGCCGAAACCCCCGCAGCGAGTTCGAGGGTGGCTCTGTGAATGGc                                                                                        >  2:130713/1‑149 (MQ=255)
                                                                cGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTCAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt                                                                                     <  1:680958/149‑1 (MQ=255)
                                                                  cACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGTCGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATgg                                                                                   <  1:690759/149‑1 (MQ=255)
                                                                  cACCCCCGCCAGAGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATgg                                                                                   <  2:603956/149‑1 (MQ=255)
                                                                   aCCCCCGCCAGTGGCGTTCGCAGTTCTGGCGCGACGTCATCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGAGACTGAGCTGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATggg                                                                                  <  2:951600/149‑1 (MQ=255)
                                                                   aCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATggg                                                                                  >  1:927077/1‑149 (MQ=255)
                                                                     ccccGCCAGTGGCGTTCGCAGTTCGTGCGCGCCGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCCGGCGACTGACCAGCTCGTTTAACGCCGACACCACCGCTGCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCg                                                                                >  2:686848/1‑149 (MQ=255)
                                                                       ccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGtt                                                                              >  2:169530/1‑149 (MQ=255)
                                                                           cAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGt                                                                          >  2:104978/1‑149 (MQ=255)
                                                                            aGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATATCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGtt                                                                         <  1:621422/149‑1 (MQ=255)
                                                                                         tttCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc                                                            <  2:638521/148‑1 (MQ=255)
                                                                                         gTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc                                                            >  2:356346/1‑149 (MQ=255)
                                                                                         gTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc                                                            <  2:793930/149‑1 (MQ=255)
                                                                                             gTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCtt                                                        <  1:826408/149‑1 (MQ=255)
                                                                                             gTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGCGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCtt                                                        <  2:853898/149‑1 (MQ=255)
                                                                                               gcgcGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTcc                                                      >  1:945525/1‑149 (MQ=255)
                                                                                                 gcgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAg                                                    >  1:932245/1‑149 (MQ=255)
                                                                                                 gcgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGATTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAg                                                    <  1:17802/149‑1 (MQ=255)
                                                                                                  cgACGTCAGCGGTAAACAACCTTTCGTTATCCACCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGTTTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGc                                                   <  2:653238/149‑1 (MQ=255)
                                                                                                   gACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGct                                                  >  2:313558/1‑149 (MQ=255)
                                                                                                         aGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTACCGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCTATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCttttt                                            >  1:92431/1‑149 (MQ=255)
                                                                                                         aGCGGTAAACAAACTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGTCCAGATCGTTTAACGCCGAAACCTCCGCTTCGATTTCGAGGGGGGCGCTGTCAATGGCAATGGGCGGTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCttttt                                            <  2:235289/149‑1 (MQ=255)
                                                                                                          gCGGTAAACAACCTTGCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTg                                           <  1:703327/149‑1 (MQ=255)
                                                                                                           cGGTAAACAACCTTTCGTTATCCTGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGc                                          <  2:702460/149‑1 (MQ=255)
                                                                                                           cGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGc                                          <  1:868198/149‑1 (MQ=255)
                                                                                                            ggTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCGCTTTTTGCa                                         >  1:867671/1‑149 (MQ=255)
                                                                                                                  caacCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCTATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCCGCTCTTTTTGCAGCTCCg                                   >  2:570518/1‑149 (MQ=255)
                                                                                                                   aacCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTACCGCCGAAACCACCGCTGCGATTTCGAGGGTGGCGCTGTGAATGGCATTGGGCGTTAAGTTGTCGGCGGTGCGCGCGTCCAGCTCTTTTTGCAGCTCCGc                                  >  2:695588/1‑149 (MQ=255)
                                                                                                                             tATCCCGCGTTCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcgggcgg                        <  2:799591/149‑1 (MQ=255)
                                                                                                                             tATCCAGCGGGCTGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTGGATGGTGGCGCTGTGAATGGCTATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcgggcgg                        <  1:790611/149‑1 (MQ=255)
                                                                                                                               tCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTGCGATTTCGAGGGTGGCGCTGTGAATGGCTATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGt                      <  1:716378/149‑1 (MQ=255)
                                                                                                                                  aGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGCTTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTGAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCGCCGCCAGCGGGCGGGTGAt                   >  2:257718/1‑149 (MQ=255)
                                                                                                                                    cGTGCTGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGTGGGCGGGTGATgc                 <  2:774745/149‑1 (MQ=255)
                                                                                                                                      ggCTGGTCAGGCGCCTGACCAGATCCTCTACCGCCGAACCCCCCGCTTCCATCTCGAGGGTGGCGCTGTGAAGGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATgcgg               <  1:902705/148‑1 (MQ=255)
                                                                                                                                                   aCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGGGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGAt  <  2:906635/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGAT  >  NC_000913/4332724‑4333019

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑