Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,331 C→T 8.0% intergenic (+42/‑80) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3310CT 8.0% 304.9 / 13.3 100intergenic (+42/‑80)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (35/57);  new base T (5/3);  total (40/60)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 101 total aligned reads displayed.

GTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGG  >  NC_000913/4612184‑4612478
                                                                                                                                                   |                                                                                                                                                    
gTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑g                                                                                                                                                    >  3:9733/1‑149 (MQ=255)
gTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑g                                                                                                                                                    <  4:784189/149‑1 (MQ=255)
 tttGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑g                                                                                                                                                    >  1:905605/1‑148 (MQ=255)
        gtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCt                                                                                                                                            <  3:328947/149‑1 (MQ=255)
         tCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑g                                                                                                                                                    <  3:948861/140‑1 (MQ=255)
          cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCtt                                                                                                                                           <  2:954998/148‑1 (MQ=255)
          cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCtt                                                                                                                                           <  2:463956/148‑1 (MQ=255)
           tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                          <  1:793471/148‑1 (MQ=255)
           tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑AGACTGGCttt                                                                                                                                          <  1:506142/148‑1 (MQ=255)
           tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCttt                                                                                                                                          <  2:445583/148‑1 (MQ=255)
            ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCtttg                                                                                                                                         <  1:953608/148‑2 (MQ=255)
             gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGc                                                                                                                                       <  4:373671/149‑1 (MQ=255)
              ccGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGc                                                                                                                                       >  1:628371/1‑148 (MQ=255)
              ccGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGc                                                                                                                                       >  1:850666/1‑148 (MQ=255)
                ggTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCgg                                                                                                                                     <  1:603339/148‑1 (MQ=255)
                  tACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGttt                                                                                                                                  <  4:620944/149‑1 (MQ=255)
                      gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTTAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCtt                                                                                                                                           <  2:102935/136‑1 (MQ=255)
                      gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTTAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCtt                                                                                                                                           >  1:102935/1‑136 (MQ=255)
                      gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTCAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGc                                                                                                                              <  4:585286/149‑1 (MQ=255)
                         gCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaa                                                                                                                           <  3:225252/149‑1 (MQ=255)
                         gCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaa                                                                                                                           <  3:728163/149‑1 (MQ=255)
                            ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaataa                                                                                                                        >  3:858067/1‑149 (MQ=255)
                             gcgcAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaataa                                                                                                                        <  2:628371/148‑1 (MQ=255)
                                cAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGTGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAAtaataa                                                                                                                     <  2:789626/148‑1 (MQ=255)
                                cAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAAtaataa                                                                                                                     <  2:669909/148‑1 (MQ=255)
                                cAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                                                                                                    <  4:62432/149‑1 (MQ=255)
                                      aaaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCttt                                                                                                                                          <  4:874682/121‑1 (MQ=255)
                                      aaaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCttt                                                                                                                                          >  3:874682/1‑121 (MQ=255)
                                       aaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGtt                                                                                                             >  4:822672/1‑149 (MQ=255)
                                       aaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGtt                                                                                                             >  4:652874/1‑149 (MQ=255)
                                          tGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACt                                                                                                          <  3:764265/149‑1 (MQ=255)
                                           gTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                         >  4:346148/1‑149 (MQ=255)
                                              tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                      >  4:814959/1‑149 (MQ=255)
                                               ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                     <  4:400051/149‑1 (MQ=255)
                                                  tcgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                                                   <  1:219554/148‑1 (MQ=255)
                                                   cgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGCTACTTTGTCAtt                                                                                                  >  2:948787/1‑148 (MQ=255)
                                                    gtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCt                                                                                                <  3:814959/149‑1 (MQ=255)
                                                       gggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAGATTGCGTTACTTTGTCATTCttt                                                                                              <  2:290286/148‑1 (MQ=255)
                                                       gggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCttt                                                                                              <  1:281209/148‑1 (MQ=255)
                                                        ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                            <  3:822672/149‑1 (MQ=255)
                                                        ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                             >  2:634213/1‑148 (MQ=255)
                                                        ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                             <  1:507015/148‑1 (MQ=255)
                                                         gATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcac                                                                                           >  3:350410/1‑149 (MQ=255)
                                                           ttattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                         <  3:697927/149‑1 (MQ=255)
                                                                   ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc                                                                                  <  2:905605/148‑1 (MQ=255)
                                                                   ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                 <  3:593353/149‑1 (MQ=255)
                                                                   ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                 <  4:152379/149‑1 (MQ=255)
                                                                       tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAACTTGCGTTACTTTATCATCCTTTCACAACGAAAACCGGaa                                                                             <  4:816359/149‑1 (MQ=255)
                                                                       tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGa                                                                              >  1:190799/1‑148 (MQ=255)
                                                                              tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                      >  3:546803/1‑149 (MQ=255)
                                                                              tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                       <  1:634213/148‑1 (MQ=255)
                                                                                  ttCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  <  3:37422/149‑1 (MQ=255)
                                                                                        ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAg                                                             <  2:568826/148‑1 (MQ=255)
                                                                                         gCCGAAAACGACCCTAGATTTCAGACGATTAAGAAAATGTTATCCAAAGCCAGTCCAGC‑GGACTGGTTTTTGCGGTTTTAGCGAATAATAAATTGCGTGACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                           >  3:800663/1‑149 (MQ=255)
                                                                                           cGAAAACGACCCTAGATTTCAGACCGTTAAGAAAATCTTATCCTAAGCCAGTCCAGC‑GGTCTGTCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACACAAACGGGTAGGCAGAGGTGGGGCAGCaaa                                                         <  4:800663/149‑1 (MQ=255)
                                                                                           cGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGCCTGGCTTTTTCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGa                                                          <  2:4480/148‑1 (MQ=255)
                                                                                           cGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGaa                                                         <  4:564759/149‑1 (MQ=255)
                                                                                             aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAAt                                                        >  2:441054/1‑148 (MQ=255)
                                                                                                aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                    <  4:858067/149‑1 (MQ=255)
                                                                                                   aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGt                                                  <  1:657140/148‑1 (MQ=255)
                                                                                                      cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                              >  3:106648/1‑149 (MQ=255)
                                                                                                      cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                              >  3:840948/1‑149 (MQ=255)
                                                                                                      cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                              >  3:342614/1‑149 (MQ=255)
                                                                                                          aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                          <  4:402998/149‑1 (MQ=255)
                                                                                                                gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAat                                     >  1:449078/1‑148 (MQ=255)
                                                                                                                  ccATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATAtt                                  <  3:894981/149‑1 (MQ=255)
                                                                                                                    aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                                 <  2:778201/148‑1 (MQ=255)
                                                                                                                        aGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  <  2:339875/111‑1 (MQ=255)
                                                                                                                        aGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  >  1:339875/1‑111 (MQ=255)
                                                                                                                          aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                   >  1:820695/1‑108 (MQ=255)
                                                                                                                          aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                   <  2:820695/108‑1 (MQ=255)
                                                                                                                          aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGaa                                                         <  3:410700/118‑1 (MQ=255)
                                                                                                                          aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGaa                                                         >  4:410700/1‑118 (MQ=255)
                                                                                                                          aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt                           <  2:448582/148‑1 (MQ=255)
                                                                                                                          aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt                           >  1:448582/1‑148 (MQ=255)
                                                                                                                               tttATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCg                      <  2:190799/148‑1 (MQ=255)
                                                                                                                                ttATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                          <  4:439457/127‑1 (MQ=255)
                                                                                                                                ttATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                          >  3:439457/1‑127 (MQ=255)
                                                                                                                                 tATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTg                                                   <  2:768989/117‑1 (MQ=255)
                                                                                                                                 tATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTg                                                   >  1:768989/1‑117 (MQ=255)
                                                                                                                                  aTCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTa                   <  2:410291/148‑1 (MQ=255)
                                                                                                                                   tCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg                 >  3:138496/1‑149 (MQ=255)
                                                                                                                                    ccAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                               <  4:497305/134‑1 (MQ=255)
                                                                                                                                    ccAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                               >  3:497305/1‑134 (MQ=255)
                                                                                                                                    caaaaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg                >  4:237684/3‑149 (MQ=255)
                                                                                                                                     cAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGc               >  3:169335/1‑149 (MQ=255)
                                                                                                                                      aaaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  3:307508/1‑149 (MQ=255)
                                                                                                                                      aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCCGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  3:239421/1‑149 (MQ=255)
                                                                                                                                       aaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg             >  3:307965/1‑149 (MQ=255)
                                                                                                                                       aaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                      <  1:864195/108‑1 (MQ=255)
                                                                                                                                       aaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                      >  2:864195/1‑108 (MQ=255)
                                                                                                                                        aGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg             >  1:318111/1‑148 (MQ=255)
                                                                                                                                          ccAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt           >  2:949511/1‑148 (MQ=255)
                                                                                                                                             gTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAg       <  3:164289/149‑1 (MQ=255)
                                                                                                                                              tCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGc      <  3:909963/149‑1 (MQ=255)
                                                                                                                                                cAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGcc     <  2:115089/148‑1 (MQ=255)
                                                                                                                                                 aGC‑GGACTGGCTTTTGCGGTTTTAGCGAAtaataa                                                                                                                     <  1:731269/35‑1 (MQ=255)
                                                                                                                                                 aGC‑GGACTGGCTTTTGCGGTTTTAGCGAAtaataa                                                                                                                     >  2:731269/1‑35 (MQ=255)
                                                                                                                                                 aGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCt    <  2:850666/148‑1 (MQ=255)
                                                                                                                                                   c‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTgg  <  1:441054/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGG  >  NC_000913/4612184‑4612478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑