Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,132,424 T→A 7.0% intergenic (+157/‑192) 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,4240TA 7.0% 499.2 / 15.8 169intergenic (+157/‑192)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.97e-01

Only 100 of 180 total aligned reads displayed.

ATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTG  >  NC_000913/4132276‑4132573
                                                                                                                                                    |                                                                                                                                                     
aTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGttt                                                                                                                                                      >  2:381723/1‑150 (MQ=255)
 tttAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttg                                                                                                                                                     >  2:1317796/1‑150 (MQ=255)
 tttAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttg                                                                                                                                                     <  2:618538/150‑1 (MQ=255)
       aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttgc                                                                                                                                               <  1:320322/150‑1 (MQ=255)
        atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                                                                                                                              <  2:182054/150‑1 (MQ=255)
          gCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATa                                                                                                                                            >  1:965407/1‑150 (MQ=255)
              ggTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCgg                                                                                                                                        <  2:696728/150‑1 (MQ=255)
              ggTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCgg                                                                                                                                        <  2:936939/150‑1 (MQ=255)
              ggTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGAGTTTACGACGAATGCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCgg                                                                                                                                        <  2:1277341/150‑1 (MQ=255)
               gTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGa                                                                                                                                       <  1:1504498/150‑1 (MQ=255)
               gTACTCATGTTTTCGGGGTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGa                                                                                                                                       <  2:1357396/150‑1 (MQ=255)
                tACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGac                                                                                                                                      <  2:1574927/150‑2 (MQ=255)
                  cTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                                                                                                                                    >  2:1245027/1‑150 (MQ=255)
                    cATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                                                                                                                                  >  1:1082348/1‑150 (MQ=255)
                    cATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                                                                                                                                  >  2:1667524/1‑150 (MQ=255)
                       gTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                                                                                                                               >  1:1213101/1‑150 (MQ=255)
                        ttttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttt                                                                                                                              >  2:838819/1‑150 (MQ=255)
                         tttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt                                                                                                                             <  2:1052127/150‑1 (MQ=255)
                          ttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                                                                                                                            <  2:549209/150‑1 (MQ=255)
                            cGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt                                                                                                                          >  1:1393626/1‑150 (MQ=255)
                              ggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc                                                                                                                        >  2:1457186/1‑150 (MQ=255)
                               gTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTATTATTCTTAcc                                                                                                                       >  2:266504/1‑150 (MQ=255)
                                  tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                    >  2:501368/1‑150 (MQ=255)
                                    tGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCt                                                                                                                  >  2:796403/1‑150 (MQ=255)
                                     ggTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg                                                                                                                 >  1:475202/1‑150 (MQ=255)
                                     ggTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg                                                                                                                 <  2:965407/150‑1 (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:306929/1‑150 (MQ=255)
                                            aaTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttt                                                                                                          <  1:1208967/150‑1 (MQ=255)
                                              tGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGATTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttttt                                                                                                        >  1:707635/1‑150 (MQ=255)
                                               gAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTc                                                                                                       >  2:119721/1‑150 (MQ=255)
                                                   tatatTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                                                                                                                            >  1:399437/1‑125 (MQ=255)
                                                      atTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTcc                                                                                                >  1:821639/1‑150 (MQ=255)
                                                         gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTTGTTTATATTTGCAGTCCGGTTTGATTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                             >  1:82978/1‑150 (MQ=255)
                                                         gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATCCCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                             >  2:1135261/1‑150 (MQ=255)
                                                         gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTGTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                             >  1:87476/1‑150 (MQ=255)
                                                            ttATCCTATGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGACCATGGGTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGAtt                                                                                          <  2:1113338/150‑1 (MQ=255)
                                                                cATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTc                                                                                      >  2:1664235/1‑150 (MQ=255)
                                                                cATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTc                                                                                      >  2:913300/1‑150 (MQ=255)
                                                                cATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTc                                                                                      >  1:1369992/1‑150 (MQ=255)
                                                                 aTAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTcc                                                                                     >  2:238453/1‑150 (MQ=255)
                                                                    ggATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCa                                                                                  >  2:829349/1‑150 (MQ=255)
                                                                    ggATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCa                                                                                  >  2:723488/1‑150 (MQ=255)
                                                                     gATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAc                                                                                 <  2:881555/150‑1 (MQ=255)
                                                                     gATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAc                                                                                 >  2:529616/1‑150 (MQ=255)
                                                                      aTTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACca                                                                               <  2:1503019/150‑1 (MQ=255)
                                                                           gCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACcatcat                                                                           >  2:1000653/1‑150 (MQ=255)
                                                                            ccGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACAATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCAtt                                                                          >  1:853759/1‑150 (MQ=255)
                                                                             cGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTg                                                                         <  2:475202/150‑1 (MQ=255)
                                                                              ggATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTgg                                                                        >  1:1036037/1‑150 (MQ=255)
                                                                                  tAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCAt                                                                    <  2:387452/150‑1 (MQ=255)
                                                                                   aaGAGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCAtt                                                                   <  2:472485/150‑1 (MQ=255)
                                                                                         ttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTctg                                                             >  1:778749/1‑148 (MQ=255)
                                                                                           acgacgAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGtt                                                           >  2:216876/1‑150 (MQ=255)
                                                                                             gacgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCTGTTGa                                                         >  1:1600431/1‑150 (MQ=255)
                                                                                              acgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGAc                                                        <  2:335119/150‑1 (MQ=255)
                                                                                              acgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGAc                                                        <  1:1587362/150‑1 (MQ=255)
                                                                                               cgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGGTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACg                                                       <  1:257304/150‑1 (MQ=255)
                                                                                                gaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGc                                                      >  2:1269441/1‑150 (MQ=255)
                                                                                                  aTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGccc                                                    >  1:947981/1‑150 (MQ=255)
                                                                                                   tCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAACATAATGAAGTTTTTTCATCTTCCCTGATTTTTACTCACCATCATTGGTCATTTTTCGGTTGACTCCCt                                                   >  1:75016/1‑150 (MQ=255)
                                                                                                     cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCtta                                                 >  1:477009/1‑149 (MQ=255)
                                                                                                      ggCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                 >  2:1528976/1‑149 (MQ=255)
                                                                                                       gCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCgg                                               <  2:1448079/150‑1 (MQ=255)
                                                                                                        cAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGc                                              <  1:1479018/150‑1 (MQ=255)
                                                                                                        cAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGc                                              <  1:1245027/150‑1 (MQ=255)
                                                                                                         aagaagCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCt                                             <  1:216876/150‑1 (MQ=255)
                                                                                                              gCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTcc                                        >  2:381302/1‑150 (MQ=255)
                                                                                                               cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTTGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                       >  2:863541/1‑150 (MQ=255)
                                                                                                               cAATAAGTACATGGTTAGTTTATATATGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACAATCCTGAAGTTTTTTCATCTTCCCCGATTTTTACTCACCATCATTGGTCATTTTTCGGTTTACTCCATTCGGCTTTTCCt                                       >  1:245765/1‑150 (MQ=255)
                                                                                                                aaTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCCTTTTTCGGTTGACGCCCTTCGGCTTTTCCtt                                      <  1:1135261/150‑1 (MQ=255)
                                                                                                                 aTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTc                                     <  1:834502/150‑1 (MQ=255)
                                                                                                                 aTAAGTACATGGGTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTc                                     <  2:1603985/150‑1 (MQ=255)
                                                                                                                    aGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATc                                  <  1:1003950/150‑1 (MQ=255)
                                                                                                                      tACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCtt                                >  2:1213028/1‑150 (MQ=255)
                                                                                                                       aCATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCttt                               >  1:1696177/1‑150 (MQ=255)
                                                                                                                             ttAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCt                         >  2:1413650/1‑150 (MQ=255)
                                                                                                                               aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                       >  1:1104757/1‑150 (MQ=255)
                                                                                                                               aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATATTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                       >  1:583208/1‑150 (MQ=255)
                                                                                                                               aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACAATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCGTCTTTACATCTgg                       >  1:1300444/1‑150 (MQ=255)
                                                                                                                                 tttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGAc                     <  2:412294/150‑1 (MQ=255)
                                                                                                                                 tttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACg                    <  1:723488/150‑1 (MQ=255)
                                                                                                                                  ttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACg                    <  2:657305/150‑1 (MQ=255)
                                                                                                                                  ttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACg                    <  1:913300/150‑1 (MQ=255)
                                                                                                                                   tataTTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGAAGt                   >  1:1259792/1‑150 (MQ=255)
                                                                                                                                    atatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                  <  1:46023/150‑1 (MQ=255)
                                                                                                                                     tatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACTCCCTTCGGCTTTTCCTACATCTTTACATCTGGACGTCt                 >  2:779744/1‑150 (MQ=255)
                                                                                                                                       tttGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCGTTTCCTTCATCTTTACATCTGGACGTCTaa               >  2:343597/1‑150 (MQ=255)
                                                                                                                                           aaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACgg           >  1:1573848/4‑150 (MQ=255)
                                                                                                                                             atCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGAt         >  1:504055/2‑150 (MQ=255)
                                                                                                                                              tctggTTTGCTTTGCATAACGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATa        <  1:989005/147‑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    <  2:707635/150‑1 (MQ=255)
                                                                                                                                                  gTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGAtg    >  1:1507475/1‑150 (MQ=255)
                                                                                                                                                   tttgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGAtgt   <  2:746658/150‑1 (MQ=255)
                                                                                                                                                    ttgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTAACCATCATTGGTCATTTTTCGGTTGATGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGAtttg  >  1:1707663/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTG  >  NC_000913/4132276‑4132573

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: