Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,077,190 C→T Q89Q (CAG→CAA yggP ← putative zinc‑binding dehydrogenase YggP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,077,1900CT95.8% 144.9 / ‑5.9 50Q89Q (CAG→CAAyggPputative zinc‑binding dehydrogenase YggP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (20/26);  minor base A (1/1);  total (22/28)
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

CCATGCGATGCTCATAAACATATTGCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGC  >  NC_000913/3076946‑3077428
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
acattcGATGCTCATAAACATATTGCGTGGGGGGATAATTGGCATGATAAGAACCAATAATGAAGCACATAGGATCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATAAGGTAATACGCAGCCCAAATTATTGACGATTTCGGGAATAATCATATAAGTAGCATTGTCGATAAAATATTCGTCGCTGTAGCCCCATGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAAAATAAAACGATGGctt                                                                                                                                                                                                                                          >  2:230774/6‑249 (MQ=255)
  aTGCGATGCACATAAACATATTGCGTGGGGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGGAACCCCATCGCCGGTTGCAATACAAAACGTTGGACttt                                                                                                                                                                                                                                        <  2:80424/251‑1 (MQ=255)
       aTGCTCATAAACATATTGCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATAGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTtat                                                                                                                                                                                                                                   >  2:118554/1‑251 (MQ=255)
        tGCTCATAAACATATTGAGTGGTGTGATAATTGGCATGCTAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTtat                                                                                                                                                                                                                                   >  2:275854/1‑250 (MQ=255)
         gCTCATAAACATATTGCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTtata                                                                                                                                                                                                                                  <  1:3626/250‑1 (MQ=255)
            cATAAACATATTGCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTg                                                                                                                                                                                                                              <  1:96802/251‑1 (MQ=255)
            cATAAAAATATTGCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGTTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTg                                                                                                                                                                                                                              <  2:95385/251‑1 (MQ=255)
                       tGCGTGGTGTTATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGAt                                                                                                                                                                                                                    <  1:382852/250‑1 (MQ=255)
                        gCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGAttt                                                                                                                                                                                                                  <  1:139097/251‑1 (MQ=255)
                         cGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAAGACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGAtttt                                                                                                                                                                                                                 <  2:398753/251‑1 (MQ=255)
                             gtgagATAATTGGCATGCTAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCACAATAAGCGCCGTGATACGGTGATACGCAGCCCAAATTAAGGGCGATTTCGGGAATAATCATATAAGTTGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTAcc                                                                                                                                                                                                             >  2:402405/5‑251 (MQ=255)
                               gtgATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCa                                                                                                                                                                                                           <  1:59651/251‑1 (MQ=255)
                                                        ataatGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAtcc                                                                                                                                                                                  <  1:301688/251‑3 (MQ=255)
                                                        ataatGCAGCACATAGGCTCTGCCAGCGAAGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGAATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACACCGGCACCTTc                                                                                                                                                                                  >  1:227615/1‑251 (MQ=255)
                                                          aaggcagcaCATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAAGTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATAGCCGCTTTGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCTCTTCGACAATAACCCCGGCACATTCAt                                                                                                                                                                                <  2:240619/248‑1 (MQ=255)
                                                           atGCAGCACATACGATCTACCAGCGAAGCACCACCAAAATAAGAGCCGTGATCCGGTAATACGCAGCCCAAATTAATGCCGATTTAGGGAATAATCATATACGTGGCATTGCCGCCAAAATCTTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATAGACGGTTGCAATTAAAAACGTTGGACTTTTTTATATTTGCCAGTGCGATTTTTACCCAATGTGACAATAACCCCTGCACATTCATg                                                                                                                                                                               >  2:287097/1‑251 (MQ=255)
                                                               cgcaCATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGAGATACGGTAATACGCAGCCCAAATTAATGGCGATGTCGGGCGTAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGGAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGACAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCg                                                                                                                                                                           >  2:140375/2‑251 (MQ=255)
                                                                  acaTAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGt                                                                                                                                                                         >  2:14614/1‑250 (MQ=255)
                                                                    aTAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTACTACGCAACCCAAATTAATGGCGATTTCGGGACTGATAATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATCCAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt                                                                                                                                                                      >  1:367702/1‑251 (MQ=255)
                                                                      aGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATACGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGCCAATAACCCCGGCACATTCATGCCCGGTAATg                                                                                                                                                                     >  1:396971/1‑250 (MQ=255)
                                                                      aGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACCTTCATGCCCGGTAATg                                                                                                                                                                     >  2:300253/1‑250 (MQ=255)
                                                                      aGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGCGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg                                                                                                                                                                     >  2:340913/1‑250 (MQ=255)
                                                                                  cGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGTCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtttg                                                                                                                                                        <  2:314472/251‑2 (MQ=255)
                                                                                               aattAAGAGCCGTGAGACGGAAAACCGCAGCCCAACATAATTGCGATTCCGAGAATAATCATATAAGTGGCATGTCCGCCAAAATATGCGGAGCTGAAGCCTGCTGTATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAGCGTTGGCCTTTTTTATATTTGCCAGTTAGATTTTTACCCACTTCGACAATATCCCCGGCACATTCACTCCCGGTAATGAGCGGTTGATTTTCTATATCGACggg                                                                                                                                           <  2:21413/248‑1 (MQ=255)
                                                                                                        ccGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGttt                                                                                                                                  <  2:350175/251‑1 (MQ=255)
                                                                                                          gTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCTGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCTGGTACGCGTTTAt                                                                                                                                <  2:148169/251‑1 (MQ=255)
                                                                                                                    aaTACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTAcc                                                                                                                      >  1:73085/1‑251 (MQ=255)
                                                                                                                       acgtaGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGACGCCAAAATATTCGTAGTTGTAGCCCGCTGAATATTCGTTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAg                                                                                                                   <  1:402405/247‑1 (MQ=255)
                                                                                                                                          cGATTTCGGGAATAAGCATATAAGTGTCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTACCCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTTCCCACTTCGACAATAACCCCGGCACATTCATGGCCGGTAATGACCGGATGATTCTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACTCCGCGTTCAAGGTc                                                                                                <  2:367702/251‑1 (MQ=255)
                                                                                                                                                    aataatCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAAc                                                                                         >  1:262667/1‑248 (MQ=255)
                                                                                                                                                     ataatCATATCAGCAGCACTACCGCAAAAATATTCGTACCTGTAACCCGCTGCATATCCGCTTGGTAACCCCATCGCCGGTTGCAACACAAAACGTTGACCTTTTTTATATTTGCCAGTGAGATATTTACCCACTTCGACACTAACCCCGGCACATTATTGCCCGGTAATGACCGGATGATTTTCTAACTAGTCGGGTACGCGTTTATGTTCACTACCGACAAACGCCGCTTTCCAGGTCGATACACAGAc                                                                                     >  2:391478/1‑251 (MQ=255)
                                                                                                                                                           atatAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTc                                                                               <  1:28041/251‑1 (MQ=255)
                                                                                                                                                                agtgtcatgtccgccAATATATTCGTAGCTGTAGCGCGCTGAATATCCGCTTGTTAACCCCATCGCCGTTTGCAATACAAAACGTTGTCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGTCAATAACCCCGGCACATTCATGCCCTGTAATAACCTGATTATTTTCTAAATCGTCGGGTACGCGTGTATGTTCTCTTCCGATTAACGCCGCGGTCCAGGTCGAGAAACAGACGCTGTCAGaaa                                                                          <  2:373079/241‑1 (MQ=255)
                                                                                                                                                                 gTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAAt                                                                         >  2:263688/1‑251 (MQ=255)
                                                                                                                                                                         gccgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAGCCCAATCGCCGGTTGCAATGCAAAACTTTGGCCTTTTTTATATTTGCAAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAATCATGCCCGGTAATGACCGGATGATTTTCTAAAGCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCTATAAACAGACGCTGTCAGAAATTacacgca                                                                 <  2:311037/251‑4 (MQ=255)
                                                                                                                                                                              cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGt                                                            <  1:14614/251‑1 (MQ=255)
                                                                                                                                                                              cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGt                                                            >  1:185314/1‑251 (MQ=255)
                                                                                                                                                                                 aaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTACATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCGGCCTGAAATTACACTCACCAGCAAt                                                         >  2:417695/1‑251 (MQ=255)
                                                                                                                                                                                   tatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTc                                                       >  1:355312/1‑251 (MQ=255)
                                                                                                                                                                                      tCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCAt                                                     <  2:262667/250‑1 (MQ=255)
                                                                                                                                                                                            cTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGACATTACACTCACCAGTAATTCATTATCGGt                                              >  2:24270/1‑251 (MQ=255)
                                                                                                                                                                                              gTAGCCCGCTGGATATCCGCTAGGTAACCCCATGTCCGTTTGCAATACAAAACGATGGCATTTTTTATATTAGCCAGTGAGATTTTTGCCCACTTCGACAATAACCCCGGCACATACAAGCCCGGTAATGACCGGATGGTTGTCTAAATCGTCGCGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTaa                                            <  1:413261/251‑1 (MQ=255)
                                                                                                                                                                                               gagaccgttGAATATCCGCTTGGTATCCCCATCGTCGGTTGCAATACAAATCGTTGGCCTTTTTGATGTTTGCCAGTGAGAGTGTTACCCACTTCAACAATAACCCCTGCACATTCATGCCCGGTAATTACCGGATGAGTGTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCAAATTCGATAAACAGACGCGGTCAGAAATTTCACTCACCAGTAATTCATTATCGGTAAt                                           <  2:152830/243‑1 (MQ=255)
                                                                                                                                                                                                      cTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATCCAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGAAGCTGTCAGAAATTACACTCACCAGAAATTCATTATCGGTAATTTCTGGc                                    >  2:80955/1‑251 (MQ=255)
                                                                                                                                                                                                            atCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGGTGGTCGTTTTTATATTTGCCAGTGAGAGTTTTACCCACTGCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGAGGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGTAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTc                               <  2:141678/250‑1 (MQ=255)
                                                                                                                                                                                                                    ggTAACCCCATCGCCGGTTGCAATACAAAACGTTGGCATTTTTTATATTTGTCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGTGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGTCAGTTCAAATAcgc                       <  2:310349/250‑1 (MQ=255)
                                                                                                                                                                                                                         ccccATCGCCGGTTGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGGTCAAATACGCGCAGAc                 <  2:391331/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ttGCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCCCTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGATACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCCGGTCGATCAACAGACGCTGTCAGACATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCt     >  1:364938/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gCAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCGCTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCtt    <  1:16514/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cAATACAAAACGTTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGTCAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGc  <  1:321884/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
CCATGCGATGCTCATAAACATATTGCGTGGTGTGATAATTGGCATGATAAGCACCAATAATGCAGCACATAGGCTCTGCCAGCGACGCCGCAGCAAAATAAGAGCCGTGATACGGTAATACGCAGCCCAAATTAATGGCGATTTCGGGAATAATCATATAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTTTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAATCGTCGGGTACGCGTTTATGTTCACTACCGAGTAACGCCGCTTTCCAGGTCGATAAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGC  >  NC_000913/3076946‑3077428

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: