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,5150CT95.2% 206.2 / ‑3.5 62M311I (ATG→ATAdctAC4 dicarboxylate/orotate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (37/22);  total (38/24)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGATACAG  >  NC_000913/3682281‑3682761
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
gACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCACCTTTAGAAGAAAGCAGCAACACGATTAACAGCGAGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCAGGGCGATAAACACCGCCGCTTTTGTCAGGgataacc                                                                                                                                                                                                                                         >  2:206897/1‑248 (MQ=255)
  ccAGGTTAGTCAGCCCCGGAGCTTCTGACAAAAAGCGGGAGATCCCGAGGATCGGCGCAACCCCCGCACCCGGCAAAGGGCCCCCCGCAGAGAGCGTCGCCCCCAGCCCGAAAAACCCACTCCCCGTACCCCCTGCCGCCCTTTAAGAAGAAAGCAGCAACACTATAAACAGCGTGATTGGGTGGACGATACCCATCTTACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGAt                                                                                                                                                                                                                                       <  1:514755/251‑1 (MQ=255)
    tggTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCCCCAGCACGATAAAGCCACTACCAGTTACCCCTGCCGCCCCTTTAGAACACAGCAGAACCACGAATAACAGCGTGATTTGGTGAGCGACATCCATCTAACTGTTAGTGGCCTGGGCTATAAACACCGCAGCTATTGTAAGGTATATCggtgg                                                                                                                                                                                                                                     >  2:675464/2‑247 (MQ=255)
         aGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGCGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCTGGTATATCGATGTGCCa                                                                                                                                                                                                                                 >  1:648176/1‑250 (MQ=255)
                           tGACATAAAGCGGTAGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGc                                                                                                                                                                                                              >  1:636807/1‑251 (MQ=255)
                            gaaaaaaggCGGCCGATACCGAGGATAACCGCCAGACCCGCTACCGGAAAATGGCCCCCCGCAGAGGGCGCCGCCGCCAGCACAAAAAAGCCCCTCCCCTTTCCCCCTGCCGCCCCTTTAGAAAAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAg                                                                                                                                                                                                               <  1:170934/241‑1 (MQ=255)
                            gACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAATGTTAAACTAGTAGcc                                                                                                                                                                                                             >  1:221088/1‑251 (MQ=255)
                             aCATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGCTTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGGTAAACACCGCCGCTATAGACAGGTATATCAATGTGCCATCAAGGTTAAACGAGTAgact                                                                                                                                                                                                            >  2:429693/1‑248 (MQ=255)
                               aTAAAGAGGACGATACCGAGGATCAGCGCCAGACCCGCTAACGGCAAATGGCCCAACGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAAAACGATTAACAGCGTGATTTGGTGGACGATATCCATATGACTGTTAGTGGCATGGGAGGTAAACACCGCCGCTCATGTCAGGACTATGTCTGTGCCATCAAGATTAAACGAATAgacttg                                                                                                                                                                                                          >  2:196724/1‑246 (MQ=255)
                                 aaaGCGGTCGATACCGCGGATCAGCGCCCGACCCGATACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTAACCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGCTTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAACGCCGCTATTGTCAGGTATATCGATGTGCCCTCAAGGTAAAACAAGTAGCCTGTCg                                                                                                                                                                                                        >  2:437271/1‑251 (MQ=255)
                                        tCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGAGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGAc                                                                                                                                                                                                 >  1:134896/1‑251 (MQ=255)
                                           atcccggggATCACCGCCAGCCCCGCTCCCGGCAAAGGGCCCACCGCAGAGAGCGTCCCCGCCAGCACGATAAACCCACTCCCCTTTCCCCCTGCCCCCCCTTTAGAAAAAAGCAGCACCACGATTAACAGCGTTATTTGGTGGACAATATCCATCTGATTTTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAg                                                                                                                                                                                              <  1:731071/244‑1 (MQ=255)
                                              ccggggacaaccgcCAGCCCCCCTACCGGAAAATGGCCCCCCGCAGAGAGCGTCGCCGCCAGCACAAAAACCCCCCTCCCCGTTCCCCCGCCCGCCCCTTAAGAAAAAAGCAGCAACACGATTAACAGCGTGATTTGGTGAACGATACCCATCTGCCTGTTAGTGGCCTGGGCGAAAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGccc                                                                                                                                                                                           <  1:782041/240‑1 (MQ=255)
                                              ccGAGGATCAGCGCCAGCCCCCCCCCCGGCAAATGCCCCCCCGCAGAGGGCGTCGCCGCCAGCACGAAAAACCCACTCCCCGTTACCCCGGCCGCCCCTTTAAAAGAAAGCAGCAACACGATTAACAGCGTAATTTGGTGGACGATATCCATCTGACTGTTAGGGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGccc                                                                                                                                                                                           <  1:722314/251‑1 (MQ=255)
                                               cGAGGATAAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCCGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGaccc                                                                                                                                                                                          >  2:537970/1‑251 (MQ=255)
                                                    aTCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCcaccac                                                                                                                                                                                     >  1:356646/1‑251 (MQ=255)
                                                      cAGCGCCAGACCCGCTCCCGGGAAATGGCCCACCGCTGAGAGCGCCCCCGCCAGCTCGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGAGATCCATCGGACTGTTAGTGGCGTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCg                                                                                                                                                                                   <  1:58924/251‑1 (MQ=255)
                                                        cccccagccccgcacccGGAAAATGGCCCCCCGCAGAGAGCGTCCCCGCCAGCACAATAAACCCACCCCCCTTTCCCCCTGCCCCCCTTTAAGAAGAAAGCAGCAACACGATAAACAGGGTGATTTGGTGGACAATACCCATCTGACTTTTAGTGGCCTGGGGGAAAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGAt                                                                                                                                                                                 <  1:255701/236‑1 (MQ=255)
                                                              cccccccaacgggaaaaggcccccccGCAGAGAGCGGCGCCCCCACCACAATAAACCCACTACCCGTTTCCCCTGCCGCCCCTTTAAAAAAAAGCAGCACCACGATTAACAGCGTGTTTTGGTGGACGATATCCATTTGACTGTTAGTGGCCGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACgg                                                                                                                                                                           <  1:631300/227‑1 (MQ=255)
                                                                      aCCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCg                                                                                                                                                                     >  1:50759/1‑249 (MQ=255)
                                                                                 cccccccGCAGAGAGCGCCGCCGCAAGCACGAAAAACCCACTCCCCTGTCCCCCTGCCGCCCTTTAAGAAAAAAGCAGCAACACGATTAACGGCGGGATTTGGTGGACGATACCCATCTGACTTTTAGTGGCCGGGGCGAAAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCt                                                                                                                                                        <  1:57616/246‑1 (MQ=255)
                                                                                  cccACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGCAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATATGCCATCCAGGTTAAACGAGTAGCCTGTAGGGATGACCAGCCCCAACCCCGATTTACGGCCACCGGGTTTCTCCCTCtt                                                                                                                                                       >  2:182681/1‑251 (MQ=255)
                                                                                     cccGCAGAGAGCGCCGCCGCCAGCACGATAAACCCCCTACCCGTCCCCCCTGCCGCCCCTTAAAAAAAAAGAAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTAAGTGGCTGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTc                                                                                                                                                    <  1:828352/250‑1 (MQ=255)
                                                                                      ccGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTaa                                                                                                                                                                                                               >  2:161201/1‑190 (MQ=255)
                                                                                           gagagCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGCTTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGctt                                                                                                                                              >  1:457440/1‑249 (MQ=255)
                                                                                            agagCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCACTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGGGCAt                                                                                                                                              >  1:864895/1‑250 (MQ=255)
                                                                                             gagCGCCGCCGCCAGCACAAAAAACCCACCCCCCGTACCCCCTGCCGCCCCTTAAGAAGAAAGCAGCAACAGGATTAACAGGGTGATTTGGTGGACAATACCCACTTGACTGTTAGTGGCCGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATAc                                                                                                                                            <  1:820579/251‑1 (MQ=255)
                                                                                                  tCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCCTCTTGTCGAGCCTACGCGGc                                                                                                                                       >  2:274456/1‑251 (MQ=255)
                                                                                                        ccggcccaaaaaaGCCACTCCCCGTTACCCCGCCCGCCCTTTAAGAAAAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACTATACCCATCTGACTTTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGcc                                                                                                                                 <  1:13190/241‑1 (MQ=255)
                                                                                                              cGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACTATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCTCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCg                                                                                                                           <  2:610353/251‑1 (MQ=255)
                                                                                                                 tAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACAAGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGACGCTATTGTCAGGTATATCGATGAGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACGAGCCCAACCACCGATTTACGACAGCCGAGTTTCTCCATCTTGTCGAGAATACGCGGCAGCGCCGACTCGGa                                                                                                                         >  1:570315/1‑250 (MQ=255)
                                                                                                                     gCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGAGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGATTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCCGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGACGAtt                                                                                                                    >  1:809789/1‑250 (MQ=255)
                                                                                                                      ccACTACCCGTTACCCCGCCCGCCCCTTAAGAAAAAAGCAGCAACACGATTAACAGCGTGATGTGGGGGACGATATCCATCTGTCTGTTAGTGGCCTGGGCGAAAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGa                                                                                                                   <  1:617970/251‑1 (MQ=255)
                                                                                                                                     cctcccgccCCTTGAGAAAAAAGCAGCAACTCGATTAACAGCTTGTTGTGGTGGACGATATCCATCTGATTTTGAGGGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAAt                                                                                                   <  1:296952/247‑1 (MQ=255)
                                                                                                                                       tGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGCAGTCCCCAGTACAATc                                                                                                  >  1:579655/1‑250 (MQ=255)
                                                                                                                                           gccCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGCTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTcccccgtccactcaac                                                                                                >  2:873333/1‑237 (MQ=255)
                                                                                                                                             cccTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCACGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGAGGAAGCGCCGACTCGGAAGATGAAGTCCCCAATACCATCAGCAGtt                                                                                           >  2:838173/1‑251 (MQ=255)
                                                                                                                                                     aaaaaaaGCAGAAACACGATAAACAGTGTGATGTGGTGGACGATACCCACCGGACTGTAGGTGGCGGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCCCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACg                                                                                    <  1:85844/247‑1 (MQ=255)
                                                                                                                                                           agcagcAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCCGGGCGATAAACACCGCCACTATTGTCAAGTATATCGAAGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCAACGATTTACGGCAGCCGAGTTTCGCCCTCTAGGCGTGCAACCGCGGCCACGCCGACTAGGAAACTGAAGTCCCCCAGACAATcagcacatctgacctgagtgg                                                                              >  2:590338/1‑235 (MQ=255)
                                                                                                                                                            acagcaCCACGATTAACAGCGGGTTTTGGTGAACGATAACAATCTGCCTGTTAGTGGCCGGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAg                                                                             <  1:342814/250‑1 (MQ=255)
                                                                                                                                                             cagcaACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGAGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGAGGAAGTCCCCAGTACAATCAGCCGTTCTTCACGGATGAAGc                                                                            >  2:678703/1‑251 (MQ=255)
                                                                                                                                                                aaccacGATTAACGGCGGGATTTGGTGGAGGATATCCTTCTTACGTTTTGTGGCTTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGa                                                                         <  1:682084/248‑1 (MQ=255)
                                                                                                                                                                    acGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCCCGGATGTAGCGGATaa                                                                      >  1:421148/1‑250 (MQ=255)
                                                                                                                                                                      gATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAAtt                                                                   >  1:60353/1‑251 (MQ=255)
                                                                                                                                                                      gATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAAtt                                                                   >  1:781794/1‑251 (MQ=255)
                                                                                                                                                                       aTTAACAGCTGGTTTTGGTGGACGATATCCTCCTGTTTGTGATGGGCCGGGGCTATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGGTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGTAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAAttt                                                                  <  1:794022/251‑1 (MQ=255)
                                                                                                                                                                          aaCAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGCGCATACGCGGCAGCGCCGGCGCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCAAGGATTGAGGGGATAAATTTGa                                                                >  2:230782/1‑250 (MQ=255)
                                                                                                                                                                            cAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGa                                                                >  1:591043/1‑248 (MQ=255)
                                                                                                                                                                                gTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCAGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACt                                                         >  1:242679/1‑251 (MQ=255)
                                                                                                                                                                                   aTTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGa                                                       >  1:142200/1‑250 (MQ=255)
                                                                                                                                                                                     ttGGTGGACGATAGCATTCTGACTTTAAGGGCCGTGGGTGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGaaa                                                     <  1:126573/250‑1 (MQ=255)
                                                                                                                                                                                        gtggACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCTAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGTTACTGAAAcacg                                                 >  1:353707/1‑248 (MQ=255)
                                                                                                                                                                                              gttcccctttttactgtgggTGGCCTGGGACAAAAACCCCGCCGCTATTGTTAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCGCCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCtt                                           <  2:464784/232‑1 (MQ=255)
                                                                                                                                                                                              gATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCtt                                           >  1:50325/1‑251 (MQ=255)
                                                                                                                                                                                               atatCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTAGAACCAGTCGCttt                                          >  2:463229/1‑251 (MQ=255)
                                                                                                                                                                                                 atCCATCTGACTGTTAGTGGCCTGGGCGATAACCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAACCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGAGCCGACTCGGAAGATGAAGTCACCCGTACCATCCGCAGTTCTTCACGGCAGTAGCGGATAAATTTGAAGATAATGAAACCAATAGCttt                                          >  2:513890/1‑249 (MQ=255)
                                                                                                                                                                                                   ccATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCg                                      >  1:354337/1‑251 (MQ=255)
                                                                                                                                                                                                                      cgtGGGCAATAAACACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGc                    <  1:873333/249‑1 (MQ=255)
                                                                                                                                                                                                                                aaaCACCGCCGCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAATTTCTCCATCTTGTCGAGCCTACGCGGCCGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATaaccacaccaccaaac          >  2:268805/1‑236 (MQ=255)
                                                                                                                                                                                                                                  acacCGCCGCTAGTGGCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCCCAAACAgg        <  1:218269/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccgccgCTATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATTAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGCACCCAATACCCGCCCCCCAAACAGGAt      >  2:73213/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ccgccgCCATTGTCAGGTATATCGCTGTGCCATCAAGGTTAAACTAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAATTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGATa     <  1:452459/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ccgctATTGTCAGGTATATCGATGAGCCATCAAGGTTAAACGAGTAGCCTGCCGGGCTGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCCGCGCGGACTCGGAAGATGCAGCCCCCAGTACAATAAAACGTTCTTCACGGATGTAACGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGGTTGAACCCAATTAAAGCACCACAAACAGGATAcgg  >  2:437187/1‑249 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
GACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGATTTACGGCAGCCGAGTTTCTCCATCTTGTCGAGCATACGCGGCAGCGCCGACTCGGAAGATGAAGTCCCCAGTACAATCAGCAGTTCTTCACGGATGTAGCGGATAAATTTGAAGATACTGAAACCAGTCGCTTTAGCGATTGAACCCAATACCAGCACCACAAACAGGATACAG  >  NC_000913/3682281‑3682761

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