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

TAAAGCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTCGATC  >  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)
          aTCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCttt                                                                                                                                       <  1:250831/149‑1 (MQ=255)
          aTCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTACTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCttt                                                                                                                                       <  4:410328/149‑1 (MQ=255)
             gCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCtttttt                                                                                                                                    >  1:825367/1‑149 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCttttt                                                                                                                                     >  4:382947/1‑139 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCttttt                                                                                                                                     <  3:382947/139‑1 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGAt                                                                                                                           >  2:642350/1‑149 (MQ=255)
                      gACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGAt                                                                                                                           >  1:405356/1‑149 (MQ=255)
                         aaCTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATaaa                                                                                                                        >  2:278064/1‑149 (MQ=255)
                           cTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAAt                                                                                                                      >  3:118403/1‑149 (MQ=255)
                              catGAACCTGATATCCTGGTACGACAACGAAACCGGCTACTCCAACAATGTTATGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAATAGCGACTTTGGTCGCTCTTTTTTTTACCTGATAAAATGaa                                                                                                                   <  2:573228/147‑1 (MQ=255)
                              cGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGaa                                                                                                                   <  4:827741/149‑1 (MQ=255)
                                tGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGt                                                                                                                 <  3:131783/149‑1 (MQ=255)
                                tGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGt                                                                                                                 <  3:38295/149‑1 (MQ=255)
                                    aCTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             <  1:500058/149‑1 (MQ=255)
                                    aCTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             >  1:745747/1‑149 (MQ=255)
                                     cTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCACTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAg                                                                                                            <  2:779911/149‑1 (MQ=255)
                                       ggTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGa                                                                                                          <  2:424905/149‑1 (MQ=255)
                                       ggTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGa                                                                                                          <  1:119841/149‑1 (MQ=255)
                                           tCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      >  2:486345/1‑149 (MQ=255)
                                           tCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      <  1:522092/149‑1 (MQ=255)
                                              tGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACAATATGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTATAGGACTGCGTc                                                                                                   >  4:269679/1‑149 (MQ=255)
                                               ggTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCa                                                                                                  <  1:250425/149‑1 (MQ=255)
                                                gTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCAt                                                                                                 <  2:825367/149‑1 (MQ=255)
                                                  aCGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTAGCTGATAAAATGAAGTTAAAGGACTGCGTCATGa                                                                                               <  3:824387/149‑1 (MQ=255)
                                                    gACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAtt                                                                                             >  2:576988/1‑149 (MQ=255)
                                                    gACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAtt                                                                                             >  3:393688/1‑149 (MQ=255)
                                                     aCAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTa                                                                                            >  2:142240/1‑149 (MQ=255)
                                                      cAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTaa                                                                                           <  2:431967/149‑1 (MQ=255)
                                                       aaCGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCa                                                                                                  <  2:351222/142‑1 (MQ=255)
                                                       aaCGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCa                                                                                                  >  1:351222/1‑142 (MQ=255)
                                                          gAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGaaa                                                                                       >  3:126622/1‑149 (MQ=255)
                                                          gAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGaaa                                                                                       <  1:101150/149‑1 (MQ=255)
                                                               cGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAtttt                                                                                  >  3:650194/1‑149 (MQ=255)
                                                                ggTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAttttt                                                                                 <  1:794164/149‑1 (MQ=255)
                                                                ggTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAttttt                                                                                 <  1:584865/149‑1 (MQ=255)
                                                                ggTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAAttttt                                                                                 <  3:126995/149‑1 (MQ=255)
                                                                 gTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTg                                                                                <  1:789419/149‑1 (MQ=255)
                                                                      tCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCtt                                                                           <  2:236266/149‑1 (MQ=255)
                                                                      tCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCttttta                                                                                                                                    <  3:98154/92‑2 (MQ=37)
                                                                      tCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCttttta                                                                                                                                    >  4:98154/1‑91 (MQ=37)
                                                                        caacaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTcc                                                                         >  1:134586/1‑149 (MQ=255)
                                                                         aacaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCg                                                                        >  2:322693/1‑149 (MQ=255)
                                                                          acaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTAc                                                                                                                                >  3:243318/1‑92 (MQ=255)
                                                                          acaaAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTAc                                                                                                                                <  4:243318/92‑1 (MQ=255)
                                                                            aaaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTc                                                                     >  2:546661/1‑149 (MQ=255)
                                                                            aaaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTc                                                                     <  1:142240/149‑1 (MQ=255)
                                                                            aaaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTc                                                                     <  2:269566/149‑1 (MQ=255)
                                                                             aaGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCa                                                                    <  2:540531/149‑1 (MQ=255)
                                                                                ttCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCg                                                                 >  4:494432/1‑149 (MQ=255)
                                                                                  cTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGaa                                                               <  1:285976/149‑1 (MQ=255)
                                                                                  cTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGaa                                                               <  2:405356/149‑1 (MQ=255)
                                                                                    ggACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACa                                                             >  3:208895/1‑149 (MQ=255)
                                                                                     gACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACaa                                                            >  3:126074/1‑149 (MQ=255)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             <  3:76770/98‑1 (MQ=255)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAAt                                                          <  3:501973/149‑1 (MQ=255)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCtttttag                                                                                                                                   <  4:823293/76‑3 (MQ=25)
                                                                                       ccTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCtttttag                                                                                                                                   >  3:823293/1‑74 (MQ=25)
                                                                                       ccTGATCGATCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTaaa                                                                                                             >  4:76770/1‑98 (MQ=255)
                                                                                        cTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGa                                                                                               <  2:751815/111‑1 (MQ=255)
                                                                                        cTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGa                                                                                               >  1:751815/1‑111 (MQ=255)
                                                                                          gATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAAtctc                                                       >  2:74817/1‑149 (MQ=255)
                                                                                          gATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAAtctc                                                       >  3:784513/1‑149 (MQ=255)
                                                                                                tCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGt                                                 >  3:368095/1‑149 (MQ=255)
                                                                                                    aTCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGc                                                                               <  1:118831/115‑1 (MQ=255)
                                                                                                    aTCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGc                                                                               >  2:118831/1‑115 (MQ=255)
                                                                                                      ctcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCttt                                                                                                                                       <  1:281285/57‑1 (MQ=21)
                                                                                                      ctcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACCGAAGTCGCTCttt                                                                                                                                       >  2:281285/1‑57 (MQ=21)
                                                                                                       tcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAt                                                                                              >  4:501652/1‑97 (MQ=255)
                                                                                                       tcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGAt                                                                                              <  3:501652/97‑1 (MQ=255)
                                                                                                       tcCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTcc                                          >  4:411129/1‑149 (MQ=255)
                                                                                                         cAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCcg                                        >  4:388273/1‑149 (MQ=255)
                                                                                                          aaaTAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCcgt                                       >  4:424213/1‑149 (MQ=255)
                                                                                                             tAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCcgtcgt                                    >  4:824231/1‑149 (MQ=255)
                                                                                                                 ttGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAAc                                >  1:409867/1‑149 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            <  1:141560/149‑1 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            >  1:54706/1‑149 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            <  2:134586/149‑1 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGa                            <  3:221744/149‑1 (MQ=255)
                                                                                                                     gaTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGAt                           >  2:634755/1‑149 (MQ=255)
                                                                                                                      aTGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGAt                           <  1:404935/149‑1 (MQ=255)
                                                                                                                             tgtgATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      <  2:479172/67‑1 (MQ=255)
                                                                                                                             tgtgATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGc                                                                                                      >  1:479172/1‑67 (MQ=255)
                                                                                                                              gtgATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGa                   <  1:760509/149‑1 (MQ=255)
                                                                                                                               tgATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGAc                  <  3:304940/149‑1 (MQ=255)
                                                                                                                               agATCTCAAAAGAGCGACTTCGGCCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGAc                  >  2:716444/2‑149 (MQ=255)
                                                                                                                                 taaaaaaaaGAGCGACCGAAGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCt                >  4:330319/5‑149 (MQ=255)
                                                                                                                                 aTCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCt                >  2:66113/1‑149 (MQ=255)
                                                                                                                                     aaaaaGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCAtt            >  3:204265/1‑149 (MQ=255)
                                                                                                                                      aaaaGAGCGACCGAAGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATtg           >  2:780785/1‑149 (MQ=255)
                                                                                                                                       aaaGAGCGACCGAAGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATtgt          >  2:72183/1‑149 (MQ=255)
                                                                                                                                        aaGAGCGACTTCGCACGATCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGTTTAAGAAAATTTTTGCCATTCAGGTCATCACACAAATCTCACCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATtgtg         >  4:715109/1‑149 (MQ=255)
                                                                                                                                         agagCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTgg        >  3:878022/1‑149 (MQ=255)
                                                                                                                                         agagCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTgg        <  1:750917/149‑1 (MQ=255)
                                                                                                                                           agCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTc      <  3:411129/149‑1 (MQ=255)
                                                                                                                                               aCTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTCGatc  >  2:471120/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAAAGCTGGTATCGCTCTGAACGACAACTTCGTGAAACTGGTATCCTGGTACGACAACGAAACCGGTTACTCCAACAAAGTTCTGGACCTGATCGCTCACATCTCCAAATAAGTTGAGATGACACTGTGATCTAAAAAGAGCGACTTCGGTCGCTCTTTTTTTTACCTGATAAAATGAAGTTAAAGGACTGCGTCATGATTAAGAAAATTTTTGCCCTTCCGGTCATCGAACAAATCTCCCCTGTCCTCTCCCGTCGTAAACTGGATGAACTGGACCTCATTGTGGTCGATC  >  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: