Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,697 C→T F334F (TTC→TTT)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,6970CT97.4% 109.5 / NA 39F334F (TTC→TTT) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/23);  minor base G (0/1);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCG  >  NC_000913/4637460‑4637940
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
gattgagccggattgttaagcgttAACCCATGAGGTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGTAAATTTTACGCGAAGGGCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGGGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGc                                                                                                                                                                                                                                         <  2:208558/228‑1 (MQ=255)
     agcagttggtttatgcggtAACCCATGAGCTAAATAGCCCACTTGCGGCGATTCGTTGCGCGGCGGAAATTTTACGCGAAGGTCCTCCACCGGAAGTGGTGGCTCGTGTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTGCTGACTTTTGGTGAAGTGGTGGCATTATATCGCCGCGTCAGCGAAgcg                                                                                                                                                                                                                                   <  2:362908/233‑1 (MQ=255)
         gTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGGCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGGGAATCGTAAGGAAGTCGTTCTGCCTGTTGTTGATGTGGCGGCATTATTTGGCCGCGGCACCGTAgcgcgcc                                                                                                                                                                                                                               >  2:138834/1‑250 (MQ=255)
             ggttatgcgttAACCCATGAGCTAAAAAGCCCACGGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCg                                                                                                                                                                                                                            <  2:259985/241‑1 (MQ=255)
              tttatgcgttAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGGGAAgagggccccggg                                                                                                                                                                                                                          >  1:40221/9‑240 (MQ=255)
                         acccATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGGTCTGACTGTTGTTGATGTGGCGGCAATATTTCGCCGCGTCAGCGAAGCGCGCCCCGAGCAGTTGGcgg                                                                                                                                                                                                                >  1:393973/4‑248 (MQ=255)
                          cccATGAGCTAAAAAGCCCACTGGCGGCGTGTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTTCCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATGATTTCGCCGCGTCAGCGAATCGCGCACCTTGCAGTTGGCAGaa                                                                                                                                                                                                              <  2:331663/249‑1 (MQ=255)
                           ccATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGaa                                                                                                                                                                                                              >  1:343515/2‑250 (MQ=255)
                                aGCTAAAAAGCCCACTGGCGGCGATTCGTTGCGCGGCGGAAATTTTACGCGAAGGGCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATGTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGaaaaaaac                                                                                                                                                                                                        <  2:251783/251‑2 (MQ=255)
                                      aaaGCCCACTTGCGGCGAGTCGTGGCGTGGCGGAAATTTTACTCGAAGGGCCGCCACCGGAAGGGGGGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGGGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGTCAGAAAAAAACATCAct                                                                                                                                                                                                  <  2:339571/251‑1 (MQ=255)
                                           ccACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGGAGTGGGGGCTCGTTTTACCGACAACATTCTGACGCAAAGTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGCTGTGGCGGCGATCTTTTGCCGCGTCCGCGAGCCGCGCCCCGTGCAGATGGCCGgaaaaaaccacaatctacc                                                                                                                                                                                             >  2:132348/1‑239 (MQ=255)
                                                  cggcgCTTCGGGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGCCAACATTCTGACGAAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTACCTGTTGTTGCTGTGGCGGCATTATTTCGACGCGGCCGCGGAGAGCGACCCGGGCCGTTGGCAGAGAAAAAAATCAatctcacagttatg                                                                                                                                                                                      >  2:503816/1‑241 (MQ=255)
                                                               cgcggcgGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTa                                                                                                                                                                         <  1:412230/250‑1 (MQ=255)
                                                                     ggAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTtg                                                                                                                                                                   <  1:378929/251‑1 (MQ=255)
                                                                            ttACGCGAAGGGTCGCCACCGAAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTTGAGAAGAGTCAGGAAGTGGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTTAGCGATGCGCGCACCGTGCAGGTGGCAGAAAAAAACAGCACTCGACATGTTATGCCCACTGAGGTTAATGTTGctgctg                                                                                                                                                             <  2:288022/250‑1 (MQ=255)
                                                                              aCGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGTAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGaa                                                                                                                                                           <  1:205956/250‑1 (MQ=255)
                                                                                                             cGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAAttt                                                                                                                           <  2:238339/251‑1 (MQ=255)
                                                                                                                               cTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACCTGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCg                                                                                                          >  2:192092/1‑250 (MQ=255)
                                                                                                                                  aCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTATTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGGCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGAttt                                                                                                      <  1:400921/251‑1 (MQ=255)
                                                                                                                                   cGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGAtttt                                                                                                     >  1:376968/1‑251 (MQ=255)
                                                                                                                                        aaTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAcccc                                                                                                <  2:376968/251‑1 (MQ=255)
                                                                                                                                               gTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGc                                                                                         <  1:358173/251‑1 (MQ=255)
                                                                                                                                                 aTGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTaatg                                                                                                                                                                      <  2:461836/172‑3 (MQ=255)
                                                                                                                                                    cAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGtcg                                                                                    <  1:188987/251‑3 (MQ=255)
                                                                                                                                                                  aaCGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCCGTTGGCAGAAAAAAACATTAATCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCggacgcataaagcttagcg                                                                      >  1:310390/1‑234 (MQ=255)
                                                                                                                                                                   aCGTTACTACGCCAGGCAAGACTGGAGAAGCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGTACCGTGCAGTTGGCAGAAAAGAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCAGCGATTTTACCCCCGAGAGCGTTCGCATAACGCTAAgcgc                                                                     <  1:503816/251‑1 (MQ=255)
                                                                                                                                                                           acGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGAATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTgg                                                             >  1:246534/1‑251 (MQ=255)
                                                                                                                                                                              ccAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATc                                                          <  1:496245/251‑1 (MQ=255)
                                                                                                                                                                               cAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCa                                                         <  2:343515/251‑1 (MQ=255)
                                                                                                                                                                                           gagaATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGGCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTcgc                                             >  2:29242/1‑249 (MQ=255)
                                                                                                                                                                                             gaATCGACAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAACCATCACTCTACATGGTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGGTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGGCGCATAACGCTAAGCGCCGAAGGGGGTCAGGAACACGTcgccc                                           >  1:24738/1‑247 (MQ=255)
                                                                                                                                                                                                 cGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTCTTTCGCCGCGTCAGCGCAGCGCGCACCGTGCAGTTGGCAGAAAAAAACCTCCCTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGGGCAGGCGCTGGGGAAATTACTGGGTAACGccctc                                                                                                           >  2:231594/1‑180 (MQ=255)
                                                                                                                                                                                                 cGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATc                                                                                                           <  1:231594/183‑1 (MQ=255)
                                                                                                                                                                                                        aaGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTg                                <  2:407711/251‑1 (MQ=255)
                                                                                                                                                                                                          gTCGTTCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGGTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTAATGGATAACGCTATCGTTTTTACCCCCGAGAGCGGTCGCATAACGCTTAGCGCCGGAGTGGTTCAGGacaacgtcgaccttcaggggctggg                              >  2:299495/1‑227 (MQ=255)
                                                                                                                                                                                                            cattCTTACTGTTGTTGATGTGGCGGAATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTGGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGAt                             <  2:365365/248‑1 (MQ=255)
                                                                                                                                                                                                                         gtggATGTGGCGGCATTATTGCGCCGCGTCAGCGATGCGCGCACCGTGCAGTTGGCTGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATATCGTTAGTGGTa                <  2:326377/247‑1 (MQ=255)
                                                                                                                                                                                                                                     gcATTATTTCGCCGCGGCAGCGACGCGCGCACCGTGCAGCTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAGGCTAAGCGCCGACGTGGATCAGGAACACGTCGCCCTTCAGGTGCGGGATACCGGGAGTGGTATTCCTGACTTcgc   >  2:37348/1‑251 (MQ=255)
                                                                                                                                                                                                                                      caatattTCGCCGCGTAAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTAGGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTCAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTAcgcg  <  2:376046/248‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
TATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCG  >  NC_000913/4637460‑4637940

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