Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,132,430 T→A 7.0% intergenic (+163/‑186) metL → / → metF Bifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,132,4300TA 7.0% 491.6 / 16.4 169intergenic (+163/‑186)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (88/69);  new base A (4/8);  total (92/77)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 187 total aligned reads displayed.

AATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAAC  >  NC_000913/4132283‑4132579
                                                                                                                                                   |                                                                                                                                                     
 atGCCCTGTACTAATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                                                                                                                                    <  1:1364731/150‑1 (MQ=255)
 atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                                                                                                                                    <  2:182054/150‑1 (MQ=255)
  tGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAt                                                                                                                                                   <  1:758496/150‑1 (MQ=255)
    cccGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAc                                                                                                                                                 >  1:657305/1‑150 (MQ=255)
     ccGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAcc                                                                                                                                                <  2:819215/150‑1 (MQ=255)
       ggTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCgg                                                                                                                                              <  2:690656/150‑1 (MQ=255)
       ggTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCgg                                                                                                                                              <  2:936939/150‑1 (MQ=255)
       ggTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGAGTTTACGACGAATGCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCgg                                                                                                                                              <  2:1277341/150‑1 (MQ=255)
        gTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGa                                                                                                                                             <  2:741269/150‑1 (MQ=255)
        gTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGa                                                                                                                                             <  2:704991/150‑1 (MQ=255)
         tACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGac                                                                                                                                            <  2:1574927/150‑2 (MQ=255)
         tACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGac                                                                                                                                            <  2:333114/150‑2 (MQ=255)
           cTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                                                                                                                                          >  2:1245027/1‑150 (MQ=255)
            tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                                                                                                                                         >  1:1580982/1‑150 (MQ=255)
             cATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                                                                                                                                        >  2:834502/1‑150 (MQ=255)
             cATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                                                                                                                                        >  2:1667524/1‑150 (MQ=255)
             cATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                                                                                                                                        >  1:1082348/1‑150 (MQ=255)
               tGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGAAGTCCGGTTTGATTTGCATACCGGAttttctt                                                                                                                                      >  1:1222941/1‑150 (MQ=255)
                       ggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc                                                                                                                              >  2:1457186/1‑150 (MQ=255)
                        gTTTATGGTTTCTAATGAAATATATTGAATTATCTTAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc                                                                                                                             >  1:990922/1‑150 (MQ=255)
                        gTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTATTATTCTTAcc                                                                                                                             >  2:266504/1‑150 (MQ=255)
                           tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATTTTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                          >  2:1050397/1‑150 (MQ=255)
                           tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                          >  2:501368/1‑150 (MQ=255)
                           tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGGTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAACAta                                                                                                                          >  2:1512541/1‑149 (MQ=255)
                            aTGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATcc                                                                                                                         <  1:1403056/150‑1 (MQ=255)
                             tGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCt                                                                                                                        >  2:796403/1‑150 (MQ=255)
                               gTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGa                                                                                                                      >  2:46023/1‑150 (MQ=255)
                               gTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGGTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGa                                                                                                                      <  1:715790/150‑1 (MQ=255)
                                 ttCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAg                                                                                                                    <  1:124995/150‑1 (MQ=255)
                                   cTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtt                                                                                                                  >  2:1479018/1‑150 (MQ=255)
                                     aaTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttt                                                                                                                <  1:1208967/150‑1 (MQ=255)
                                        taaaTATATTGAATTATCATATGATTAGTCCGGATTAAGCGTTTACGACGTATCCGGCTAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTc                                                                                                             <  1:253169/149‑1 (MQ=255)
                                            tatatTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                                                                                                                                  >  1:399437/1‑125 (MQ=255)
                                            tatatTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCt                                                                                                         >  2:54137/1‑150 (MQ=255)
                                               atTTAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTcc                                                                                                      <  1:785018/150‑1 (MQ=255)
                                                ttGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTccc                                                                                                     <  1:1379249/150‑1 (MQ=255)
                                                  gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATCCCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                                   >  2:1135261/1‑150 (MQ=255)
                                                  gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTGTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                                   >  1:87476/1‑150 (MQ=255)
                                                     ttATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGAtt                                                                                                <  1:306929/150‑1 (MQ=255)
                                                         cATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTc                                                                                            >  2:1664235/1‑150 (MQ=255)
                                                         cATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTc                                                                                            >  2:913300/1‑150 (MQ=255)
                                                          aTAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTcc                                                                                           >  2:238453/1‑150 (MQ=255)
                                                             ggATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCa                                                                                        <  1:466887/150‑1 (MQ=255)
                                                             ggATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCa                                                                                        >  2:717261/1‑150 (MQ=255)
                                                              gATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAc                                                                                       >  2:529616/1‑150 (MQ=255)
                                                              gATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAc                                                                                       <  2:881555/150‑1 (MQ=255)
                                                                     ccGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACAATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCAtt                                                                                >  1:853759/1‑150 (MQ=255)
                                                                       ggATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTgg                                                                              >  1:1036037/1‑150 (MQ=255)
                                                                           tAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCAt                                                                          <  1:1493122/150‑1 (MQ=255)
                                                                            aaGAGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCAtt                                                                         <  2:472485/150‑1 (MQ=255)
                                                                              gCGTTTACGACGAATCCGGCAAGAAGCAATATGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCAtttt                                                                       >  2:478141/1‑150 (MQ=255)
                                                                               cGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCAttttt                                                                      >  1:1593682/1‑150 (MQ=255)
                                                                                      gacgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGa                                                               >  2:1337240/1‑150 (MQ=255)
                                                                                           aTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGccc                                                          >  1:947981/1‑150 (MQ=255)
                                                                                            tCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAACATAATGAAGTTTTTTCATCTTCCCTGATTTTTACTCACCATCATTGGTCATTTTTCGGTTGACTCCCt                                                         >  1:75016/1‑150 (MQ=255)
                                                                                              ctgCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTTCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGAATTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                       <  2:560180/148‑1 (MQ=255)
                                                                                              cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCtta                                                       >  1:477009/1‑149 (MQ=255)
                                                                                               ggCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                       <  1:1528976/149‑1 (MQ=255)
                                                                                               ggCAAGAAGCAATAACTACATGGTTAGTTTATATTTGCAGTCAGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCg                                                      <  1:945269/150‑1 (MQ=255)
                                                                                                 cAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGc                                                    <  2:1261107/150‑1 (MQ=255)
                                                                                                  aagaagCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCt                                                   <  1:216876/150‑1 (MQ=255)
                                                                                                     aagCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCtttt                                                <  2:1546560/150‑1 (MQ=255)
                                                                                                       gCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTcc                                              >  2:381302/1‑150 (MQ=255)
                                                                                                        cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                             >  2:1200478/1‑150 (MQ=255)
                                                                                                        cAATAAGTACATGGTTAGTTTATATATGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACAATCCTGAAGTTTTTTCATCTTCCCCGATTTTTACTCACCATCATTGGTCATTTTTCGGTTTACTCCATTCGGCTTTTCCt                                             >  1:245765/1‑150 (MQ=255)
                                                                                                         aaTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCCTTTTTCGGTTGACGCCCTTCGGCTTTTCCtt                                            <  1:1135261/150‑1 (MQ=255)
                                                                                                         aaTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCtt                                            >  2:1436768/1‑150 (MQ=255)
                                                                                                         aaTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTAGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCtt                                            >  2:205570/1‑150 (MQ=255)
                                                                                                          aTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTc                                           <  1:834502/150‑1 (MQ=255)
                                                                                                             aGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATc                                        <  1:1003950/150‑1 (MQ=255)
                                                                                                             aGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATc                                        >  1:986241/1‑150 (MQ=255)
                                                                                                               tACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCtt                                      >  2:1213028/1‑150 (MQ=255)
                                                                                                                aCATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCttt                                     >  1:1696177/1‑150 (MQ=255)
                                                                                                                   tGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACa                                  >  2:338726/1‑150 (MQ=255)
                                                                                                                   tGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACa                                  >  2:366470/1‑150 (MQ=255)
                                                                                                                        aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                             >  1:1104757/1‑150 (MQ=255)
                                                                                                                        aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATATTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                             >  1:583208/1‑150 (MQ=255)
                                                                                                                          tttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGAc                           <  2:412294/150‑1 (MQ=255)
                                                                                                                          tttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACg                          <  1:723488/150‑1 (MQ=255)
                                                                                                                           ttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACg                          <  1:913300/150‑1 (MQ=255)
                                                                                                                            tataTTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGAAGt                         >  1:1259792/1‑150 (MQ=255)
                                                                                                                             atatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                        <  2:1082348/150‑1 (MQ=255)
                                                                                                                             atatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                        <  1:46023/150‑1 (MQ=255)
                                                                                                                              tatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACTCCCTTCGGCTTTTCCTACATCTTTACATCTGGACGTCt                       >  2:779744/1‑150 (MQ=255)
                                                                                                                                aagaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTATGTTGACGCCCTTCGGCTTTTCCTTAATCTTTAAATCTGGACGTCTaa                     >  1:1709310/8‑150 (MQ=255)
                                                                                                                                aagaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATACTGAAGTTTTTTCATCTTCCCTGATTTTTGCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCGTTTCCTTCATCTTTACATCTGGACGTCTaa                     >  1:1248880/8‑150 (MQ=255)
                                                                                                                                  tGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAAc                   <  2:1369992/150‑1 (MQ=255)
                                                                                                                                    cAGTCCGGTTTGCTTTGCATACCGGATGTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACgg                 >  1:829561/1‑150 (MQ=255)
                                                                                                                                      atCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGAt               >  1:504055/2‑150 (MQ=255)
                                                                                                                                       tCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATa              <  2:863500/150‑1 (MQ=255)
                                                                                                                                       tCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATa              <  1:1269441/150‑1 (MQ=255)
                                                                                                                                        ccGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCACTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAg             >  1:38782/1‑150 (MQ=255)
                                                                                                                                         cGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGa            >  1:427111/1‑150 (MQ=255)
                                                                                                                                           gTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGAtg          >  1:1507475/1‑150 (MQ=255)
                                                                                                                                              tgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGc       <  2:1583710/150‑1 (MQ=255)
                                                                                                                                              tgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGc       <  1:1664235/150‑1 (MQ=255)
                                                                                                                                              tgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGc       >  2:279057/1‑150 (MQ=255)
                                                                                                                                                ctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGcac     >  2:388867/1‑150 (MQ=255)
                                                                                                                                                   tgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGcacaac  <  1:381723/150‑1 (MQ=255)
                                                                                                                                                   tgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGcacaac  >  2:943318/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAAC  >  NC_000913/4132283‑4132579

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