Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,539,241 T→C 5.0% K115E (AAA→GAA)  cysM ← cysteine synthase B (O‑acetylserine sulfhydrolase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,539,2410TC 5.0% 250.5 / 2.5 79K115E (AAA→GAA) cysMcysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base C (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CGCCGGTTTGCTGCCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACTGGTGGCTTC  >  NC_000913/2539094‑2539388
                                                                                                                                                   |                                                                                                                                                   
cgccgGTTTGCTGCCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTg                                                                                                                                                    <  6:140473/149‑1 (MQ=255)
cgccgGTTTGCTGCCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTg                                                                                                                                                    <  1:49187/149‑1 (MQ=255)
    ggTTTGCTGCCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGCTTATCGGGATTCTTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGAGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGa                                                                                                                                                <  2:286583/149‑1 (MQ=255)
     gTTTGCTGCCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCCCGATTCGCCATCTCCCGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGa                                                                                                                                                >  1:208528/1‑148 (MQ=255)
          ctgcCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCGTGCTCTTTGGTGACAAGa                                                                                                                                           <  5:158627/148‑1 (MQ=255)
          ctgcCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCGCTTTGGTGACAAGaa                                                                                                                                          <  8:141379/149‑1 (MQ=255)
            gcCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAAt                                                                                                                                         >  4:107370/1‑148 (MQ=255)
             ccAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGAAAAGAAt                                                                                                                                         >  7:142170/1‑147 (MQ=255)
              cAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAg                                                                                                                                      <  3:107370/149‑1 (MQ=255)
                      cGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCAc                                                                                                                              >  7:215009/1‑149 (MQ=255)
                        gCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCa                                                                                                                            <  8:215009/149‑1 (MQ=255)
                           cAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATag                                                                                                                         >  6:274460/1‑148 (MQ=255)
                            aGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAg                                                                                                                        <  5:117821/149‑1 (MQ=255)
                                gtggtgTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCAc                                                                                                                     <  7:164143/148‑1 (MQ=255)
                                 tggtgTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTCGGTGACAAGAATCAGTTCCGCACCATAAGCACg                                                                                                                    >  5:109158/1‑148 (MQ=255)
                                 tggggTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGc                                                                                                                   <  7:150378/149‑1 (MQ=255)
                                    tgtAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATc                                                                                                                >  2:125329/1‑149 (MQ=255)
                                      tAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCg                                                                                                                                 <  7:163025/130‑1 (MQ=255)
                                      tAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCg                                                                                                                                 >  8:163025/1‑130 (MQ=255)
                                      tAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATcgc                                                                                                              <  3:376773/149‑1 (MQ=255)
                                           cgcATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCAc                                                                                                         >  3:233777/1‑149 (MQ=255)
                                            gcATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACg                                                                                                        <  2:271249/149‑1 (MQ=255)
                                                   ggTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATAGCGCGCACCTTCCATGCCCTGCTCTTCGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTcc                                                                                                 <  6:109158/149‑1 (MQ=255)
                                                       aTCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTgg                                                                                              <  8:35601/148‑1 (MQ=255)
                                                       aTCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGc                                                                                             >  2:119353/1‑149 (MQ=255)
                                                          gggATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCa                                                                                          >  8:270211/1‑149 (MQ=255)
                                                           ggATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCAt                                                                                         >  2:316465/1‑149 (MQ=255)
                                                           ggATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCAt                                                                                         <  5:99856/149‑1 (MQ=255)
                                                           ggATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCAt                                                                                         >  1:69581/1‑149 (MQ=255)
                                                                attGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCAtgttg                                                                                     <  6:327273/148‑1 (MQ=255)
                                                                  tGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCg                                                                                  >  2:355079/1‑149 (MQ=255)
                                                                   gAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCgg                                                                                 >  2:109660/1‑149 (MQ=255)
                                                                       tGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCa                                                                                                                      >  3:268123/1‑108 (MQ=255)
                                                                       tGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCa                                                                                                                      <  4:268123/108‑1 (MQ=255)
                                                                         aTCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCt                                                                                                <  2:155878/128‑1 (MQ=255)
                                                                         aTCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCt                                                                                                >  1:155878/1‑128 (MQ=255)
                                                                         aTCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGTCTCATGTTGTCGGGCAt                                                                             <  3:3044/147‑1 (MQ=255)
                                                                                     tCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTcc                                                                                                 <  5:321004/115‑1 (MQ=255)
                                                                                     tCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTcc                                                                                                 >  6:321004/1‑115 (MQ=255)
                                                                                     tCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCAt                                                                 <  1:355079/147‑1 (MQ=255)
                                                                                     tCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCGCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATg                                                                <  1:263018/148‑1 (MQ=255)
                                                                                       cTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCTTAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGa                                                             >  5:366541/1‑149 (MQ=255)
                                                                                            tcACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGcc                                                        >  4:193182/2‑149 (MQ=255)
                                                                                            ccACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGcc                                                        >  3:181363/1‑149 (MQ=255)
                                                                                             cACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAAt                                                                                                                                         <  8:134135/67‑1 (MQ=255)
                                                                                             cACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGTAGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCt                                                       >  6:40487/1‑149 (MQ=255)
                                                                                             cACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGACAGCCt                                                       >  4:60357/1‑149 (MQ=255)
                                                                                              aCGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCtt                                                      >  1:11452/1‑149 (MQ=255)
                                                                                                gATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTc                                                    <  3:193182/149‑1 (MQ=255)
                                                                                                     gCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAgcgc                                               >  4:341739/1‑149 (MQ=255)
                                                                                                       cagctcCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATACGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCgg                                             <  3:343989/146‑1 (MQ=255)
                                                                                                         tctcCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCAtgttg                                                                                     <  4:339005/107‑1 (MQ=255)
                                                                                                         tctcCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCAtgttg                                                                                     >  3:339005/1‑107 (MQ=255)
                                                                                                         tctcCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGc                                            >  1:58665/1‑148 (MQ=255)
                                                                                                                  ccAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGc                                   >  7:315489/1‑148 (MQ=255)
                                                                                                                   cAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCCTAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCa                                 <  6:366541/149‑1 (MQ=255)
                                                                                                                   cAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGc                                   >  3:78319/1‑147 (MQ=255)
                                                                                                                          cgcgcACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTGCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAAt                          >  2:289612/1‑149 (MQ=255)
                                                                                                                               aCCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCt                                                       >  2:132509/1‑115 (MQ=255)
                                                                                                                               aCCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCt                                                       <  1:132509/115‑1 (MQ=255)
                                                                                                                               aCCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCgg                     >  5:245462/1‑149 (MQ=255)
                                                                                                                                ccTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGgt                    <  4:233777/149‑1 (MQ=255)
                                                                                                                                ccTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGgt                    >  1:229752/1‑149 (MQ=255)
                                                                                                                                 cTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACCGCGTTCCTGCCTCATGATGTCGGGCAACAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGgtg                   >  1:7767/1‑149 (MQ=255)
                                                                                                                                  ttCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCa                                 >  7:247295/1‑134 (MQ=255)
                                                                                                                                  ttCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCa                                 <  8:247295/134‑1 (MQ=255)
                                                                                                                                  ttCCATGCCCTGCTCTTCGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTTGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAgc                               >  5:166550/1‑136 (MQ=255)
                                                                                                                                  ttCCATGCCCTGCTCTTCGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTTGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAgc                               <  6:166550/136‑1 (MQ=255)
                                                                                                                                    ccATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTa                >  7:324643/1‑149 (MQ=255)
                                                                                                                                     cATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTa                <  1:119353/148‑1 (MQ=255)
                                                                                                                                         cccTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACt           <  2:229752/149‑1 (MQ=255)
                                                                                                                                          ccTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCgg                     <  3:16163/138‑1 (MQ=255)
                                                                                                                                          ccTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCgg                     >  4:16163/1‑138 (MQ=255)
                                                                                                                                          ccTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACtg          >  8:31207/1‑149 (MQ=255)
                                                                                                                                          ccTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACtg          >  8:248191/1‑149 (MQ=255)
                                                                                                                                          catGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACtg          >  6:176611/3‑149 (MQ=255)
                                                                                                                                              ctctTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATTCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACtggtgg      <  1:94884/149‑1 (MQ=255)
                                                                                                                                              ctctTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACtggtgg      <  7:31207/149‑1 (MQ=255)
                                                                                                                                              ctctTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACtggt        <  1:109660/147‑1 (MQ=255)
                                                                                                                                                  ttGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACTGGTGGCTTc  >  4:173618/1‑149 (MQ=255)
                                                                                                                                                   tGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACTGGTGGCTTc  <  2:11452/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCCGGTTTGCTGCCAGATTTCCGGCCCAGTGGTGGTGTAATGCGCATAAGGGTTATCGGGATTATTGAACTGATCGAGCAGCTTTCCTTCGCCACGATTCGCCATCTCCAGCGCCAGATCGCGCGCACCTTCCATGCCCTGCTCTTTGGTGACAAGAATCAGTTCCGCACCATAAGCACGCATCGCCGCACGGCGTTCCTGGCTCATGTTGTCGGGCATCAGCAATTTCATGCGATAGCCTTTCAGCGCGGCAATCATTGCCAGCGCAATGCCGGTGTTACCACTGGTGGCTTC  >  NC_000913/2539094‑2539388

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