Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,868,859 Δ1 bp 22.2% intergenic (+17/‑96) yeaG → / → yeaH protein kinase, endogenous substrate unidentified; autokinase/UPF0229 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,868,8590C.22.2% 241.7 / 70.1 90intergenic (+17/‑96)yeaG/yeaHprotein kinase, endogenous substrate unidentified; autokinase/UPF0229 family protein
Reads supporting (aligned to +/- strand):  ref base C (39/31);  new base . (9/11);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGA  >  NC_000913/1868760‑1868958
                                                                                                   |                                                                                                    
aCCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTac                                                                                                      <  2:1943009/100‑1 (MQ=255)
aCCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTac                                                                                                      >  2:2304563/1‑100 (MQ=255)
 ccGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTaca                                                                                                     <  2:76075/100‑1 (MQ=255)
    tATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACAt                                                                                                  <  2:874346/100‑1 (MQ=255)
     atgatgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATg                                                                                                 >  2:2468963/1‑100 (MQ=255)
       gatgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGtt                                                                                              >  1:910773/1‑100 (MQ=255)
        atgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTa                                                                                             >  2:372370/1‑100 (MQ=255)
         tgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTAc                                                                                            >  2:576941/1‑100 (MQ=255)
           naaaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCg                                                                                          >  1:1531075/2‑100 (MQ=255)
             aaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCggg                                                                                        <  2:2563121/100‑1 (MQ=255)
              tctggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGc                                                                                       <  1:2098237/97‑1 (MQ=255)
              aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGc                                                                                       >  1:486337/1‑100 (MQ=255)
              aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGc                                                                                       <  2:448368/100‑1 (MQ=255)
              aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCggg                                                                                        <  1:1425417/99‑1 (MQ=255)
              aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGc                                                                                       <  1:2360157/100‑1 (MQ=255)
                aGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGcca                                                                                     <  2:1313878/100‑2 (MQ=39)
                 ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTa                                                                                    <  2:2085919/100‑1 (MQ=255)
                 ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGccag                                                                                    <  2:2104057/100‑3 (MQ=39)
                 ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGcaga                                                                                    >  2:2098251/1‑97 (MQ=255)
                  gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTAc                                                                                   <  2:562952/100‑1 (MQ=255)
                  gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGccagg                                                                                   <  2:1453311/100‑4 (MQ=38)
                  gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGccag                                                                                    <  2:2360815/99‑3 (MQ=38)
                   cTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACCTGTTACCGGGCCTaaa                                                                                  >  1:1249426/1‑98 (MQ=255)
                    tACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGTGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGccagggt                                                                                 <  2:2493163/100‑6 (MQ=37)
                    tACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGccagggt                                                                                 <  2:132415/100‑6 (MQ=38)
                     acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACCAc                                                                                >  1:2031577/1‑100 (MQ=255)
                     acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAAc                                                                                <  2:2132023/100‑1 (MQ=255)
                     acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTCGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAACATGTGTACCGGGccagggt                                                                                 <  2:783663/99‑6 (MQ=37)
                      cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACg                                                                               >  2:91665/1‑100 (MQ=255)
                      cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCACGGTACACATGTTACCGGGCCTACAACg                                                                               >  1:1064488/1‑100 (MQ=255)
                      cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGAGCGTCCGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACCACg                                                                               >  1:1149676/1‑100 (MQ=255)
                        cccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGAc                                                                             >  2:1016016/1‑100 (MQ=255)
                        cccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGAc                                                                             >  2:531575/1‑100 (MQ=255)
                         ccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACa                                                                            >  2:1925694/1‑100 (MQ=255)
                             aaaCAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGa                                                                        >  1:1867197/1‑100 (MQ=255)
                              aaCAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGaa                                                                       >  1:34414/1‑100 (MQ=255)
                                 aGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCg                                                                    <  1:714808/100‑1 (MQ=255)
                                  ggTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCg                                                                    >  2:2359081/1‑99 (MQ=255)
                                  ggTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCg                                                                    <  1:1128414/99‑1 (MQ=255)
                                    tGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTgg                                                                 <  1:333702/100‑1 (MQ=255)
                                     gCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTggg                                                                <  2:1460080/100‑1 (MQ=255)
                                      cGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGc                                                               <  1:667650/100‑1 (MQ=255)
                                       gTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGcc                                                              <  1:1316312/100‑1 (MQ=255)
                                       gTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGcc                                                              <  2:699804/100‑1 (MQ=255)
                                         ttACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTg                                                            <  1:2473489/100‑1 (MQ=255)
                                          tACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGcc                                                              >  1:2024308/1‑97 (MQ=255)
                                           aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAGCGACAGCGAACCGTGGGCCTgag                                                          <  1:774688/100‑1 (MQ=255)
                                            cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTGACCGGGCCTACAACGACAGCGAACCGTGGGCCTgaga                                                         <  1:1389094/100‑1 (MQ=255)
                                            cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGATACCGGGCCTACAACGACAGCGAACCGTGGGCATgag                                                          >  1:688329/1‑99 (MQ=255)
                                              gtgCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAg                                                       >  2:949295/1‑100 (MQ=255)
                                                 cGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCgg                                                    >  2:989774/1‑100 (MQ=255)
                                                  gAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGc                                                   <  2:1637298/100‑1 (MQ=255)
                                                   aaTGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGc                                                   <  1:2072693/99‑1 (MQ=255)
                                                   aaTGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCa                                                  <  1:944057/100‑1 (MQ=255)
                                                    aTGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCaa                                                 >  1:1987875/1‑100 (MQ=255)
                                                      ggTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAaca                                               >  2:85142/1‑100 (MQ=255)
                                                      ggTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAaca                                               >  2:492509/1‑100 (MQ=255)
                                                        tATTTGCGCGTAGGTAAATCGTCTCAACAACCCGGGCCCGGTACACACGTTACAGGGCCTACTACGATAGCGAACCGTGGGCCTGAGATGCGGCAac                                                <  1:1152545/97‑1 (MQ=255)
                                                        tATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAAcaca                                             <  2:486339/100‑1 (MQ=255)
                                                         aTTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAg                                            <  2:1064492/100‑1 (MQ=255)
                                                           tttcgcgTACGTAAATCGTCTTAACAACCCTGGCCCTGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGACAAGCGGCAACACAGGc                                          >  1:130637/4‑100 (MQ=255)
                                                            tGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCg                                         >  1:343040/1‑100 (MQ=255)
                                                             gcgcGTACGAAAATCGGCTTAACAACCCTGGCGCGGTACACATGATACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGGGGCAACACAgg                                           >  1:863986/1‑97 (MQ=255)
                                                                    cgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATAc                                 >  2:642951/1‑100 (MQ=255)
                                                                     gtaAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACa                                >  2:2241710/1‑100 (MQ=255)
                                                                     gtaAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACa                                >  2:488216/1‑100 (MQ=255)
                                                                     gtaAATCGTCTTAACAACCCTGGCACGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACa                                <  1:1025376/100‑1 (MQ=255)
                                                                      taAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAg                               <  1:1086117/100‑1 (MQ=255)
                                                                       aaaTCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGt                              >  2:828887/1‑100 (MQ=255)
                                                                           cGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGc                          >  1:1012590/1‑100 (MQ=255)
                                                                            gTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCa                         <  2:1993035/100‑1 (MQ=255)
                                                                            gTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCa                         <  1:1762085/100‑1 (MQ=255)
                                                                             tCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATCCAGTTGGCaa                        >  2:1171803/1‑100 (MQ=255)
                                                                              cTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCaaa                       >  1:1394375/1‑100 (MQ=255)
                                                                               ttAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCaact                      >  1:1571304/1‑98 (MQ=255)
                                                                                tAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGCCCGCGAACCGTGGGCCTGAGAAGCGGCCACACAGGCGTCGCATACAGTTGGCaaagg                     >  1:347899/1‑98 (MQ=255)
                                                                                tAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATg                     >  2:792941/1‑100 (MQ=255)
                                                                                 aacaacCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgt                    <  1:1058034/100‑1 (MQ=255)
                                                                                  acaacCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgta                   <  2:2510039/100‑1 (MQ=255)
                                                                                     acCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgtagta                >  2:2031246/1‑100 (MQ=255)
                                                                                     acCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgtagta                <  1:942591/100‑1 (MQ=255)
                                                                                       gtaggCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACg              >  2:2265158/4‑100 (MQ=38)
                                                                                         aggCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgttcaatcag            >  2:1949632/2‑92 (MQ=38)
                                                                                           gCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  1:344538/1‑100 (MQ=255)
                                                                                           gCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  2:763892/1‑100 (MQ=255)
                                                                                           gCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  2:245938/1‑100 (MQ=255)
                                                                                           gCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  2:1922072/1‑100 (MQ=255)
                                                                                           gCCCGGTAACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  2:918545/1‑100 (MQ=255)
                                                                                            cccGGTAACATGTGTACCGGGCCTACAACGACCGCGAACCGTGGGCCTGAGACGCGGCATCACAGGCGTAGCATACAGCTGCCACATGTAGTACggggg          >  1:2205765/1‑99 (MQ=255)
                                                                                             ccGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGc        <  2:2140856/100‑1 (MQ=255)
                                                                                               ggTACACATGTTACCGGGCCTACAACGCCAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCat      >  2:680203/1‑100 (MQ=255)
                                                                                                gTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCata     <  2:34415/100‑1 (MQ=255)
                                                                                                gTACACATGTGACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCata     <  1:1978934/100‑1 (MQ=255)
                                                                                                  acacATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATg   >  2:613843/1‑100 (MQ=255)
                                                                                                  acacATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATg   <  2:2532374/100‑1 (MQ=255)
                                                                                                   cacaTGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGa  >  1:1220018/1‑100 (MQ=255)
                                                                                                   cacaTGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGa  >  2:34439/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
ACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGTTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGA  >  NC_000913/1868760‑1868958

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