Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,682,515 C→T M311I (ATG→ATA dctA ← C4 dicarboxylate/orotate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,682,5150CT98.2% 185.7 / ‑3.8 56M311I (ATG→ATAdctAC4 dicarboxylate/orotate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (40/15);  total (40/16)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

GTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGATAC  >  NC_000913/3682275‑3682759
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGt                                                                                                                                                                                                                                            >  1:463975/1‑251 (MQ=255)
  tACCGACCAGGTTAGTCAGCGCACGCGCTTCTGACGTAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGATACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGACAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGATAGTGGCCTGGGCGATAACACCCGCCGCTATTCTCAGGta                                                                                                                                                                                                                                           >  2:565879/1‑250 (MQ=255)
     cGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGCTTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACCCCGCCGCTATTGTCAGGGatatt                                                                                                                                                                                                                                       >  2:141597/1‑250 (MQ=255)
         cAGGTTAGTCATCGCACGCGCTTCTGACATAAAGAGGTCGATACCGAGGATCAGCGCCCGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGATACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGAGATAAACACCGCCGCTAATGTCGGGTATATCGAtg                                                                                                                                                                                                                                   >  2:82256/1‑251 (MQ=255)
                gTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATc                                                                                                                                                                                                                            >  1:699491/1‑251 (MQ=255)
                  cAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATACACACCGCCGCTCTTGTCAGGTATATCGATGTGCCATCa                                                                                                                                                                                                                           >  1:27498/1‑250 (MQ=255)
                      gcACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGtt                                                                                                                                                                                                                      >  1:668474/1‑251 (MQ=255)
                                    cATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATACAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGGATATCGATTTGCCATCCAGGTTCAAc                                                                                                                                                                                                                  >  2:204356/1‑241 (MQ=255)
                                       aaaGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCg                                                                                                                                                                                                     >  1:645701/1‑251 (MQ=255)
                                          gCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCAAGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGAATATCGATGTGCCATCAAGGTTAAACGTGTAGCCTGTCggg                                                                                                                                                                                                   >  1:492387/1‑250 (MQ=255)
                                                       aGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCccac                                                                                                                                                                                     >  1:642880/1‑251 (MQ=255)
                                                           ccagcgcaaccccccctcccgaaaaaTGCCCCACCGCAGAGGGCGCCGCCGCCAGCACAAAAAAGCCCCCACCCGTTCCCCCGCCCCCCCCTTTAGAAAAAAGCAGCAACACGATTAACAGCGTGATGTGGTGGACGAAACCCATCTGTCTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCaccacc                                                                                                                                                                                 <  1:227276/228‑1 (MQ=255)
                                                                   aGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACg                                                                                                                                                                         >  1:668328/1‑251 (MQ=255)
                                                                                 cAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTtct                                                                                                                                                           >  1:192328/1‑251 (MQ=255)
                                                                                     tGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCACCCCCACCACCGAATTACGGCAGCCGCGTTTCTCCAt                                                                                                                                                       >  2:570442/1‑251 (MQ=255)
                                                                                       gCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTCCCCGTTCCCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCTATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCt                                                                                                                                                     <  1:237870/251‑1 (MQ=255)
                                                                                         ccACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCCCCGATTTACGGCAGCCGCGTTTCTCCATTTTg                                                                                                                                                   >  2:142465/1‑251 (MQ=255)
                                                                                              gCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGAtt                                                                                                                                                                             >  1:221131/1‑220 (MQ=255)
                                                                                                agagagCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGc                                                                                                                                             >  1:221225/1‑250 (MQ=255)
                                                                                                      cGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGCTGTGCCATCAAGGTTAAACGAGTAGCCTGTCGCGATGACCAACCCCACCACCGATTTACGGCAGCCGGGTTTCTCCATCTTGTCGAGCATAcgcc                                                                                                                                      >  2:331076/1‑250 (MQ=255)
                                                                                                      cGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGAGGACGATATACATCTGACTGTTAGAGGCCTGGGGGATAAAAACCGCAGCTATTGTAAGGCATATAGATGTGCCATAAAGGATAAAGGAGGAGAATGTAGGGATGAACAGCCACAACACCGAATTACGGCAGCAGAGTTTATCCAtatagaagagcataag                                                                                                                                        >  1:459142/1‑234 (MQ=255)
                                                                                                        ccgccgcAAGAACGAAAAAGCCACTCCCCGTTACCCCGGCCGCCCCTTAAGAAGAAAGCAGAAACACGATTAACAGCGTGATTGGGGGGACGATATCCATCTGACTGTTAGTGGCCGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGc                                                                                                                                    <  1:279267/250‑1 (MQ=255)
                                                                                                         cgccgcCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGc                                                                                                                                    >  1:7264/1‑250 (MQ=255)
                                                                                                                        caaGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGGCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGGGTTTCTCCATCTTGTCGAGCATACGCGGCaagatcggcagagagg                                                                                                                    >  2:19509/2‑236 (MQ=255)
                                                                                                                        aaacccccgaCCCGTTCCCCCGGCCGCCCCTTAAGAAGAAAGCAGCAACACGATAAACAGGGTGATTTGGTGGACGATATCCATCTGACTGTTAGGGGCCTGGGCAAAAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCa                                                                                                                                   <  1:19509/227‑1 (MQ=255)
                                                                                                                        aaaaccactcccggtccccccGGCCGCCCTTTTAGAAAAAAGCAGCAACACGATTAACAGCGTGATTGGGGGGACGATATCCATTTGATTGTTGGGGGCTTGGGCGAAAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAg                                                                                                                    <  1:119104/234‑1 (MQ=255)
                                                                                                                        aaaGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATg                                                                                                                                                                                                >  2:374485/1‑175 (MQ=255)
                                                                                                                                       ccccccGGCCGCCCCTTAAGAAAAAGGCAGCAACCCGATTAACAGCGTGTTTTGGTGGACGATATCCTTCTGACTGTTAGTGTCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGt                                                                                                      <  1:141597/248‑1 (MQ=255)
                                                                                                                                        cccccGGCCCCCCCTTAAGAAAAAAGCACAAACACGATTACCAGGGGAATTTGGTGAACGATACCCATCTGACTGTTAGTGGCCTGGGCGAAAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTAc                                                                                                    <  1:142465/250‑1 (MQ=255)
                                                                                                                                              gccgccCCTTTAAAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACGGTTAGTGGCCGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATc                                                                                                <  1:331076/249‑1 (MQ=255)
                                                                                                                                              cccgccCTTTTAGAAAAAAGCAGCAACTCTATTAACAGCGTGATTTGGTGGACGATACCAATCTGCCTGTGAGGGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCTGAAGATGAAGTCCCCAGTACAATcag                                                                                              <  1:557554/250‑1 (MQ=255)
                                                                                                                                                 gccCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATcagc                                                                                             >  1:2131/1‑249 (MQ=255)
                                                                                                                                                               aaGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGCGCATACGCGGCAGCGCCGCATCGGAAGATGAAGTCCCCAGTACAATCAGCAGttcttccc                                                                                   >  2:30792/1‑243 (MQ=255)
                                                                                                                                                               aaGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGGTGAAGTCCCCAGTACAATCAGCAGTTCTTCACggctgg                                                                             >  2:209522/1‑247 (MQ=255)
                                                                                                                                                               aaGCAGCAACACGATTAACAGCGTGATTGGGGGGACGATACCCATCTGACTGTTAGGGGCCGGGGAGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCAc                                                                                   <  1:30792/245‑1 (MQ=255)
                                                                                                                                                               aaGCAGCAACACGATAAACGGCGGGATTTGGGGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGAAAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGt                                                                             <  1:14587/251‑1 (MQ=255)
                                                                                                                                                                      aaCACGATTACCAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGCTGTGCCATCACGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAACTGAAGTCCCCCGTACAATCAGCAGCTCTTCACGGATGAAGCggg                                                                       >  2:449536/1‑249 (MQ=255)
                                                                                                                                                                       ccgcgttaaaaggcttgttttGGTGAACGTTACCCATCTGGTTGTTAGGGGCTGGGGCGAAAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATa                                                                     <  1:427864/233‑1 (MQ=255)
                                                                                                                                                                           gATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAAtt                                                                 >  1:620966/1‑251 (MQ=255)
                                                                                                                                                                           gATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACAACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCCGCGCCGCCTCGGAAGAAGAAGTCCCCAGTCCACTCAGCAGTTCTTCACGGATGTAGCGGCTAAAtt                                                                 >  2:81273/1‑251 (MQ=255)
                                                                                                                                                                             ttAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGt                                                                                                                                                                                                              >  2:433627/1‑108 (MQ=255)
                                                                                                                                                                                  aGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCCGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTCGCGGATAAATTTGAAGa                                                           >  2:165268/1‑250 (MQ=255)
                                                                                                                                                                                  aGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTa                                                                            >  1:175452/1‑233 (MQ=255)
                                                                                                                                                                                   ccttccgttctgggACAATATCCATCTGATTGTTGGTGCCCGGGGAGAAAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCAACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATa                                                         <  1:319552/239‑1 (MQ=255)
                                                                                                                                                                                         tttGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGaa                                                    >  1:429131/1‑250 (MQ=255)
                                                                                                                                                                                              tggACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCAACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATa                                                         >  2:319552/1‑240 (MQ=255)
                                                                                                                                                                                                  cGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGGGCCATCACGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGCCGAGCCTACGCGGCAGCGCCGACTCGGAAGATGCAGTCCCCAGTACAATCAGCAGTTATTCACGGCTGTAGCGGATAAATTTGAAGCTACTGAAACCCGTCGCt                                          >  1:474660/1‑251 (MQ=255)
                                                                                                                                                                                                      cccctcctagttGTTGGTGGCCTGGGGTAAAACCACCGCCGCTATTTGCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAg                                      <  2:18712/240‑1 (MQ=255)
                                                                                                                                                                                                          aTCTGCCTGTTAGTGGCCGGGGCGATAAACAACGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCAACACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGAAGTAGAGGATAATTTTGAAGATACTGAAACCAGTCGCTTTAGCGAt                                  >  1:153519/1‑251 (MQ=255)
                                                                                                                                                                                                                  gTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAAccc                           >  1:451165/1‑250 (MQ=255)
                                                                                                                                                                                                                   gtggtgctgtgcgCGATAAAGTCCGCCGCTGTTGTCAGGGATATCGATGTGCCATCAAGGGTAAACGAGTAGCCTTTCGGTATGACCAGCCTCCCCACCGATTTACGGCAGCAGAGTTTCGCCATCTGGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCTCGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCaa                         <  2:155934/239‑1 (MQ=255)
                                                                                                                                                                                                                         gCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCCAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTAAAATCAGCAGTTCTTCAAGGGTGTAGCGGATAAATTTTAAGATAATGAAACCAGTCGCTTTAGCGATTGAAACCCATAcc                     >  2:467809/1‑249 (MQ=255)
                                                                                                                                                                                                                                cgaaaacctccgccgCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCGCCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACcaca            <  1:82256/242‑1 (MQ=255)
                                                                                                                                                                                                                                 gatcaactccgccgCTATTGTCAGGTATATCGATGTGACATCAAGGTTAACCGAGTAGCCTATCGGGATGACCAGCCCCACCACCGATTTACGACAGCCGAGTTTCTCCATCTTGTCGAGCATCCGCGCCAGCGCCGACTCTGAAGATGAAGTCCCCCGTACAATCAGCAGTTCTTCACGGCAGTAGCGGATAAATTTGCAGATACTTAACACAGTCGCTTTAGCGATTGCACCCAATACCAGCACcacac           >  2:540441/9‑250 (MQ=255)
                                                                                                                                                                                                                                 gATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACcaca            >  1:42613/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cacagccgctATTGTCAGGTATTACGAGTTGCACTCAAGGGTAAACGAGTAGCCTGTCGGGATGACCAGTCACACCACCGAGGTACGGCAGCCGAGTTTCTACATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTGCAATCAGCAGTTCGTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTATCGATTGAACCCAATACCAGCTCCACAAACAgg      <  2:192855/246‑1 (MQ=255)
                                                                                                                                                                                                                                        acCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGAt    >  1:347587/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ccgctATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGATAc  >  1:496719/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGATAC  >  NC_000913/3682275‑3682759

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