Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,718,285 G→T 6.0% intergenic (+1/+49) ghrB → / ← yiaF glyoxylate reductase/DUF3053 domain‑containing protein YiaF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,718,2850GT 6.0% 545.9 / 12.9 160intergenic (+1/+49)ghrB/yiaFglyoxylate reductase/DUF3053 domain‑containing protein YiaF
Reads supporting (aligned to +/- strand):  ref base G (78/72);  new base T (7/3);  total (85/75)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 161 total aligned reads displayed.

CAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGT  >  NC_000913/3718136‑3718433
                                                                                                                                                     |                                                                                                                                                    
cAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAg                                                                                                                                                      >  1:653124/1‑150 (MQ=255)
 aaTGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGc                                                                                                                                                     <  1:491251/150‑1 (MQ=255)
 aaTGGCCAACGTCGTCGCAGTACCACATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGc                                                                                                                                                     >  2:352482/1‑150 (MQ=255)
      ccAACGGCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACCCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGTAAAGGTTGAGAAGATCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGa                                                                                                                                                <  2:1265993/150‑1 (MQ=255)
       cAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGAc                                                                                                                                               <  2:2097892/150‑1 (MQ=255)
          cgtcgtCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGCTTAGTCCGCGAcgtg                                                                                                                                            <  2:387802/150‑4 (MQ=255)
            tcgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGt                                                                                                                                          <  1:208617/150‑1 (MQ=255)
            tcgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGt                                                                                                                                          <  2:1756838/150‑1 (MQ=255)
            tcgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGt                                                                                                                                          >  1:1930518/1‑150 (MQ=255)
             cgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTa                                                                                                                                         <  2:1901402/150‑1 (MQ=255)
             cgtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTa                                                                                                                                         >  2:724225/1‑150 (MQ=255)
              gtcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTag                                                                                                                                        <  1:2155864/150‑1 (MQ=255)
               tcgCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTaga                                                                                                                                       <  2:1258536/150‑1 (MQ=255)
                     tACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAg                                                                                                                                 <  2:1225585/150‑1 (MQ=255)
                      aCCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGc                                                                                                                                >  2:1200333/1‑150 (MQ=255)
                       cctcATCTTGGATCTGCCACCCATGAGAGTCGTTATCGCATGGCCGCCTGTGCCCGGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTCTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGcc                                                                                                                               <  1:920273/147‑1 (MQ=255)
                        cGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGccc                                                                                                                              <  2:2468998/150‑1 (MQ=255)
                        cGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGccc                                                                                                                              >  1:735191/1‑150 (MQ=255)
                         gCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCg                                                                                                                             <  2:1053258/150‑1 (MQ=255)
                              ttGGATCTGCCACCCATGAGACGCGTTATGGCATGGTCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAAt                                                                                                                        >  1:397796/1‑150 (MQ=255)
                              ttGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAAt                                                                                                                        >  1:1515294/1‑150 (MQ=255)
                                 gATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGAt                                                                                                                           >  2:514484/1‑144 (MQ=255)
                                 gATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGc                                                                                                                     >  1:1115118/1‑150 (MQ=255)
                                     tGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCgg                                                                                                                 >  2:2382609/1‑150 (MQ=255)
                                     tGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCgg                                                                                                                 >  1:318419/1‑150 (MQ=255)
                                        cACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCt                                                                                                              >  1:671586/1‑150 (MQ=255)
                                             aTGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTAc                                                                                                         <  1:1371268/150‑1 (MQ=255)
                                               gagaCGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTActc                                                                                                       >  1:2450091/1‑150 (MQ=255)
                                                  aCGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCttt                                                                                                    <  2:1515294/150‑1 (MQ=255)
                                                   cgcgTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTa                                                                                                   >  1:1304475/1‑150 (MQ=255)
                                                    gcgTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTAt                                                                                                  <  1:875812/150‑1 (MQ=255)
                                                    gcgTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTAt                                                                                                  >  2:1899520/1‑150 (MQ=255)
                                                      gTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTg                                                                                                >  1:221828/1‑150 (MQ=255)
                                                       ttATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTgg                                                                                               <  2:1061097/150‑1 (MQ=255)
                                                       ttATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTgg                                                                                               <  2:320260/150‑1 (MQ=255)
                                                         atggcatggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGt                                                                                             <  1:152212/150‑1 (MQ=255)
                                                           ggcatggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTg                                                                                           >  1:660450/1‑150 (MQ=255)
                                                            gcatggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTGGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGc                                                                                          <  1:1971958/150‑1 (MQ=255)
                                                              atggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAg                                                                                        <  2:372759/150‑1 (MQ=255)
                                                              atggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAg                                                                                        <  2:1355334/150‑1 (MQ=255)
                                                               tggcCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGt                                                                                       <  1:352482/150‑1 (MQ=255)
                                                                 gccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTaa                                                                                     >  2:981919/1‑150 (MQ=255)
                                                                 gccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTaa                                                                                     >  1:1956955/1‑150 (MQ=255)
                                                                  ccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAc                                                                                    <  1:1394076/150‑1 (MQ=255)
                                                                  ccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAc                                                                                    <  1:782507/150‑1 (MQ=255)
                                                                  ccgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAc                                                                                    <  2:802641/150‑1 (MQ=255)
                                                                   cgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAct                                                                                   >  1:1055959/1‑150 (MQ=255)
                                                                   cgccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAAct                                                                                   <  2:922886/150‑1 (MQ=255)
                                                                     ccTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCAGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAActgc                                                                                 >  1:1114708/1‑150 (MQ=255)
                                                                        gtgCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGt                                                                              <  2:188369/150‑1 (MQ=255)
                                                                        gagCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGt                                                                              >  2:741565/3‑150 (MQ=255)
                                                                           ccGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGt                                                                           >  1:1304189/1‑150 (MQ=255)
                                                                            cGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGCGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTc                                                                          <  1:28764/150‑1 (MQ=255)
                                                                                gATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAgg                                                                      >  2:463308/1‑150 (MQ=255)
                                                                                  tAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGcc                                                                    >  1:1181013/1‑150 (MQ=255)
                                                                                    aTTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCAGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCAGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTg                                                                  >  1:2084354/1‑150 (MQ=255)
                                                                                      ttgattgaTGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGAt                                                                <  2:1853229/150‑1 (MQ=255)
                                                                                      ttgattgaTGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGAt                                                                >  2:297052/1‑150 (MQ=255)
                                                                                       tgattgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGAtt                                                               <  1:1329609/150‑1 (MQ=255)
                                                                                        gattgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTa                                                              >  2:59747/1‑150 (MQ=255)
                                                                                         attgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTaa                                                             >  1:454936/1‑150 (MQ=255)
                                                                                           tgatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAAc                                                           >  2:1592155/1‑150 (MQ=255)
                                                                                            gatGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACg                                                          >  1:2025768/1‑150 (MQ=255)
                                                                                              tGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGcc                                                        <  2:2214910/150‑1 (MQ=255)
                                                                                              tGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGcc                                                        >  1:213506/1‑150 (MQ=255)
                                                                                                cGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTg                                                      >  2:1978224/1‑150 (MQ=255)
                                                                                                 gTTACAAGTAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGGTGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTTCTGTAGTCCAGGCCTGATTAAACGCCTGa                                                     <  2:2084354/150‑1 (MQ=255)
                                                                                                 gTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGa                                                     >  2:1004431/1‑150 (MQ=255)
                                                                                                 gTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGa                                                     >  1:1684448/1‑150 (MQ=255)
                                                                                                   tACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATg                                                   >  1:2173408/1‑150 (MQ=255)
                                                                                                   tACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATg                                                   >  1:440280/1‑150 (MQ=255)
                                                                                                    aCAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGc                                                  <  1:575788/150‑1 (MQ=255)
                                                                                                       aGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCtgt                                               >  2:191687/1‑150 (MQ=255)
                                                                                                            aGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCgg                                          <  2:1477257/150‑1 (MQ=255)
                                                                                                             ggTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGt                                         >  2:768769/1‑150 (MQ=255)
                                                                                                              gTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTa                                        >  1:524136/1‑150 (MQ=255)
                                                                                                                tGAGACGAACTGTGTGAATCCGCACGTCGCCGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCGCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAAt                                      <  2:2346426/150‑1 (MQ=255)
                                                                                                                tGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAAt                                      >  1:1401365/1‑150 (MQ=255)
                                                                                                                   gaagaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgc                                   >  1:2169962/1‑150 (MQ=255)
                                                                                                                    aagaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgcg                                  <  1:724225/150‑1 (MQ=255)
                                                                                                                     agaaCTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGgcgc                                 <  2:221828/150‑1 (MQ=255)
                                                                                                                           gtgtgAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAg                           <  1:1332568/150‑1 (MQ=255)
                                                                                                                           gtgtgAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAg                           <  1:896523/150‑1 (MQ=255)
                                                                                                                           gtgtgAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAg                           >  1:197134/1‑150 (MQ=255)
                                                                                                                             gtgAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGtt                         >  2:2091044/1‑150 (MQ=255)
                                                                                                                               gAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTg                       >  2:929569/1‑150 (MQ=255)
                                                                                                                                aaTCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGt                      <  2:1930716/150‑1 (MQ=255)
                                                                                                                                  ctacgcagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTta                    >  2:874189/8‑150 (MQ=255)
                                                                                                                                     gCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTatat                 <  2:302966/150‑1 (MQ=255)
                                                                                                                                     cgcagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTatat                 >  2:2189463/5‑150 (MQ=255)
                                                                                                                                      gcagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATAtt                >  1:2277478/4‑150 (MQ=255)
                                                                                                                                       cagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATAtta               >  2:217063/3‑149 (MQ=255)
                                                                                                                                       cagTCGCGGCTTAGTCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTc               >  2:1212006/3‑150 (MQ=255)
                                                                                                                                          tCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACt            >  1:886002/1‑150 (MQ=255)
                                                                                                                                            gcgGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTg          >  1:1103844/1‑150 (MQ=255)
                                                                                                                                              ggACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGcc        <  2:2412424/150‑1 (MQ=255)
                                                                                                                                                aCTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCctg      >  2:1727987/1‑150 (MQ=255)
                                                                                                                                                   aaGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCctgctg   >  2:440246/1‑150 (MQ=255)
                                                                                                                                                   aaGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCctgctg   <  1:2162213/150‑1 (MQ=255)
                                                                                                                                                    aGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCctgctgg  >  1:2289241/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAATGGCCAACGTCGTCGCAGTACCGCATATTGGATCTGCCACCCATGAGACGCGTTATGGCATGGCCGCCTGTGCCGTGGATAATTTGATTGATGCGTTACAAGGAAAGGTTGAGAAGAACTGTGTGAATCCGCACGTCGCGGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGT  >  NC_000913/3718136‑3718433

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