Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 670,395 G→A L60F (CTT→TTT)  nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913670,3950GA98.0% 160.4 / ‑6.3 49L60F (CTT→TTT) nadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/25);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCC  >  NC_000913/670148‑670634
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
gcCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCTTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCGTACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACAt                                                                                                                                                                                                                                               <  2:139888/251‑1 (MQ=255)
gcCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACAt                                                                                                                                                                                                                                               >  1:179145/1‑251 (MQ=255)
                     ggACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTg                                                                                                                                                                                                                            <  2:29502/249‑1 (MQ=255)
                     ggACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGtt                                                                                                                                                                                                                          >  1:86706/1‑251 (MQ=255)
                       aCGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCg                                                                                                                                                                                                                        >  2:93973/1‑251 (MQ=255)
                             aCAGACGATCAAATGTGCATTGTCGAGTATCTTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAACCGCCAGCGGCACGTCCGGTCCTTGTTCCGGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCg                                                                                                                                                                                                                  <  2:179145/251‑1 (MQ=255)
                                   gATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACAAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGTACGTCAGGTCCTTGTTCCTGCCGCCACTCTTTCCGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAATCTGCGTTTACGCTGCACGCTGTTCGCTTCCGGCTggg                                                                                                                                                                                                            >  2:71543/1‑251 (MQ=255)
                                    aTCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTgggg                                                                                                                                                                                                           <  1:136921/251‑1 (MQ=255)
                                           gtgCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACAAATAATAAACGCCAGCGTCACGTCCGGTCTTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGTATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATg                                                                                                                                                                                                     <  1:220673/250‑1 (MQ=255)
                                                 tGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCgg‑aggg                                                                                                                                                                                             >  1:29265/1‑250 (MQ=255)
                                                            gTTTCGTATTCGTACCAGGTCGAAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTGGGCAATAGAC‑GTGTAAAAAATGGGTTTACGCGGCACGGTGTTCGCGTCCGGCTGGGGACG‑CGCGGGAGGAACATAATTAgg                                                                                                                                                                                  >  2:389405/1‑251 (MQ=255)
                                                                         aCCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCGGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGAc                                                                                                                                                                      <  1:326927/250‑1 (MQ=255)
                                                                         aCCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGAc                                                                                                                                                                      <  1:23707/250‑1 (MQ=255)
                                                                            aGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCCGT‑TCAAACATGTGTTTACGCTGCCCGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCcgc                                                                                                                                                                  >  2:336780/1‑251 (MQ=255)
                                                                                    aaGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGAcc                                                                                                                                                          <  1:450106/251‑1 (MQ=255)
                                                                                       gTCAGCAGTGAACCCTGACCAATAATAAACGCCAGCGGCACGTCGGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAAt                                                                                                                                                       <  2:129866/251‑1 (MQ=255)
                                                                                         cagcagTGAATCCGGACCAATAATAAACGCCAGCGGCACGTCTGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCa                                                                                                                                                     <  2:400603/251‑1 (MQ=255)
                                                                                                  aTCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCa                                                                                                                                            <  2:93159/251‑1 (MQ=255)
                                                                                                           aataatAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTcc                                                                                                                                    >  2:451651/1‑250 (MQ=255)
                                                                                                            ataataAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCAc                                                                                                                                  <  2:92879/251‑1 (MQ=255)
                                                                                                            agaaaaaaCGCCAGCGGCACGTCCGGTCCTTGTTCCGGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCAc                                                                                                                                  <  2:179971/246‑1 (MQ=255)
                                                                                                              aataaACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACgg                                                                                                                                <  1:93973/251‑1 (MQ=255)
                                                                                                                 aaaCGCCCGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGtt                                                                                                                             <  2:384328/251‑1 (MQ=255)
                                                                                                                 aaaCGCCAGGGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGtt                                                                                                                             <  1:430524/251‑1 (MQ=255)
                                                                                                                  atcGCCAGCGCCACGTCCGGTCCTTGTTCCTGCCGCCACTCTGTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATTGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGttt                                                                                                                            <  1:430094/249‑1 (MQ=255)
                                                                                                                    cGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTa                                                                                                                          <  1:425751/251‑1 (MQ=255)
                                                                                                                          cGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTgggg                                                                                                                                                                                                           >  2:148878/1‑165 (MQ=255)
                                                                                                                          cGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTgggg                                                                                                                                                                                                           <  1:148878/165‑1 (MQ=255)
                                                                                                                               cGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTTGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCAt                                                                                                                 >  2:46481/1‑249 (MQ=255)
                                                                                                                                 tCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑GCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCAACGGGTTTTAGATGAACATAGTg                                                                                                             >  2:136429/1‑251 (MQ=255)
                                                                                                                                  ccGGTCATTGTTCATACCGACCCTCTTTCAGTGTTTGCGCAGCGTAAGAGGGGCCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATACTGGCTTCTCGCCAACCAACAGG‑TCAAAAATGTGTTTACGCTGCACGCTCTTCCCCTCCGGCTTGGGACGATGCGG‑GGGAACATTATTAGGGATGATTGTGACCCGCTTCCGACCAATCAATTTCGCCCCCGTTTCCACGGGGTTTATATTACCATANTGc                                                                                                            >  2:404922/1‑251 (MQ=255)
                                                                                                                                                    gccACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTgccg                                                                                          <  1:366322/251‑1 (MQ=255)
                                                                                                                                                                  gtTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTg                                                                            >  1:380076/1‑251 (MQ=255)
                                                                                                                                                                    ttGCGCAGTGTAAGAGGTGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGt                                                                           <  2:61171/250‑1 (MQ=255)
                                                                                                                                                                    ttGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCGTAGTGCACCGTATCAAAGGTTCCGCCAAACATAGCCTGt                                                                           <  1:395797/250‑1 (MQ=255)
                                                                                                                                                                           gtgtAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGGGCCGCCAAACAGAGCCTGTAAAGAttt                                                                   >  1:81479/1‑251 (MQ=255)
                                                                                                                                                                              tAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCat                                                                >  2:113266/1‑251 (MQ=255)
                                                                                                                                                                                    ggggCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGAGGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGGGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAAc                                                          >  1:205513/1‑251 (MQ=255)
                                                                                                                                                                                                       cgcgTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑GCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTcc                                         <  2:343363/249‑1 (MQ=255)
                                                                                                                                                                                                         cgTTCCTCAAGAGTAAATAATGGCTTGTCGGCAATAGCCCGT‑TTAAACATGTGTTTTAGCTGCACGCTGTTCGCTTCCGGCTGGGGACGCCGGGG‑AGGAACATTTTTAGGGATGATTGTGACCCGCGTCAGCCCAATCAACTTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGAGCCGCAAAACAGAGCCTGCAAAGATTTCACATCAACCGTCGCTAAATACGTCCGcc                                      >  1:245694/1‑250 (MQ=255)
                                                                                                                                                                                                          gTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTCAAGATTTCATATCAACCGTCGATACATACGTCCGcc                                      >  2:445613/1‑249 (MQ=255)
                                                                                                                                                                                                                         ataatGGCTTGTCGTCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGATGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAAATGACCATAGTGCACCGGATCAAAGGTGCCGCCCAACAGAGCCTGTAAAGATTTCCTATCAACCGCCGACAAATACGTCCGCCAGGGGTTTCTGGcaca                      >  2:308240/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATa                                                  >  1:345268/1‑212 (MQ=255)
                                                                                                                                                                                                                                   gTCGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATa                                                  <  2:345268/212‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAGCCGTCGATAAATGCGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATa          <  2:86706/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGCAATCGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATa          >  1:112762/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cGGCAATCGACAGG‑TCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGTGACGATGCGG‑AGATACATTATTAGGGATGATTGTGACCCGCGTCAGAACAGTCAAATTCGCCAGCGTTTCCACGGGTTGTAGATGACCATAGTGAACCGGATCAAAGGTACCGCCAAACAGAGAATGTAAAGATTTCATAGCAACCGTCGGGAAATACGTCCGCAAGGGGTTTATGGCACAACAGAAGAGATa          >  2:152048/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cGCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGAACCGCGTCAGACCAATCAAATTCGCCAGCGTTGCCACGGGTTTTAGATGACCATAGTGCACCGCATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCCGGGGTTTATGGCCCACCAGAAGAGATAACCCTTcc  >  2:391352/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCCAGT‑TCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTcc  >  2:301862/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
GCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGT‑TCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGG‑AGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCC  >  NC_000913/670148‑670634

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