Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,863,801 T→G 10.2% intergenic (+35/‑49) gapA → / → yeaD glyceraldehyde‑3‑phosphate dehydrogenase A/D‑hexose‑6‑phosphate epimerase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,863,8010TG10.2% 303.0 / 25.0 98intergenic (+35/‑49)gapA/yeaDglyceraldehyde‑3‑phosphate dehydrogenase A/D‑hexose‑6‑phosphate epimerase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (46/42);  new base G (6/4);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

TAAAGCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTCGATC  >  NC_000913/1863655‑1863946
                                                                                                                                                  |                                                                                                                                                  
tAAAGCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCg                                                                                                                                                  >  2:404935/1‑149 (MQ=255)
tAAAGCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCg                                                                                                                                                  >  2:119841/1‑149 (MQ=255)
    gCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCg                                                                                                                                              >  4:156618/1‑149 (MQ=255)
          aTCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑ttt                                                                                                                                       <  1:250831/149‑1 (MQ=255)
          aTCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTACTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑ttt                                                                                                                                       <  4:410328/149‑1 (MQ=255)
             gCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑tttttt                                                                                                                                    >  1:825367/1‑149 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑ttttt                                                                                                                                     >  4:382947/1‑139 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑ttttt                                                                                                                                     <  3:382947/139‑1 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGAt                                                                                                                           >  2:642350/1‑149 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGAt                                                                                                                           >  1:405356/1‑149 (MQ=255)
                         aaCTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATaaa                                                                                                                        >  2:278064/1‑149 (MQ=255)
                           cTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAAt                                                                                                                      >  3:118403/1‑149 (MQ=255)
                              catGAACCTGATATCCTGGTACGACAACGAAACCGGCTACTCCAACAATGTTATGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAATAGCGACTTTGGTCGCTC‑TTTTTTTTACCTGATAAAATGaa                                                                                                                   <  2:573228/147‑1 (MQ=255)
                              cGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGaa                                                                                                                   <  4:827741/149‑1 (MQ=255)
                                tGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGt                                                                                                                 <  3:131783/149‑1 (MQ=255)
                                tGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGt                                                                                                                 <  3:38295/149‑1 (MQ=255)
                                    aCTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             <  1:500058/149‑1 (MQ=255)
                                    aCTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             >  1:745747/1‑149 (MQ=255)
                                     cTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCACTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAg                                                                                                            <  2:779911/149‑1 (MQ=255)
                                       ggTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGa                                                                                                          <  2:424905/149‑1 (MQ=255)
                                       ggTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGa                                                                                                          <  1:119841/149‑1 (MQ=255)
                                           tCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      >  2:486345/1‑149 (MQ=255)
                                           tCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      <  1:522092/149‑1 (MQ=255)
                                              tGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACAATATGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTATAGGACTGCGTc                                                                                                   >  4:269679/1‑149 (MQ=255)
                                               ggTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCa                                                                                                  <  1:250425/149‑1 (MQ=255)
                                                gTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCAt                                                                                                 <  2:825367/149‑1 (MQ=255)
                                                  aCGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTAGCTGATAAAATGAAGTTAAAGGACTGCGTCATGa                                                                                               <  3:824387/149‑1 (MQ=255)
                                                    gACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAtt                                                                                             >  2:576988/1‑149 (MQ=255)
                                                    gACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAtt                                                                                             >  3:393688/1‑149 (MQ=255)
                                                     aCAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTa                                                                                            >  2:142240/1‑149 (MQ=255)
                                                      cAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTaa                                                                                           <  2:431967/149‑1 (MQ=255)
                                                       aaCGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCa                                                                                                  <  2:351222/142‑1 (MQ=255)
                                                       aaCGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCa                                                                                                  >  1:351222/1‑142 (MQ=255)
                                                          gAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGaaa                                                                                       >  3:126622/1‑149 (MQ=255)
                                                          gAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGaaa                                                                                       <  1:101150/149‑1 (MQ=255)
                                                               cGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAtttt                                                                                  >  3:650194/1‑149 (MQ=255)
                                                                ggTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAttttt                                                                                 <  1:794164/149‑1 (MQ=255)
                                                                ggTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAttttt                                                                                 <  1:584865/149‑1 (MQ=255)
                                                                ggTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAttttt                                                                                 <  3:126995/149‑1 (MQ=255)
                                                                 gTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTg                                                                                <  1:789419/149‑1 (MQ=255)
                                                                      tCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCtt                                                                           <  2:236266/149‑1 (MQ=255)
                                                                      tCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑ttttta                                                                                                                                    <  3:98154/92‑2 (MQ=37)
                                                                      tCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑ttttta                                                                                                                                    >  4:98154/1‑91 (MQ=37)
                                                                        caacaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTcc                                                                         >  1:134586/1‑149 (MQ=255)
                                                                         aacaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCg                                                                        >  2:322693/1‑149 (MQ=255)
                                                                          acaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTAc                                                                                                                                >  3:243318/1‑92 (MQ=255)
                                                                          acaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTAc                                                                                                                                <  4:243318/92‑1 (MQ=255)
                                                                            aaaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTc                                                                     >  2:546661/1‑149 (MQ=255)
                                                                            aaaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTc                                                                     <  1:142240/149‑1 (MQ=255)
                                                                            aaaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTc                                                                     <  2:269566/149‑1 (MQ=255)
                                                                             aaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCa                                                                    <  2:540531/149‑1 (MQ=255)
                                                                                ttCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCg                                                                 >  4:494432/1‑149 (MQ=255)
                                                                                  cTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGaa                                                               <  1:285976/149‑1 (MQ=255)
                                                                                  cTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGaa                                                               <  2:405356/149‑1 (MQ=255)
                                                                                    ggACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACa                                                             >  3:208895/1‑149 (MQ=255)
                                                                                     gACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACaa                                                            >  3:126074/1‑149 (MQ=255)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             <  3:76770/98‑1 (MQ=255)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAAt                                                          <  3:501973/149‑1 (MQ=255)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑tttttag                                                                                                                                   <  4:823293/76‑3 (MQ=25)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑tttttag                                                                                                                                   >  3:823293/1‑74 (MQ=25)
                                                                                       ccTGATCGATCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             >  4:76770/1‑98 (MQ=255)
                                                                                        cTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGa                                                                                               <  2:751815/111‑1 (MQ=255)
                                                                                        cTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGa                                                                                               >  1:751815/1‑111 (MQ=255)
                                                                                          gATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAAtctc                                                       >  2:74817/1‑149 (MQ=255)
                                                                                          gATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAAtctc                                                       >  3:784513/1‑149 (MQ=255)
                                                                                                tCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGt                                                 >  3:368095/1‑149 (MQ=255)
                                                                                                    aTCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGc                                                                               <  1:118831/115‑1 (MQ=255)
                                                                                                    aTCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGc                                                                               >  2:118831/1‑115 (MQ=255)
                                                                                                      ctcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑ttt                                                                                                                                       <  1:281285/57‑1 (MQ=21)
                                                                                                      ctcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTC‑ttt                                                                                                                                       >  2:281285/1‑57 (MQ=21)
                                                                                                       tcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAt                                                                                              >  4:501652/1‑97 (MQ=255)
                                                                                                       tcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAt                                                                                              <  3:501652/97‑1 (MQ=255)
                                                                                                       tcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTcc                                          >  4:411129/1‑149 (MQ=255)
                                                                                                         cAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCcg                                        >  4:388273/1‑149 (MQ=255)
                                                                                                          aaaTAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCcgt                                       >  4:424213/1‑149 (MQ=255)
                                                                                                             tAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCcgtcgt                                    >  4:824231/1‑149 (MQ=255)
                                                                                                                 ttGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAAc                                >  1:409867/1‑149 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            <  1:141560/149‑1 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            >  1:54706/1‑149 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            <  2:134586/149‑1 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            <  3:221744/149‑1 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑‑TTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGAt                           >  2:634755/1‑149 (MQ=255)
                                                                                                                      aTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGAt                           <  1:404935/149‑1 (MQ=255)
                                                                                                                             tgtgATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      <  2:479172/67‑1 (MQ=255)
                                                                                                                             tgtgATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      >  1:479172/1‑67 (MQ=255)
                                                                                                                              gtgATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGa                   <  1:760509/149‑1 (MQ=255)
                                                                                                                               tgATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGAc                  <  3:304940/149‑1 (MQ=255)
                                                                                                                               agATCTCAAAAGAGCGACTTCGGCCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGAc                  >  2:716444/2‑149 (MQ=255)
                                                                                                                                 taaaaaaaaGAGCGACCGAAGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCt                >  4:330319/5‑149 (MQ=255)
                                                                                                                                 aTCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCt                >  2:66113/1‑149 (MQ=255)
                                                                                                                                     aaaaaGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCAtt            >  3:204265/1‑149 (MQ=255)
                                                                                                                                      aaaaGAGCGACCGAAGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATtg           >  2:780785/1‑149 (MQ=255)
                                                                                                                                       aaaGAGCGACCGAAGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATtgt          >  2:72183/1‑149 (MQ=255)
                                                                                                                                        aaGAGCGACTTCGCACGATC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGTTTAAGAAAATTTTTGCCATTCAGGTCATCACACAAATCTCACCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATtgtg         >  4:715109/1‑149 (MQ=255)
                                                                                                                                         agagCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTgg        >  3:878022/1‑149 (MQ=255)
                                                                                                                                         agagCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTgg        <  1:750917/149‑1 (MQ=255)
                                                                                                                                           agCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTc      <  3:411129/149‑1 (MQ=255)
                                                                                                                                               aCTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTCGatc  >  2:471120/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAAAGCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTC‑TTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTCGATC  >  NC_000913/1863655‑1863946

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