Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,861,428 G→T 5.7% intergenic (‑96/+274) ydjL ← / ← yeaC putative Zn‑dependent NAD(P)‑binding oxidoreductase/DUF1315 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,861,4280GT 5.7% 313.9 / 8.4 88intergenic (‑96/+274)ydjL/yeaCputative Zn‑dependent NAD(P)‑binding oxidoreductase/DUF1315 family protein
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base T (2/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-02

TTGTAGCCGCCAAAGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGGAGAAGTC  >  NC_000913/1861289‑1861570
                                                                                                                                           |                                                                                                                                              
ttGTAGCCGCCAAAGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACGCAGATCATTATCTgtgt                                                                                                                                        <  2:1024122/148‑1 (MQ=255)
          cAAAGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTc                                                                                                                              >  2:645334/1‑148 (MQ=255)
            aaGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTc                                                                                                                              <  4:717991/146‑1 (MQ=255)
            aaGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTc                                                                                                                              >  3:717991/1‑146 (MQ=255)
             aGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGa                                                                                                                           <  1:645334/148‑1 (MQ=255)
                   tGCCAAACCGTGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGtt                                                                                                                                       <  2:779340/130‑1 (MQ=255)
                   tGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGa                                                                                                                                               >  1:905711/1‑122 (MQ=255)
                   tGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGtt                                                                                                                                       >  1:779340/1‑130 (MQ=255)
                   tGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCGGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGa                                                                                                                                               <  2:905711/122‑1 (MQ=255)
                      taaaCCGAGCTAGTCCTTTCGTTGCGTCGTCCTATGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATTTTTTTTGAGGTCAACGCTGTATCTGTAAACCACAGATAATGACGTGCGTTTTACAACTCAGATCAAAATTa                                                                                                                  <  2:648925/147‑1 (MQ=255)
                      cAACCCGAGCCAGTGCCTTCATTGCGTCATCCTCTTTGTTTCTTTTTGCTTAGATGGCCATCATACTTACCCACATTATTTTTGTTGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTa                                                                                                                  <  4:685851/148‑1 (MQ=255)
                       aaaCCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATGATGATCTGCGTTTTACAACTCAGATCACAATTAc                                                                                                                 >  3:685851/1‑148 (MQ=255)
                         aCCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTTGTTTTGCTTAGATCGCCATCATACTTACCCACATGATTTTTGATGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTatga                                                                                                               >  1:648925/1‑145 (MQ=255)
                                   tGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATg                                                                                                                                                >  4:925999/1‑105 (MQ=255)
                                   tGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATg                                                                                                                                                <  3:925999/105‑1 (MQ=255)
                                     cTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGtttt                                                                                                                                     >  2:766257/1‑114 (MQ=255)
                                     cTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGtttt                                                                                                                                     <  1:766257/114‑1 (MQ=255)
                                       ttCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTggg                                                                                                 >  3:373374/1‑148 (MQ=255)
                                        tCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATc                                                                                                                                             <  4:1002560/103‑1 (MQ=255)
                                        tCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATc                                                                                                                                             >  3:1002560/1‑103 (MQ=255)
                                           ttGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGc                                                                                                                                          <  1:438045/103‑1 (MQ=255)
                                           ttGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGc                                                                                                                                          >  2:438045/1‑103 (MQ=255)
                                           ttGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTg                                                                                             <  4:373374/148‑1 (MQ=255)
                                                   tCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCg                                                                                                                                         >  3:221914/1‑96 (MQ=255)
                                                   tCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCg                                                                                                                                         <  4:221914/96‑1 (MQ=255)
                                                             ttAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCa                                                                                                                             <  2:204602/98‑1 (MQ=255)
                                                             ttAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGACAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCa                                                                                                                             >  1:204602/1‑98 (MQ=255)
                                                                  tttGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGt                                                                                                                                        <  2:385376/82‑1 (MQ=255)
                                                                  tttGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGt                                                                                                                                        >  1:385376/1‑82 (MQ=255)
                                                                  tttGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCgggg                                                                      >  1:24617/1‑148 (MQ=255)
                                                                       ttAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATcaca                                                                                                                      <  4:798269/95‑1 (MQ=255)
                                                                       ttAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATcaca                                                                                                                      >  3:798269/1‑95 (MQ=255)
                                                                        tAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATg                                                                                                                                                >  4:1139128/1‑68 (MQ=255)
                                                                        tAGATGGCCATCATACCTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATg                                                                                                                                                <  3:1139128/68‑1 (MQ=255)
                                                                         aGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCaaaa                                                               <  2:24617/148‑1 (MQ=255)
                                                                           aTGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGtt                                                                                           <  3:139600/118‑1 (MQ=255)
                                                                           aTGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGtt                                                                                           >  4:139600/1‑118 (MQ=255)
                                                                                  tcatACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAAt                                                      >  3:848536/1‑148 (MQ=255)
                                                                                        ttACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGAt                                                                                 <  1:225965/115‑1 (MQ=255)
                                                                                        ttACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGAt                                                                                 >  2:225965/1‑115 (MQ=255)
                                                                                         tACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGa                                                                                                                                               <  2:564369/52‑1 (MQ=255)
                                                                                         tACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGa                                                                                                                                               >  1:564369/1‑52 (MQ=255)
                                                                                         tACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGt                                                                                     >  1:743856/1‑110 (MQ=255)
                                                                                         tACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGt                                                                                     <  2:743856/110‑1 (MQ=255)
                                                                                            ccACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCg                                                   >  1:821713/1‑141 (MQ=255)
                                                                                            ccACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCg                                                   <  2:821713/141‑1 (MQ=255)
                                                                                               caTGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAGTGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTa                                         >  3:265449/1‑148 (MQ=255)
                                                                                                   aTTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTc                                                                                               <  3:349533/90‑1 (MQ=255)
                                                                                                   aTTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTc                                                                                               >  4:349533/1‑90 (MQ=255)
                                                                                                      tttGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCgg                                                                        <  1:26465/110‑1 (MQ=255)
                                                                                                      tttGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCgg                                                                        >  2:26465/1‑110 (MQ=255)
                                                                                                          aGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATa                              >  4:761232/1‑148 (MQ=255)
                                                                                                          aGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATa                              >  2:1099537/1‑148 (MQ=255)
                                                                                                             tCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGa                                                                                  >  2:81513/1‑93 (MQ=255)
                                                                                                             tCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGa                                                                                  <  1:81513/93‑1 (MQ=255)
                                                                                                                aCATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATa                              >  3:409578/1‑142 (MQ=255)
                                                                                                                aCATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATa                              <  4:409578/142‑1 (MQ=255)
                                                                                                                 cATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCaaaa                                                               <  3:1126563/108‑1 (MQ=255)
                                                                                                                 cATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCaaaa                                                               >  4:1126563/1‑108 (MQ=255)
                                                                                                                   ttAAATCTGTAAAACATAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATgg                                                                                                      <  2:964034/67‑1 (MQ=255)
                                                                                                                   ttAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGTAGAAAATgg                                                                                                      >  1:964034/1‑67 (MQ=255)
                                                                                                                   ttAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTgaga                     >  4:767098/1‑148 (MQ=255)
                                                                                                                      aaTCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTg                                                                                                   >  1:310701/1‑67 (MQ=255)
                                                                                                                      aaTCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTg                                                                                                   <  2:310701/67‑1 (MQ=255)
                                                                                                                        tCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAGTGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTaa                <  4:265449/148‑1 (MQ=255)
                                                                                                                        tCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTaa                >  4:275119/1‑148 (MQ=255)
                                                                                                                         cTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCaaa                                >  4:256907/1‑131 (MQ=255)
                                                                                                                         cTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCaaa                                <  3:256907/131‑1 (MQ=255)
                                                                                                                         cTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAg               <  3:767098/148‑1 (MQ=255)
                                                                                                                           gTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTa                 <  3:587663/144‑1 (MQ=255)
                                                                                                                           gTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTa                 >  4:587663/1‑144 (MQ=255)
                                                                                                                                acgcaGATCATTATCTGTGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCgggg                                                                      >  4:947293/4‑86 (MQ=38)
                                                                                                                                acgcaGATCATTATCTGTGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCgggg                                                                      <  3:947293/83‑1 (MQ=38)
                                                                                                                                 cgcaGATCATTATCTGTGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAgg         <  4:2237/144‑1 (MQ=255)
                                                                                                                                 cgcaGATCATTATCTGTGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAgg         >  3:2237/3‑146 (MQ=255)
                                                                                                                                   caGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGGAGaa     <  3:275119/148‑1 (MQ=255)
                                                                                                                                     gATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGGAGAAGt   >  1:153490/1‑148 (MQ=255)
                                                                                                                                      aTAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTgag                      >  3:536623/1‑128 (MQ=255)
                                                                                                                                      aTAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTgag                      >  1:34433/1‑128 (MQ=255)
                                                                                                                                      aTAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTgag                      <  2:34433/128‑1 (MQ=255)
                                                                                                                                      aTAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTgag                      <  4:536623/128‑1 (MQ=255)
                                                                                                                                      aTAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGGAGAAGTc  >  3:638541/1‑148 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCaaaa                                                               <  3:870854/86‑1 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCaaaa                                                               >  4:870854/1‑86 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGagag                                               <  4:389950/102‑1 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGagag                                               >  3:389950/1‑102 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTg                        >  4:668371/1‑125 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTg                        <  3:668371/125‑1 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGgaga      >  1:198711/1‑143 (MQ=255)
                                                                                                                                       tAATGATCTGCGTTTTACAACGCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGgaga      <  2:198711/143‑1 (MQ=255)
                                                                                                                                        aaTGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTAcc                                       >  2:507108/1‑109 (MQ=255)
                                                                                                                                        aaTGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTAcc                                       <  1:507108/109‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
TTGTAGCCGCCAAAGGCCTTGCCAAACCGAGCCAGTGCTTTCATTGCGTCATCCTCTGTGTTTAGTTTTGCTTAGATGGCCATCATACTTACCCACATGATTTTTGAGGTCAACATTAAATCTGTAAAACACAGATAATGATCTGCGTTTTACAACTCAGATCACAATTACGCAGAAAATGGCTGGGTCTGTTACAGGTTGATGGAAGGCGGGGCGCAAAAAGAGCAAATTCGAGAGGGATTACCGGCCAAATAGGTCTGAGACTTAAGTTCAGGAGAAGTC  >  NC_000913/1861289‑1861570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: