Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 76,000 T→C K434E (AAA→GAA)  sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091376,0000TC97.5% 120.5 / ‑4.1 41K434E (AAA→GAA) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (19/20);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

ATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTG  >  NC_000913/75754‑76244
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
gacgttgccccatgtgagcgtgaccaggttatgttttggcagcgccaggttggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAg                                                                                                                                                                                                                                                   <  1:155069‑M2/81‑1 (MQ=255)
 acgttgccccatgtgagcgtgaccaggttatgttttggcagcgccaggttggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCg                                                                                                                                                                                                                                                 >  1:83949‑M2/169‑251 (MQ=255)
      gccccatgtgagcgtgaccaggttatgttttggcagcgccaggttggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGAcc                                                                                                                                                                                                                                            <  1:308699‑M2/88‑1 (MQ=255)
          catgtgagcgtgaccaggttatgttttggcagcgccaggttggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttctccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGTCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCtggg                                                                                                                                                                                                                                        >  2:429835‑M2/160‑250 (MQ=255)
                       ccaggttatgctttggcagcgccaggttggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCCGATATCGCTCTCGCTCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAAt                                                                                                                                                                                                                           >  2:199667‑M2/147‑251 (MQ=255)
                              gTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGACAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAACTGTGCGAGAACAGAGAAGTCCAGCGGCAGGGTAAAGATGGCGCTGTTTAGCCAGATATCACTTTCGATCTATCCTGTATGCCACTGATCGTAGTCGATCTATTTTGTTTTCCGCGGTACCTGGCGACATGACAGAATCtgatgc                                                                                                                                                                                                                    >  2:137933/1‑251 (MQ=255)
                                        agcgccaggttggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCgg                                                                                                                                                                                                          <  2:110375‑M2/122‑1 (MQ=255)
                                            ccaggttggcctctaatacctggcgtttgagatcttctaacatggtgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAAt                                                                                                                                                                                                      <  2:157992‑M2/126‑1 (MQ=255)
                                                   ggcctctaattcctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGa                                                                                                                                                                                               <  1:180700‑M2/133‑1 (MQ=255)
                                                   ggcctctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCg                                                                                                                                                                                                <  1:445770‑M2/132‑1 (MQ=255)
                                                     cctctaatacctggggtttgatatcttctaacattttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTGAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGAGCCCGGCAATCACCCGATg                                                                                                                                                                                             <  2:157203‑M2/135‑1 (MQ=255)
                                                        ctaatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGCTCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGCTTTCCAGCGTGCCCTGGTGCCTTGCCAGAATCTGCTGCATGATCCCGGCACTCCCCCGATGCTc                                                                                                                                                                                          >  2:195628‑M2/114‑251 (MQ=255)
                                                          aatacctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCt                                                                                                                                                                                           <  1:282039‑M2/137‑1 (MQ=255)
                                                             acctggcgtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACtgt                                                                                                                                                                                     <  1:412342‑M2/143‑1 (MQ=255)
                                                              cctggggtttgagatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACtgtg                                                                                                                                                                                    <  2:83949‑M2/144‑1 (MQ=255)
                                                                         agatcttctaacatgttgactccttcgtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaa                                                                                                                                                                          >  2:178143‑M2/97‑250 (MQ=255)
                                                                                                   gtgccggatgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTTCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGcc                                                                                                                                                >  1:430590‑M2/71‑250 (MQ=255)
                                                                                                       cggatgcgctttgcttatccggcctacaaaagcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCg                                                                                                                                           <  1:80916‑M2/185‑1 (MQ=255)
                                                                                                           tgcgctttgcttatccggcctacaaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGATGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGAAAAAAGGTCAGGAGGAGGCTTTCCAGGCCAGCCGGttt                                                                                                                                       >  2:271899‑M2/63‑251 (MQ=255)
                                                                                                                       atccggcctacaaaagcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTTTTAAGCCAGATATCTCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTTTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTa                                                                                                                            <  2:293713‑M2/200‑1 (MQ=255)
                                                                                                                                 caaaatcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGCTTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCGTGATCCCGTCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGCTGAGGCTTACCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGAGGGc                                                                                                                  >  1:139979‑M2/41‑250 (MQ=255)
                                                                                                                                     atcgcagcgtgtaggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGACAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGATCTTTATTGTGCGGGCAtggg                                                                                                             >  2:299426‑M2/37‑250 (MQ=255)
                                                                                                                                                 aggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACggggg                                                                                                  >  1:171929‑M2/25‑250 (MQ=255)
                                                                                                                                                  ggcctgataagacgcgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGATTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACggggg                                                                                                  >  1:82413‑M2/24‑249 (MQ=255)
                                                                                                                                                            gacgcgacagcgtAAAGTTGGCGCTGTTAAGCCATATATCGCGCTCGATCGCCCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTGTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATaa                                                                                      <  1:328971‑M2/238‑1 (MQ=255)
                                                                                                                                                                cgccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCATCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCACCGGGGGAGCAGTCCATAAGCCg                                                                                  >  2:11980‑M2/10‑251 (MQ=255)
                                                                                                                                                                  ccagcgtAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                                                                                <  1:178143‑M2/244‑1 (MQ=255)
                                                                                                                                                                      cgtAAAGTTGGCGCTGTTAAGCTAGTTATCGCTCCCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAAcc                                                                            <  1:106129‑M2/248‑1 (MQ=255)
                                                                                                                                                                            gTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGtt                                                                      <  1:195628/251‑1 (MQ=255)
                                                                                                                                                                                 cgcTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTTCATGAGCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTgg                                                                 <  1:16647/251‑1 (MQ=255)
                                                                                                                                                                                               atatCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGc                                                   >  1:246735/1‑251 (MQ=255)
                                                                                                                                                                                                             ctctcCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTc                                                                                                                                    <  2:175054/156‑1 (MQ=255)
                                                                                                                                                                                                             ctctcCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTc                                                                                                                                    >  1:175054/1‑156 (MQ=255)
                                                                                                                                                                                                                  cttcaTGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGtt                                >  1:223009/5‑251 (MQ=255)
                                                                                                                                                                                                                  cttcaTGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGtt                                >  1:255400/5‑251 (MQ=255)
                                                                                                                                                                                                                        gcctcTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTTGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGaa                          <  2:88820/247‑1 (MQ=255)
                                                                                                                                                                                                                        gCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGaa                          <  1:429835/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tagctgATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGCACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCt               >  2:72597/6‑251 (MQ=255)
                                                                                                                                                                                                                                        gATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCCGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGCAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTAccc          >  2:65944/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGAACAGCTTAGTTGGAGAAAGAACATAGCTTACCCAGTCCCTg  >  2:320847/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gATTTGCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTg  <  1:137933/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTG  >  NC_000913/75754‑76244

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