Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,132,429 T→A 7.3% intergenic (+162/‑187) 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,4290TA 7.3% 356.5 / 14.2 123intergenic (+162/‑187)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (46/68);  new base A (5/4);  total (51/72)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

Only 100 of 133 total aligned reads displayed.

AATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAA  >  NC_000913/4132280‑4132578
                                                                                                                                                     |                                                                                                                                                     
aataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctt                                                                                                                                                       >  2:554341/1‑150 (MQ=255)
aataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctt                                                                                                                                                       >  1:907842/1‑150 (MQ=255)
  taatGCCCTGTACTCATGTTTTAGGGTTTATGTTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttg                                                                                                                                                     <  1:340355/150‑1 (MQ=255)
    atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                                                                                                                                   <  2:506634/150‑1 (MQ=255)
      gCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATa                                                                                                                                                 >  2:309879/1‑150 (MQ=255)
         cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGTCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                                                                                                                              <  1:505686/150‑1 (MQ=255)
         cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                                                                                                                              >  2:1628525/1‑150 (MQ=255)
         cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                                                                                                                              >  2:1232045/1‑150 (MQ=255)
            tACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAt                                                                                                                                           <  1:1604350/150‑1 (MQ=255)
            tACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGac                                                                                                                                           <  1:1347079/150‑2 (MQ=255)
             aCTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGact                                                                                                                                          <  2:1826887/150‑3 (MQ=255)
              cTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                                                                                                                                         >  2:1028072/1‑150 (MQ=255)
               tCATGTTTTCGTGTTTATGGTTTCTAATGAAATATATTGATTTATCATAGGTTTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                                                                                                                                        <  1:195087/150‑1 (MQ=255)
               tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAATAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                                                                                                                                        >  1:1599612/1‑150 (MQ=255)
               tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                                                                                                                                        <  1:1515152/150‑1 (MQ=255)
                 aTGTTTTCTGGTTTATGGTTTCTAATGAAATATATTGATTTATCCTAGGATTAGGGCGGATTTAGCGTTTACGACGTATCCGGCAAGAAGCACTAAGTACATGGTTTGCTTATATTTGCTGTCCGGTTTGCTTTGCATACCGGTttttct                                                                                                                                      <  2:1654237/150‑1 (MQ=255)
                 aTGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                                                                                                                                      <  2:1697849/150‑1 (MQ=255)
                 aTGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                                                                                                                                      <  2:1555050/150‑1 (MQ=255)
                      ttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                                                                                                                                 <  1:1333196/150‑1 (MQ=255)
                       tCGGGTTTATGGTTTCTTATGAAATATATTGTTTTATGATTGGATTAGGCCGGATTAATCGTTTACGACGAATCCGGCAAGTAGCAATAAGTACATGGTTAGTTTTTATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttct                                                                                                                                <  2:1175938/150‑1 (MQ=255)
                        cGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTCCGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt                                                                                                                               <  1:491261/150‑1 (MQ=255)
                         gggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa                                                                                                                              <  2:570698/150‑1 (MQ=255)
                          ggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGATGCAATAAGTCCATGGTTAGTTGATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTaa                                                                                                                             >  2:1060465/1‑149 (MQ=255)
                              tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                         >  1:1513151/1‑150 (MQ=255)
                              tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                         >  1:584879/1‑150 (MQ=255)
                                 ggTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACAACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg                                                                                                                      >  2:757831/1‑150 (MQ=255)
                                       tAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttt                                                                                                                >  2:36077/1‑150 (MQ=255)
                                       tAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttt                                                                                                                >  2:750818/1‑150 (MQ=255)
                                         aTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttttt                                                                                                              <  1:658961/150‑1 (MQ=255)
                                           gAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTc                                                                                                            >  1:737444/1‑150 (MQ=255)
                                           gAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTc                                                                                                            >  2:1701045/1‑150 (MQ=255)
                                                  atTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTcc                                                                                                     >  2:85497/1‑150 (MQ=255)
                                                  atTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTcc                                                                                                     >  1:1494316/1‑150 (MQ=255)
                                                   ttGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTccc                                                                                                    <  1:548050/150‑1 (MQ=255)
                                                    tGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCt                                                                                                   <  2:1624899/150‑1 (MQ=255)
                                                     gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                                  >  1:108918/1‑150 (MQ=255)
                                                      aaTTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGa                                                                                                 <  2:1172396/150‑1 (MQ=255)
                                                           tCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGAttttt                                                                                            <  2:900410/150‑1 (MQ=255)
                                                             aTAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTcc                                                                                          >  1:126492/1‑150 (MQ=255)
                                                                 gATTAGGCCGGATTAAGAGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAc                                                                                      <  2:418259/150‑1 (MQ=255)
                                                                        ccGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCAtt                                                                               >  1:482510/1‑150 (MQ=255)
                                                                         cGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTg                                                                              <  2:1277068/150‑1 (MQ=255)
                                                                           gATTTAGCGTTTACGACGAATCCAGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGATTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGt                                                                            <  1:1379233/150‑1 (MQ=255)
                                                                            aTTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTc                                                                           <  1:66284/150‑1 (MQ=255)
                                                                                   gTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTc                                                                    <  2:1550332/150‑1 (MQ=255)
                                                                                   gTTTACAACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTc                                                                    <  1:757831/150‑1 (MQ=255)
                                                                                    tttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTTGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCg                                                                   >  2:515700/1‑150 (MQ=255)
                                                                                    tttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCg                                                                   >  1:898182/1‑150 (MQ=255)
                                                                                    tttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCg                                                                   <  1:353803/150‑1 (MQ=255)
                                                                                     ttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCgg                                                                  >  2:1435974/1‑150 (MQ=255)
                                                                                     ttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCgg                                                                  <  2:676723/150‑1 (MQ=255)
                                                                                      tACGACGTATACGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGt                                                                 <  2:484845/150‑1 (MQ=255)
                                                                                      tACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGt                                                                 <  2:850730/150‑1 (MQ=255)
                                                                                       acgacgAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGtt                                                                <  1:1757309/150‑1 (MQ=255)
                                                                                           cgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACg                                                            <  1:889790/150‑1 (MQ=255)
                                                                                             aaTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCAACTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGcc                                                          >  2:764254/1‑150 (MQ=255)
                                                                                              aTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGccc                                                         <  1:1535197/150‑1 (MQ=255)
                                                                                                ccGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCtt                                                       >  2:902522/1‑150 (MQ=255)
                                                                                                 cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                      >  1:1281554/1‑150 (MQ=255)
                                                                                                 cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                      <  2:835081/150‑1 (MQ=255)
                                                                                                 cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTGTTCGGTTGACGCCCTTc                                                      >  2:1242389/1‑150 (MQ=255)
                                                                                                 cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTATCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                      >  2:1116634/1‑150 (MQ=255)
                                                                                                         agCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGTCTTTTc                                              >  1:1441077/1‑150 (MQ=255)
                                                                                                          gCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTcc                                             <  1:1285554/150‑1 (MQ=255)
                                                                                                           cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                            >  2:1107118/1‑150 (MQ=255)
                                                                                                           cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                            <  1:1339746/150‑1 (MQ=255)
                                                                                                           cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                            <  2:584879/150‑1 (MQ=255)
                                                                                                             aTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTc                                          <  1:1028072/150‑1 (MQ=255)
                                                                                                              tAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCa                                         <  1:1214033/150‑1 (MQ=255)
                                                                                                              tAAGTACATGGTTAGTTTATATTTGCAGTCCGGATTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCa                                         <  1:241270/150‑1 (MQ=255)
                                                                                                               aaGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCAt                                        <  2:1826697/150‑1 (MQ=255)
                                                                                                                  tACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCtt                                     <  1:610092/150‑1 (MQ=255)
                                                                                                                  tACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCtt                                     <  2:179063/150‑1 (MQ=255)
                                                                                                                     aTGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTAc                                  <  1:1232045/150‑1 (MQ=255)
                                                                                                                        gTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATc                               <  2:1322509/150‑1 (MQ=255)
                                                                                                                          tAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTg                             <  1:674858/150‑1 (MQ=255)
                                                                                                                          tAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTATTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTg                             <  1:252017/150‑1 (MQ=255)
                                                                                                                           aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                            <  1:1604829/150‑1 (MQ=255)
                                                                                                                            gTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACAATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGa                           >  1:1046299/1‑150 (MQ=255)
                                                                                                                            gTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAGTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATATTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCGTCATCTTTACATCTGGa                           >  1:55735/1‑150 (MQ=255)
                                                                                                                                atatTTTCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATACTGAAGTTTTATCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                       <  1:668860/150‑1 (MQ=255)
                                                                                                                                atatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                       <  1:1201456/150‑1 (MQ=255)
                                                                                                                                   aagaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTaa                    >  1:1637193/8‑150 (MQ=255)
                                                                                                                                    agaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTaaa                   >  2:690399/7‑150 (MQ=255)
                                                                                                                                    agaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTCTCCTCAACATCATTGGTCATTTTTCGGTTGACGCCCTTCGGATTTTCCTTCATCTTTACATCTGGACGTCTaaa                   >  2:1135570/7‑150 (MQ=255)
                                                                                                                                      gCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACg                 <  1:452478/150‑1 (MQ=255)
                                                                                                                                       cAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACgg                <  2:1094633/150‑1 (MQ=255)
                                                                                                                                       aaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACgg                >  1:1789020/4‑150 (MQ=255)
                                                                                                                                        aGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGa               <  1:5236/150‑1 (MQ=255)
                                                                                                                                         atCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGAt              >  2:1309091/2‑150 (MQ=255)
                                                                                                                                           ccGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGTTTTTTCCTCACCATCATTGGTCATTTTTAGTTTGACGACCTTAGGCTTTTCCTTCATCTTTACATATGTACTTATAAACGGATaa            >  1:361315/1‑149 (MQ=255)
                                                                                                                                           ccGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAg            >  2:569557/1‑150 (MQ=255)
                                                                                                                                            cGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGa           <  1:309879/150‑1 (MQ=255)
                                                                                                                                            cGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCATTCGGATTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGa           >  1:213469/1‑150 (MQ=255)
                                                                                                                                              gTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGAtg         <  2:963251/150‑1 (MQ=255)
                                                                                                                                                 tgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGc      <  2:773014/150‑1 (MQ=255)
                                                                                                                                                  gctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGca     <  1:999977/150‑1 (MQ=255)
                                                                                                                                                  gctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAAAGGATAGATGTGca     >  1:1144275/1‑150 (MQ=255)
                                                                                                                                                    tttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCAGGACGTCTAAACGGATAGATGTGcaca   <  2:1513151/150‑1 (MQ=255)
                                                                                                                                                     ttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGcacaa  >  1:312165/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAA  >  NC_000913/4132280‑4132578

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: ‑