Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,644 C→T S33N (AGC→AAC)  yjdN ← PF06983 family protein YjdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,6440CT93.8% 150.7 / NA 48S33N (AGC→AAC) yjdNPF06983 family protein YjdN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/24);  minor base G (0/3);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGG  >  NC_000913/4325403‑4325890
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
ttCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTtat                                                                                                                                                                                                                                                 >  1:833/1‑251 (MQ=255)
         tCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCAGTCCTGAGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGtt                                                                                                                                                                                                                                        >  2:435185/1‑251 (MQ=255)
           cggtAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCAGCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATGCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATGTTATAGAGCAGTTcc                                                                                                                                                                                                                                      <  2:30351/248‑1 (MQ=255)
                  ggCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCCCCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCa                                                                                                                                                                                                                               <  1:367631/251‑1 (MQ=255)
                  ggCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGCTGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGCATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGccc                                                                                                                                                                                                                                >  1:409260/1‑250 (MQ=255)
                              aaaCCAGCTGTTTCCTTCTTCGACCCGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACATACGTTGat                                                                                                                                                                                                                   <  2:89479/251‑1 (MQ=255)
                                        tttCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCaa                                                                                                                                                                                                          <  1:143866/250‑1 (MQ=255)
                                             ttcttcGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAACTTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATACTAGGCAATCGCGt                                                                                                                                                                                                    >  2:110394/1‑251 (MQ=255)
                                             ttcttcGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGt                                                                                                                                                                                                    >  1:158733/1‑251 (MQ=255)
                                                   gACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAAc                                                                                                                                                                                              >  1:1651/1‑251 (MQ=255)
                                                    aCCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAAc                                                                                                                                                                                              >  1:64690/1‑250 (MQ=255)
                                                            cGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGATCAGTTACCgg                                                                                                                                                                                     <  2:312334/251‑1 (MQ=255)
                                                               aTCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCa                                                                                                                                                                                   >  1:321940/1‑250 (MQ=255)
                                                                  gAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAag                                                                                                                                                                               <  1:201311/251‑1 (MQ=255)
                                                                   aGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCACTCGCGTCGGAACAGTTACCGGCAAAaga                                                                                                                                                                              <  1:235300/251‑1 (MQ=255)
                                                                        cAGCGTAAAGCCGTAGTAGCTGGCTTTTCCGGACGGCATGGCATCGCTCATCATGGTGTCGCTGCCGGCAATGGGCACATGGGGGTGTGTGATGGCGGTATCGGGAAATTGCATTCCGGGAGGGCAGTGCTCGGCGCGGTCCTGCGCTGATTATGGGATTGCGCCGATGGTGAGTTTATAGTGTGGTGCCGGGCCTTTCGTACGTTGAGAATGGGCAATGGCGGGGGAACATTTACCTGCAGAATAGAGGt                                                                                                                                                                         <  2:176400/251‑1 (MQ=255)
                                                                                 gCCGGAGTAGCTGGCTTTTCCTGACGGCGTGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTa                                                                                                                                                                <  2:258346/251‑1 (MQ=255)
                                                                                 gCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTa                                                                                                                                                                <  1:227799/251‑1 (MQ=255)
                                                                                        tAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAAt                                                                                                                                                                                                         <  1:92854/203‑1 (MQ=255)
                                                                                        tAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAAt                                                                                                                                                                                                         >  2:92854/1‑203 (MQ=255)
                                                                                                         cGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCAGTTCGCCGAAGTTGATGTTATAGAGCAGTTCCGCGCCCAACGGACGTTGATAATAGGCTATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCTTGA‑TGGTGACCTCAGTTa                                                                                                                                       <  1:152256/251‑1 (MQ=255)
                                                                                                          ggcatggcatCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGt                                                                                                                                                                                           >  1:389228/1‑199 (MQ=255)
                                                                                                          ggcatggcatCGCTCATCAGGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCTATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGt                                                                                                                                                                                           <  2:389228/199‑1 (MQ=255)
                                                                                                                     gCTCATCAGGATGTCGCGTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑t                                                                                                                          <  2:151712/251‑1 (MQ=255)
                                                                                                                        catcatGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGTCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGGTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTa                                                                                                                        <  1:148319/250‑1 (MQ=255)
                                                                                                                        catcatGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAAGTGCATTCCGGAAGGGCAGTTCTCGTCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGTAGTTGATTTTATAGGGCAGTTCCGCGCCCAGCGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAATAGAGGCATGGACTTAACGGCAGGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑tt                                                                                                                         <  1:270535/249‑1 (MQ=255)
                                                                                                                          tcatGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAgg                                                                                                                     >  1:357883/1‑250 (MQ=255)
                                                                                                                                  gcGCTTCCGGCAATGCGCACATGGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCg                                                                                                             <  2:29192/250‑1 (MQ=255)
                                                                                                                                                          gcatgagCGATGGCGGAATCGGGAAATTGCATTCCGGACGGGCAGTTCTCGGCGCTGTCGTGCGCTGATTTTGGCATTTCTCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTGACCGTCAAAAGAGGGGTCGGGAGCCAACGGCATGA‑CGGTGAACTCAGTTAAGGGACGCCAGGTTAGGT‑TAGCTCGAAATTAATACAGATGTGcggcggg                                                                                     >  2:16459/7‑250 (MQ=255)
                                                                                                                                                           cGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGcc                                                                                    <  1:191823/251‑1 (MQ=255)
                                                                                                                                                                              gggAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGt                                                                 >  1:80841/1‑251 (MQ=255)
                                                                                                                                                                                gAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCg                                                               >  1:125706/1‑251 (MQ=255)
                                                                                                                                                                                  aaTTGCATTCCGGAAGGGCTGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGCAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGTACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGCGAAGCCAG‑TTAAGTGTAGTTCGGAATTTATACAGATGAGAGGCGGCCTGATAAGAcgcgttagcgacgc                                                              >  2:256302/1‑240 (MQ=255)
                                                                                                                                                                                   aTTGCATTCCGGCAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGCTTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGCGGAAGGGACTTAACGGCATGA‑TGGTGACCTCCGTTAAGAGAAGCCAG‑TTAAGGGTAGTTCGACATTTATACAGATGAGAGGCGGCCTGATacgaaacccgagcgtagcat                                                            >  2:316378/1‑232 (MQ=255)
                                                                                                                                                                                     tGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGACCAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGGCGCATc                                                           <  1:299181/250‑1 (MQ=255)
                                                                                                                                                                                         ttCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCACAAGAGAGGTAGGGACTTAACGGCATGA‑TGGAGACCTCAGTTAAGAGCAGCCAG‑TTAAGTGTAGTTCGAAATTTAAACAGATGAGAGGCGGCCTGATAAGACGCGAGCGCGTCGccaccggc                                                       >  2:254572/1‑243 (MQ=255)
                                                                                                                                                                                            cGGAAGTGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGGTGATTTTATAGAGCAGTTCCGCGCGCGTCGTACGTTGATAATAGGCAGTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGAGTGG‑GACCGCAGGTAAGAGGAGCCAG‑TTGAGTTTAGTTCTAAAGTTATAGAGAGGAGAGGCTGCGGGATAAGACGCGCGAGCGTCGCATCAGGCAGTc                                                   <  1:408136/251‑1 (MQ=255)
                                                                                                                                                                                            cGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGGTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGTCCTGAGAAGACGCGCGAGCGTCGCATCAGGCAGTc                                                   <  2:434549/251‑1 (MQ=255)
                                                                                                                                                                                               aaGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGc                                                <  1:237372/251‑1 (MQ=255)
                                                                                                                                                                                                               gcTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcg                                <  2:295192/251‑1 (MQ=255)
                                                                                                                                                                                                                 tGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGTGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGGCGCATCAGGCAGTCGGCACTGGTGCCGGATgcggc                              >  2:279179/1‑251 (MQ=255)
                                                                                                                                                                                                                  gcccTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcgg                               <  2:245761/247‑1 (MQ=255)
                                                                                                                                                                                                                             aTTTTGGCATGTCTCCGAAGTTGAGTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAGGAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTGAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACTCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCGGGATGCGGGGTAAActactta                  <  2:410617/251‑7 (MQ=255)
                                                                                                                                                                                                                               tttGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAt                 <  2:158733/250‑1 (MQ=255)
                                                                                                                                                                                                                                ttGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCTGAACAGTTACCGGCAAAAGAGAGTTAGGGACGTACCGGCatg                                                                                                                                                        <  2:178146/116‑1 (MQ=255)
                                                                                                                                                                                                                                ttGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCatg                                                                                                                                                        >  1:178146/1‑116 (MQ=255)
                                                                                                                                                                                                                                         cGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACgg      >  1:133097/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGGAGGGACTTAACGGCACGA‑TGGTGACTTCAGTTAAGCGCAGCCAG‑TTACGTGTAGTTCGAAATTTATACAGCTGAGAGGCGGACTGATCAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGCTGAGGCGTAACCGCCTTATCCGGCCTACGGCggg  >  1:439008/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
TTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGA‑TGGTGACCTCAGTTAAGAGAAGCCAG‑TTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGG  >  NC_000913/4325403‑4325890

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