Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,869,015 C→T 6.0% R21C (CGT→TGT)  yeaH → UPF0229 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,869,0150CT 6.0% 211.3 / 2.8 68R21C (CGT→TGT) yeaHUPF0229 family protein
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base T (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGCGTGGC  >  NC_000913/1868874‑1869161
                                                                                                                                             |                                                                                                                                                  
cTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGtttttt                                                                                                                                             <  2:14808/149‑1 (MQ=255)
cTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGttttt                                                                                                                                              >  2:351597/1‑148 (MQ=255)
   cAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTa                                                                                                                                            <  3:255331/147‑1 (MQ=255)
       gACAGCGAACCGTGGGCCTGAGAAGCGCCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAACCAACGTGAATCGCCCGCGTTTTTTACGCCga                                                                                                                                      >  6:165903/1‑148 (MQ=255)
        aCAGCGAACCGTGGGCCGGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAACGCATGGTGAATCGCCAGCGTTTTTTACGCCGtt                                                                                                                                     <  5:332937/149‑1 (MQ=255)
           gcgcACCGTGGGCATGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATTTGACCTGGTTTATTGACCGGCCGCTGAACGGCAAAAACAAAACCTAGGTGAATCGCCAGAGTTTTTTACGCCGTtat                                                                                                                                   >  4:110729/1‑148 (MQ=255)
            cGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGAATACATTTGGCAAATTTAGAAATGGGGTCATATCACATGTTTTATTGACCGGCGTCTTCACGGCAAACACAAAATCATGGTTAAACGCCAGCGTTTTTTACGCATTtata                                                                                                                                  >  2:316240/1‑148 (MQ=255)
                       ccTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGCACGGGAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAAtt                                                                                                                       >  3:142721/1‑148 (MQ=255)
                        cTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGTCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTACGCCGTTATAAAGCGCAAATTaaa                                                                                                                    <  1:225003/149‑1 (MQ=255)
                          gagaAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAAc                                                                                                                   <  4:64311/149‑1 (MQ=255)
                           agaAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTaaa                                                                                                                    <  5:165903/147‑1 (MQ=255)
                           agaAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTaaa                                                                                                                    <  1:118732/147‑1 (MQ=255)
                             aaGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGt                                                                                                                <  8:128564/149‑1 (MQ=255)
                               gCGGCAACACAGGCGTAGGATAGAGTTGGCAAATGTAGTACGGGCGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGAAAAAACAAAAGCATGTTGAATCGCCATAGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTc                                                                                                               <  3:285829/148‑1 (MQ=255)
                               gCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGCCCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTc                                                                                                               >  5:261096/1‑148 (MQ=255)
                                       acaGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGa                                                                                                      <  1:351597/149‑1 (MQ=255)
                                         aGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAgg                                                                                                    >  8:249855/1‑149 (MQ=255)
                                              tAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCAt                                                                                                >  6:280798/1‑148 (MQ=255)
                                                 cATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATtaat                                                                                            <  2:343213/149‑1 (MQ=255)
                                                  aTACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATtaat                                                                                            >  2:13061/1‑148 (MQ=255)
                                                  aTACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATtaat                                                                                            >  3:259261/1‑148 (MQ=255)
                                                   tACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCAAGCGTTTATTACGCCGTTATAAAGCGCAAATAAAACAGTc                                                                                                               >  4:252309/1‑128 (MQ=255)
                                                        ttGGCAAATGTAGTACGGGGGGCAGATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGtt                                                                                     <  3:9051/149‑1 (MQ=255)
                                                           tcAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCg                                                                                   <  7:375527/147‑1 (MQ=255)
                                                                              cATATGACCTGGTTTTTTGACCTGCGTCTGATCGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCTTTATAAAGCGCAATTTAAAGAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGTCTTCGACAgcgg                                                               <  6:138074/149‑1 (MQ=255)
                                                                                tatGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAgcggc                                                              >  8:354142/1‑148 (MQ=255)
                                                                                tatGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCAATAATAAGCGTTCGGTGACTGACGTCGACAgcggc                                                              >  1:333771/1‑148 (MQ=255)
                                                                                     ccTGGTTTATTGACGGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATcc                                                        <  3:268955/149‑1 (MQ=255)
                                                                                     ccTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATcc                                                        <  4:60581/149‑1 (MQ=255)
                                                                                      cTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGTGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCg                                                       >  3:120858/1‑149 (MQ=255)
                                                                                      cTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGTGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCg                                                       >  3:125771/1‑149 (MQ=255)
                                                                                      cTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATCGTGAATCGCCAGCGTTATTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATAAATAAGCGTTCGGAGACTGACGTCGACAGCGGCGAATCCg                                                       >  1:336513/1‑149 (MQ=255)
                                                                                      cTGGTTTATTGACCGGCGTCTGAACGGCAAAAAAAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAAAAGTCGATCTCCGCGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCg                                                       >  6:52708/1‑149 (MQ=255)
                                                                                         gTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTAt                                                    <  4:135152/149‑1 (MQ=255)
                                                                                          tttATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATc                                                   >  1:211399/1‑149 (MQ=255)
                                                                                             aTTGACCGGGGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACTCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCa                                                 <  4:37982/148‑1 (MQ=255)
                                                                                             aTTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCa                                                 <  7:249855/148‑1 (MQ=255)
                                                                                              ttGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCAtt                                               <  1:13061/149‑1 (MQ=255)
                                                                                              ttGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCAtt                                               <  5:280798/149‑1 (MQ=255)
                                                                                               tGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTc                                              >  2:29705/1‑149 (MQ=255)
                                                                                                 aCCGGCGTCTGAACGGCAAAAACAAAATAATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTAGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTcc                                             >  2:238151/1‑148 (MQ=255)
                                                                                                 aCCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTc                                                                       >  2:64919/1‑122 (MQ=255)
                                                                                                 aCCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTc                                                                       <  1:64919/122‑1 (MQ=255)
                                                                                                 aCCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTccc                                            >  3:12403/1‑149 (MQ=255)
                                                                                                 aCCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTcc                                             >  3:71272/1‑148 (MQ=255)
                                                                                                       gTCTGAACGGCAAAAACAAAAGCCTGGTGAAGCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGaa                                      <  2:288899/149‑1 (MQ=255)
                                                                                                                gCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGc                             >  1:306544/1‑149 (MQ=255)
                                                                                                                gCAAAAAAAAAAGCATGGTGAATCGCAAGCGTTTTTTAAGACGTTATAAAGCGCAAATTAAACAGTCGTTCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGc                             <  1:316240/149‑1 (MQ=255)
                                                                                                                  aaaaaCAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCg                            >  7:173331/1‑148 (MQ=255)
                                                                                                                      aCAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGCGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGATcc                        >  2:350029/1‑148 (MQ=255)
                                                                                                                      aCAAAAGCATGGTGAATCGCCAGCGTATATTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGTTTACTAACGTAGAAAGCGGCGAATCAGTAACCAAACCAACGGaa                                      >  5:18069/1‑134 (MQ=255)
                                                                                                                        aaaaGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGAt                     >  8:309082/1‑149 (MQ=255)
                                                                                                                         aaaGCATGGTGAATCGCCAGTGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATg                    <  4:125771/149‑1 (MQ=255)
                                                                                                                         aaaGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAAt                                                          <  1:203893/111‑1 (MQ=255)
                                                                                                                         aaaGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAAt                                                          >  2:203893/1‑111 (MQ=255)
                                                                                                                          aaGCATGGTGAATCGCCAGTTTTTTTTAAGCCGTTATAATGTGCTAATTAAACTGTCGATCTCCGAGGCCTTTAATTAGCGTTGGGTGACTGTCGTCGACAGCGGCGAATCCGTATCCATTCCCACGTATGATATTAGCGTACCGATg                    <  4:120858/148‑1 (MQ=255)
                                                                                                                               tGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTtcat              >  7:330846/1‑149 (MQ=255)
                                                                                                                                ggTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTtcat              >  1:186259/1‑148 (MQ=255)
                                                                                                                                  tGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAg           >  5:158425/1‑149 (MQ=255)
                                                                                                                                    aaTCGCCCGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGCGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCTAACCGATGTTTCATCAggg         <  6:261096/149‑1 (MQ=255)
                                                                                                                                      tCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCTACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGc       >  6:124652/1‑149 (MQ=255)
                                                                                                                                      tCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGc       >  4:226658/1‑149 (MQ=255)
                                                                                                                                      tCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGc       >  7:352841/1‑149 (MQ=255)
                                                                                                                                         ccAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTa                                                     >  4:124287/1‑100 (MQ=255)
                                                                                                                                         ccAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTa                                                     <  3:124287/100‑1 (MQ=255)
                                                                                                                                          cAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGTATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGCGTgg   >  3:17657/1‑149 (MQ=255)
                                                                                                                                          cAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTCGCGAACCGATGTTTCATCAGGGGCGTgg   <  4:259261/149‑1 (MQ=255)
                                                                                                                                           aGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGCGTGGc  <  7:321559/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGGTGAATCGCCAGCGTTTTTTACGCCGTTATAAAGCGCAAATTAAACAGTCGATCTCCGAGGCCATTAATAAGCGTTCGGTGACTGACGTCGACAGCGGCGAATCCGTATCCATTCCCACGGAAGATATTAGCGAACCGATGTTTCATCAGGGGCGTGGC  >  NC_000913/1868874‑1869161

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