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

GCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCC  >  NC_000913/670148‑670634
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
gcCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCTTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCGTACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACAt                                                                                                                                                                                                                                               <  2:139888/251‑1 (MQ=255)
gcCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACAt                                                                                                                                                                                                                                               >  1:179145/1‑251 (MQ=255)
                     ggACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTg                                                                                                                                                                                                                            <  2:29502/249‑1 (MQ=255)
                     ggACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGtt                                                                                                                                                                                                                          >  1:86706/1‑251 (MQ=255)
                       aCGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCg                                                                                                                                                                                                                        >  2:93973/1‑251 (MQ=255)
                             aCAGACGATCAAATGTGCATTGTCGAGTATCTTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAACCGCCAGCGGCACGTCCGGTCCTTGTTCCGGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCg                                                                                                                                                                                                                  <  2:179145/251‑1 (MQ=255)
                                   gATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACAAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGTACGTCAGGTCCTTGTTCCTGCCGCCACTCTTTCCGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAATCTGCGTTTACGCTGCACGCTGTTCGCTTCCGGCTggg                                                                                                                                                                                                            >  2:71543/1‑251 (MQ=255)
                                    aTCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTgggg                                                                                                                                                                                                           <  1:136921/251‑1 (MQ=255)
                                           gtgCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACAAATAATAAACGCCAGCGTCACGTCCGGTCTTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGTATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATg                                                                                                                                                                                                     <  1:220673/250‑1 (MQ=255)
                                                 tGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCggaggg                                                                                                                                                                                             >  1:29265/1‑250 (MQ=255)
                                                            gTTTCGTATTCGTACCAGGTCGAAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTGGGCAATAGACGTGTAAAAAATGGGTTTACGCGGCACGGTGTTCGCGTCCGGCTGGGGACGCGCGGGAGGAACATAATTAgg                                                                                                                                                                                  >  2:389405/1‑251 (MQ=255)
                                                                         aCCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCGGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGAc                                                                                                                                                                      <  1:326927/250‑1 (MQ=255)
                                                                         aCCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGAc                                                                                                                                                                      <  1:23707/250‑1 (MQ=255)
                                                                            aGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCCGTTCAAACATGTGTTTACGCTGCCCGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCcgc                                                                                                                                                                  >  2:336780/1‑251 (MQ=255)
                                                                                    aaGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGAcc                                                                                                                                                          <  1:450106/251‑1 (MQ=255)
                                                                                       gTCAGCAGTGAACCCTGACCAATAATAAACGCCAGCGGCACGTCGGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAAt                                                                                                                                                       <  2:129866/251‑1 (MQ=255)
                                                                                         cagcagTGAATCCGGACCAATAATAAACGCCAGCGGCACGTCTGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCa                                                                                                                                                     <  2:400603/251‑1 (MQ=255)
                                                                                                  aTCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCa                                                                                                                                            <  2:93159/251‑1 (MQ=255)
                                                                                                           aataatAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTcc                                                                                                                                    >  2:451651/1‑250 (MQ=255)
                                                                                                            ataataAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCAc                                                                                                                                  <  2:92879/251‑1 (MQ=255)
                                                                                                            agaaaaaaCGCCAGCGGCACGTCCGGTCCTTGTTCCGGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCAc                                                                                                                                  <  2:179971/246‑1 (MQ=255)
                                                                                                              aataaACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACgg                                                                                                                                <  1:93973/251‑1 (MQ=255)
                                                                                                                 aaaCGCCCGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGtt                                                                                                                             <  2:384328/251‑1 (MQ=255)
                                                                                                                 aaaCGCCAGGGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGtt                                                                                                                             <  1:430524/251‑1 (MQ=255)
                                                                                                                  atcGCCAGCGCCACGTCCGGTCCTTGTTCCTGCCGCCACTCTGTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATTGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGttt                                                                                                                            <  1:430094/249‑1 (MQ=255)
                                                                                                                    cGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTa                                                                                                                          <  1:425751/251‑1 (MQ=255)
                                                                                                                          cGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTgggg                                                                                                                                                                                                           >  2:148878/1‑165 (MQ=255)
                                                                                                                          cGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTgggg                                                                                                                                                                                                           <  1:148878/165‑1 (MQ=255)
                                                                                                                               cGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTTGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCAt                                                                                                                 >  2:46481/1‑249 (MQ=255)
                                                                                                                                 tCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTGCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCAACGGGTTTTAGATGAACATAGTg                                                                                                             >  2:136429/1‑251 (MQ=255)
                                                                                                                                  ccGGTCATTGTTCATACCGACCCTCTTTCAGTGTTTGCGCAGCGTAAGAGGGGCCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATACTGGCTTCTCGCCAACCAACAGGTCAAAAATGTGTTTACGCTGCACGCTCTTCCCCTCCGGCTTGGGACGATGCGGGGGAACATTATTAGGGATGATTGTGACCCGCTTCCGACCAATCAATTTCGCCCCCGTTTCCACGGGGTTTATATTACCATANTGc                                                                                                            >  2:404922/1‑251 (MQ=255)
                                                                                                                                                    gccACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTgccg                                                                                          <  1:366322/251‑1 (MQ=255)
                                                                                                                                                                  gtTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTg                                                                            >  1:380076/1‑251 (MQ=255)
                                                                                                                                                                    ttGCGCAGTGTAAGAGGTGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGt                                                                           <  2:61171/250‑1 (MQ=255)
                                                                                                                                                                    ttGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCGTAGTGCACCGTATCAAAGGTTCCGCCAAACATAGCCTGt                                                                           <  1:395797/250‑1 (MQ=255)
                                                                                                                                                                           gtgtAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGGGCCGCCAAACAGAGCCTGTAAAGAttt                                                                   >  1:81479/1‑251 (MQ=255)
                                                                                                                                                                              tAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCat                                                                >  2:113266/1‑251 (MQ=255)
                                                                                                                                                                                    ggggCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGAGGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGGGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAAc                                                          >  1:205513/1‑251 (MQ=255)
                                                                                                                                                                                                       cgcgTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTGCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTcc                                         <  2:343363/249‑1 (MQ=255)
                                                                                                                                                                                                         cgTTCCTCAAGAGTAAATAATGGCTTGTCGGCAATAGCCCGTTTAAACATGTGTTTTAGCTGCACGCTGTTCGCTTCCGGCTGGGGACGCCGGGGAGGAACATTTTTAGGGATGATTGTGACCCGCGTCAGCCCAATCAACTTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGAGCCGCAAAACAGAGCCTGCAAAGATTTCACATCAACCGTCGCTAAATACGTCCGcc                                      >  1:245694/1‑250 (MQ=255)
                                                                                                                                                                                                          gTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTCAAGATTTCATATCAACCGTCGATACATACGTCCGcc                                      >  2:445613/1‑249 (MQ=255)
                                                                                                                                                                                                                         ataatGGCTTGTCGTCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGATGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAAATGACCATAGTGCACCGGATCAAAGGTGCCGCCCAACAGAGCCTGTAAAGATTTCCTATCAACCGCCGACAAATACGTCCGCCAGGGGTTTCTGGcaca                      >  2:308240/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATa                                                  >  1:345268/1‑212 (MQ=255)
                                                                                                                                                                                                                                   gTCGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATa                                                  <  2:345268/212‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAGCCGTCGATAAATGCGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATa          <  2:86706/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGCAATCGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATa          >  1:112762/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cGGCAATCGACAGGTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGTGACGATGCGGAGATACATTATTAGGGATGATTGTGACCCGCGTCAGAACAGTCAAATTCGCCAGCGTTTCCACGGGTTGTAGATGACCATAGTGAACCGGATCAAAGGTACCGCCAAACAGAGAATGTAAAGATTTCATAGCAACCGTCGGGAAATACGTCCGCAAGGGGTTTATGGCACAACAGAAGAGATa          >  2:152048/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cGCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGAACCGCGTCAGACCAATCAAATTCGCCAGCGTTGCCACGGGTTTTAGATGACCATAGTGCACCGCATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCCGGGGTTTATGGCCCACCAGAAGAGATAACCCTTcc  >  2:391352/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCCAGTTCAAACATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTcc  >  2:301862/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
GCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCC  >  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: