Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,062,354 T→A 6.2% intergenic (+19/+37) asnV → / ← ldtA tRNA‑Asn/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,062,3540TA 6.2% 391.1 / 11.9 113intergenic (+19/+37)asnV/ldtAtRNA‑Asn/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (50/56);  new base A (3/4);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 115 total aligned reads displayed.

CCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGC  >  NC_000913/2062210‑2062498
                                                                                                                                                |                                                                                                                                                
ccTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTa                                                                                                                                              >  3:181821/1‑149 (MQ=255)
    tGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAg                                                                                                                                          <  4:613168/149‑1 (MQ=255)
          cACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttt                                                                                                                                    <  4:572816/149‑1 (MQ=255)
           aCCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCtttt                                                                                                                                    <  2:728727/148‑1 (MQ=255)
           aCCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCttttt                                                                                                                                   <  4:53151/149‑1 (MQ=255)
            ccATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTg                                                                                                                                  >  3:761976/1‑149 (MQ=255)
             cATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCATAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc                                                                                                                                 >  3:609329/1‑149 (MQ=255)
             cATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCttttttg                                                                                                                                 <  4:933568/149‑3 (MQ=17)
              aTGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCGGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCt                                                                                                                                <  4:670115/149‑1 (MQ=255)
              aTGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc                                                                                                                                 <  2:106646/148‑1 (MQ=255)
               tGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCtt                                                                                                                               <  4:916711/149‑1 (MQ=255)
               tGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCt                                                                                                                                <  2:514620/148‑1 (MQ=255)
               tGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCt                                                                                                                                >  1:899680/1‑148 (MQ=255)
                  aCTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCt                                                                                                                            >  4:177094/1‑149 (MQ=255)
                    tGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTaa                                                                                                                          >  4:639546/1‑149 (MQ=255)
                     gCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCttt                                                                                                                                                 >  2:374164/1‑125 (MQ=255)
                      ccATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAAtt                                                                                                                        >  4:87542/1‑149 (MQ=255)
                        aTTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTAc                                                                                                                      >  3:921434/1‑149 (MQ=255)
                         ttAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTAcc                                                                                                                     >  4:136472/1‑149 (MQ=255)
                         ttAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTAc                                                                                                                      >  1:265580/1‑148 (MQ=255)
                          tAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCa                                                                                                                    <  4:712524/149‑1 (MQ=255)
                           aaTATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCa                                                                                                                    <  1:66981/148‑1 (MQ=255)
                            atatGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAACCCTGCTTTCTa                                                                                                                                              >  4:467169/1‑121 (MQ=255)
                            atatGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAACCCTGCTTTCTa                                                                                                                                              <  3:467169/121‑1 (MQ=255)
                              atGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGc                                                                                                                <  4:232866/149‑1 (MQ=255)
                                  gcCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt                                                                                                                                  >  4:122347/1‑126 (MQ=14)
                                  gcCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt                                                                                                                                  <  3:122347/127‑2 (MQ=14)
                                   ccccGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCtt                                                                                                            >  2:647235/1‑148 (MQ=255)
                                    cccGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTaa                                                                                                          >  3:104596/1‑149 (MQ=255)
                                      cGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTaaa                                                                                                         <  1:291902/148‑1 (MQ=255)
                                         tCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACa                                                                                                      >  2:633977/1‑148 (MQ=255)
                                          cacaCGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTCAAAACATc                                                                                                    <  4:84899/149‑1 (MQ=255)
                                           acacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCt                                                                                                   >  4:712317/1‑149 (MQ=255)
                                              cGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTc                                                                                                <  3:348361/149‑1 (MQ=255)
                                               gATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCt                                                                                               <  4:93466/149‑1 (MQ=255)
                                                     tctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAAcc                                                                                         >  3:173026/1‑149 (MQ=255)
                                                      ctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCa                                                                                        <  4:384800/149‑1 (MQ=255)
                                                        gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACa                                                                                                      <  1:423747/133‑1 (MQ=255)
                                                        gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACa                                                                                                      >  2:423747/1‑133 (MQ=255)
                                                           gTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACt                                                                                   >  3:793758/1‑149 (MQ=255)
                                                             tCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTaa                                                                                 <  4:921434/149‑1 (MQ=255)
                                                             tCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTaa                                                                                 <  3:620091/149‑1 (MQ=255)
                                                             tCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTa                                                                                  <  1:181862/148‑1 (MQ=255)
                                                              cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAAt                                                                                <  3:55098/149‑1 (MQ=255)
                                                                  cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGc                                                                            <  4:604451/149‑1 (MQ=255)
                                                                  cGGGAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGc                                                                            >  4:500156/1‑149 (MQ=255)
                                                                   ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGca                                                                           <  4:30209/149‑1 (MQ=255)
                                                                     tAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGcac                                                                          >  1:270187/1‑148 (MQ=255)
                                                                     tAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGcac                                                                          <  1:495246/148‑1 (MQ=255)
                                                                      aGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGcaca                                                                         <  2:955495/148‑1 (MQ=255)
                                                                        taCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAgg                                                                       >  1:226675/2‑148 (MQ=255)
                                                                          cggcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGa                                                                                            >  1:541730/1‑125 (MQ=255)
                                                                          cggcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGa                                                                                            <  2:541730/125‑1 (MQ=255)
                                                                                 cTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGa                                                              <  2:969294/148‑1 (MQ=255)
                                                                                  tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATc                                                            >  4:23721/1‑149 (MQ=255)
                                                                                  tGTTAATCCGTATGTCACTGGTTCGAGGCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACGGATTTGCCCAGGCATTCCCGATc                                                            >  4:635776/1‑149 (MQ=255)
                                                                                    ttAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCg                                                           >  2:717905/1‑148 (MQ=255)
                                                                                     tAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGa                                                          >  1:60181/1‑148 (MQ=255)
                                                                                      aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACg                                                        >  3:294849/1‑149 (MQ=255)
                                                                                       aTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACg                                                        <  1:680641/148‑1 (MQ=255)
                                                                                            tATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCaa                                                  >  4:578637/1‑149 (MQ=255)
                                                                                             aTGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCaa                                                  <  1:91190/148‑1 (MQ=255)
                                                                                                tCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgca                                              <  3:136472/149‑1 (MQ=255)
                                                                                                 cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc                                                                                                                                 <  1:593447/65‑1 (MQ=255)
                                                                                                 cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc                                                                                                                                 >  2:593447/1‑65 (MQ=255)
                                                                                                  aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcag                                             >  2:128431/1‑148 (MQ=255)
                                                                                                   cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcagc                                            >  1:796964/1‑148 (MQ=255)
                                                                                                    tGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAt                                          >  3:296945/1‑149 (MQ=255)
                                                                                                      gTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAtt                                         <  1:535358/148‑1 (MQ=255)
                                                                                                       ttCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTg                                       <  4:891783/149‑1 (MQ=255)
                                                                                                        tCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTgc                                      <  4:756351/149‑1 (MQ=255)
                                                                                                            gTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGcgcg                                   <  2:112510/148‑1 (MQ=255)
                                                                                                             tCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGa                                  >  1:131247/1‑148 (MQ=255)
                                                                                                              ccAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGAtt                                <  4:466688/149‑1 (MQ=255)
                                                                                                              ccAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGAtt                                <  3:431367/149‑1 (MQ=255)
                                                                                                                aGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTa                              >  3:100428/1‑149 (MQ=255)
                                                                                                                 gTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTa                              <  1:505263/148‑1 (MQ=255)
                                                                                                                 gTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACATTGCAACGCAGCATTTGCGCGATTTAc                             <  3:639546/149‑1 (MQ=255)
                                                                                                                  tCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa                            >  4:162754/1‑149 (MQ=255)
                                                                                                                   cAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa                            <  1:703811/148‑1 (MQ=255)
                                                                                                                      aggagCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCaa                        <  3:254224/149‑1 (MQ=255)
                                                                                                                       ggagCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCaa                        <  2:596623/148‑1 (MQ=255)
                                                                                                                            cAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAActtctt                  <  3:487324/149‑1 (MQ=255)
                                                                                                                            cAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAActtct                   <  1:499039/148‑1 (MQ=255)
                                                                                                                             aaaTTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAAGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGAGTTACATCAACTTCTTg                 >  4:570590/1‑149 (MQ=255)
                                                                                                                              aaTTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTg                 <  2:491069/148‑1 (MQ=255)
                                                                                                                              aaTTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGc                <  3:162754/149‑1 (MQ=255)
                                                                                                                                   aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGtt           <  4:294849/149‑1 (MQ=255)
                                                                                                                                   aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGtt           >  3:770878/1‑149 (MQ=255)
                                                                                                                                   aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGt            <  2:347025/148‑1 (MQ=255)
                                                                                                                                   aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGt            <  1:688101/148‑1 (MQ=255)
                                                                                                                                    aaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTg          <  4:175590/149‑1 (MQ=255)
                                                                                                                                    aaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTg          >  4:453371/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTg          >  2:440352/1‑148 (MQ=255)
                                                                                                                                      aaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgc                                               <  2:686973/110‑1 (MQ=255)
                                                                                                                                      aaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgc                                               >  1:686973/1‑110 (MQ=255)
                                                                                                                                      aaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGa         >  1:322804/1‑148 (MQ=255)
                                                                                                                                       aaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATa       <  3:647106/149‑1 (MQ=255)
                                                                                                                                            tGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACg   >  2:412222/1‑148 (MQ=255)
                                                                                                                                             gCTTTCCAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGc  >  1:961396/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGC  >  NC_000913/2062210‑2062498

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