Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,715 Δ1 bp 6.6% intergenic (+22/‑145) sppA → / → ansA protease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,850,7150A. 6.6% 420.1 / 22.3 122intergenic (+22/‑145)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base A (52/61);  new base . (4/4);  total (57/65)
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.69e-01

Only 100 of 123 total aligned reads displayed.

GGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATG  >  NC_000913/1850567‑1850862
                                                                                                                                                    |                                                                                                                                                    
ggCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACa                                                                                                                                                      >  1:378271/1‑149 (MQ=255)
    aTGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGGCCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtc                                                                                                                                                  <  1:1159989/149‑1 (MQ=255)
        tACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtc‑gtcg                                                                                                                                             <  4:473788/149‑1 (MQ=255)
          ccTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGGCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGc                                                                                                                                           <  3:803776/149‑1 (MQ=255)
           cTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGcc                                                                                                                                          <  4:691247/149‑1 (MQ=255)
           cTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGcc                                                                                                                                          <  4:656479/149‑1 (MQ=255)
           cTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGAACTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGcc                                                                                                                                          <  3:135435/149‑1 (MQ=255)
               aCCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACat                                                                                                                                      >  1:736640/1‑149 (MQ=255)
                 cGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGcc                                                                                                                                          >  1:1074724/1‑143 (MQ=255)
                 cGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGcc                                                                                                                                          <  2:1074724/143‑1 (MQ=255)
                  gCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACattat                                                                                                                                   <  1:707443/149‑1 (MQ=255)
                   cTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGGGGCCGACAGATC‑GGCGGCCACattatg                                                                                                                                  <  2:1558293/149‑2 (MQ=255)
                     ggCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttt                                                                                                                                <  4:1404514/149‑1 (MQ=255)
                      gCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAttttt                                                                                                                               <  3:898508/149‑1 (MQ=255)
                      gCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAttttt                                                                                                                               >  2:405121/1‑149 (MQ=255)
                        cTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTa                                                                                                                             <  1:1079321/149‑1 (MQ=255)
                        cTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTa                                                                                                                             >  1:1028599/1‑149 (MQ=255)
                           ggTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGt                                                                                                                          <  3:168860/149‑1 (MQ=255)
                              aGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcga                                                                                                                       <  2:1390744/149‑1 (MQ=255)
                               gCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                                                                      <  1:405121/149‑1 (MQ=255)
                                   ctACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAAt                                                                                                                  <  2:1163910/149‑1 (MQ=255)
                                   ctACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAAt                                                                                                                  <  2:168117/149‑1 (MQ=255)
                                      ccGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCct                                                                                                               <  2:445671/149‑1 (MQ=255)
                                        gTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACTTCGACGAATCctct                                                                                                             >  3:597473/1‑149 (MQ=255)
                                         ttAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCtt                                                                                                            >  3:1086142/1‑149 (MQ=255)
                                         ttAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCtt                                                                                                            >  4:1431697/1‑149 (MQ=255)
                                         ttAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCtt                                                                                                            >  1:1553194/1‑149 (MQ=255)
                                            aaagtgaaagtgACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTccc                                                                                                         <  2:734499/149‑1 (MQ=255)
                                             aagtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                                                                              <  1:892064/143‑1 (MQ=255)
                                              agtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGc                                                                                                       <  2:1047697/149‑1 (MQ=255)
                                                tgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTg                                                                                                     >  2:921624/1‑149 (MQ=255)
                                                tgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTg                                                                                                     >  1:996572/1‑149 (MQ=255)
                                                tgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTg                                                                                                     >  2:354088/1‑149 (MQ=255)
                                                 gaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGt                                                                                                    <  3:475732/149‑1 (MQ=255)
                                                    agtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                                                                 <  2:789547/149‑1 (MQ=255)
                                                      tgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                                                                               <  4:1240146/149‑1 (MQ=255)
                                                       gaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                                                                              <  2:736640/149‑1 (MQ=255)
                                                           aGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCat                                                                                          >  3:1163997/1‑149 (MQ=255)
                                                                 ccGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTcc                                                                                    <  4:597473/149‑1 (MQ=255)
                                                                   gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCtt                                                                                  <  2:996572/149‑1 (MQ=255)
                                                                     gcttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTAt                                                                                >  3:245766/3‑149 (MQ=255)
                                                                       ttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATcc                                                                              <  2:100711/149‑1 (MQ=255)
                                                                         aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                                                                            >  1:353487/1‑149 (MQ=255)
                                                                            gACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATc                                                                                                                 >  2:996449/1‑109 (MQ=255)
                                                                            gACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGtt                                                                                                   >  2:1538964/1‑123 (MQ=255)
                                                                              cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGt                                                                       >  4:129070/1‑149 (MQ=255)
                                                                               ccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTa                                                                      >  1:787776/1‑149 (MQ=255)
                                                                                cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                                                                     >  4:611309/1‑149 (MQ=255)
                                                                                 acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                                                                    >  2:1504594/1‑149 (MQ=255)
                                                                                 acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                                                                    <  2:1204788/149‑1 (MQ=255)
                                                                                 acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                                                                    <  4:33504/149‑1 (MQ=255)
                                                                                    aaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAgg                                                                 <  3:915498/149‑1 (MQ=255)
                                                                                     aaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGt                                                                <  3:1431697/149‑1 (MQ=255)
                                                                                       ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATc                                                                                                                 <  1:939418/98‑1 (MQ=255)
                                                                                       ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATc                                                                                                                 >  2:939418/1‑98 (MQ=255)
                                                                                             tGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCa                                                                           >  2:1100920/1‑130 (MQ=255)
                                                                                                gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCt                                                     <  4:1086142/149‑1 (MQ=255)
                                                                                                gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCt                                                     >  2:988065/1‑149 (MQ=255)
                                                                                                 ttttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTg                                                    >  2:1397212/1‑149 (MQ=255)
                                                                                                     gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTc                                                >  1:28622/1‑149 (MQ=255)
                                                                                                     gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTc                                                <  3:129070/149‑1 (MQ=255)
                                                                                                      ccTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCACTCCCCCTGCCTCa                                               >  2:1153945/1‑149 (MQ=255)
                                                                                                         gACTTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATA‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGtt                                                                                                   >  2:712077/1‑94 (MQ=255)
                                                                                                         gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGtt                                                                                                   <  1:712077/94‑1 (MQ=255)
                                                                                                            cTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGtata                                         >  1:113244/1‑149 (MQ=255)
                                                                                                                 ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACtt                                    >  4:914127/1‑149 (MQ=255)
                                                                                                                  cAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACttt                                   >  1:1238989/1‑149 (MQ=255)
                                                                                                                   aaCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACtttt                                  <  1:354088/149‑1 (MQ=255)
                                                                                                                     cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGc                                <  1:1504594/149‑1 (MQ=255)
                                                                                                                     cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGc                                >  1:1017937/1‑149 (MQ=255)
                                                                                                                      aTGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGct                               >  2:779906/1‑149 (MQ=255)
                                                                                                                           ttAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTc                          >  4:853482/1‑149 (MQ=255)
                                                                                                                           ttAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTc                          <  3:956975/149‑1 (MQ=255)
                                                                                                                              aGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTGTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGat                       <  4:3753/149‑1 (MQ=255)
                                                                                                                                tCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTGTCAGGTGCGCTCCCCCTGCCTCACGtatata                                       >  1:1109964/1‑131 (MQ=255)
                                                                                                                                tCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTGTCAGGTGCGCTCCCCCTGCCTCACGtatata                                       <  2:1109964/131‑1 (MQ=255)
                                                                                                                                tCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGatat                     >  1:211848/1‑149 (MQ=255)
                                                                                                                                  ttGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                   <  1:1197526/149‑1 (MQ=255)
                                                                                                                                  ttGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                   >  1:1234806/1‑149 (MQ=255)
                                                                                                                                  ttGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                   >  4:1468619/1‑149 (MQ=255)
                                                                                                                                  tcgTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                   >  2:83177/3‑149 (MQ=255)
                                                                                                                                   tGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATTAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcat                  >  1:202340/1‑149 (MQ=255)
                                                                                                                                        tGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCttt                           >  2:1002584/1‑135 (MQ=255)
                                                                                                                                        tGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCttt                           <  1:1002584/135‑1 (MQ=255)
                                                                                                                                        tGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCata             >  4:821208/1‑149 (MQ=255)
                                                                                                                                          aGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatc           <  3:1145350/149‑1 (MQ=255)
                                                                                                                                           gTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGt                                                                <  4:257206/95‑1 (MQ=255)
                                                                                                                                           gTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGt                                                                >  3:257206/1‑95 (MQ=255)
                                                                                                                                           gTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatca          >  2:1235873/1‑149 (MQ=255)
                                                                                                                                            tGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattc                >  1:1520631/1‑142 (MQ=255)
                                                                                                                                            tGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattc                <  2:1520631/142‑1 (MQ=255)
                                                                                                                                             ggCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGtat                                          <  2:87500/115‑1 (MQ=255)
                                                                                                                                             ggCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGa                        <  1:799855/133‑1 (MQ=255)
                                                                                                                                             ggCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGa                        >  2:799855/1‑133 (MQ=255)
                                                                                                                                              gCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCAc                                              >  4:1202655/1‑110 (MQ=255)
                                                                                                                                              gCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCAc                                              <  3:1202655/110‑1 (MQ=255)
                                                                                                                                                cGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaacatc     <  1:1397212/149‑1 (MQ=255)
                                                                                                                                                  aCAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCAt   <  1:921624/149‑1 (MQ=255)
                                                                                                                                                  aCAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCAt   >  1:983647/1‑149 (MQ=255)
                                                                                                                                                   cAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATg  <  2:353487/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATG  >  NC_000913/1850567‑1850862

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