Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,313:1	(C)7→8	100%	pseudogene (675/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,306	1	.	C	92.9%	269.6 / 14.8	84	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (44/33);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597208‑1597400                                                                                                      |                                                                                                   gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTaacca‑a‑‑c                                                                                                  <  1:144374/100‑1 (MQ=255)  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTaacca‑a‑‑c‑c                                                                                                >  2:2971137/1‑100 (MQ=255)   gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c                                                                                              <  2:2064645/100‑1 (MQ=255)   gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c                                                                                              >  2:337541/1‑100 (MQ=255)      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGACACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c‑cc                                                                                           >  2:3086813/1‑99 (MQ=255)       gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c‑cccc                                                                                         >  2:3132467/1‑100 (MQ=255)          tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c‑ccccctc                                                                                      >  1:1914963/1‑97 (MQ=255)           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c‑ccccctcg                                                                                     <  1:2663927/100‑5 (MQ=255)            ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c‑ccccct                                                                                       >  2:869880/1‑95 (MQ=255)             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATCCCGGTGGAAGTGCCACCGATAAT‑CGTATTTTCAGCCGTAACTA‑AC‑C‑C‑C‑CCCCTcgcg                                                                                    <  2:125830/100‑1 (MQ=255)              gtcatgaCAGAGTCTTGATTCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTACCCA‑A‑‑c‑c‑c‑ccccctcccca                                                                                  >  2:3217018/4‑93 (MQ=255)              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGa                                                                                   >  1:1207104/1‑100 (MQ=255)              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGa                                                                                   >  2:909033/1‑100 (MQ=255)              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑c‑c‑c‑cccccctcgcg                                                                                  <  1:2965072/100‑8 (MQ=255)               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGaa                                                                                  <  1:147035/100‑1 (MQ=255)                gatgaCAGAGTCTTGCTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAAt                                                                                 <  2:1718370/100‑1 (MQ=255)                gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAAt                                                                                 <  1:3784379/100‑1 (MQ=255)                 atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGCATa                                                                                >  2:1729616/1‑100 (MQ=255)                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAg                                                                               <  2:1289956/100‑1 (MQ=255)                    aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTc                                                                             >  1:3110954/1‑100 (MQ=255)                     gagggTCTTGACCGTGCTCATTGAGGGTTCCGATACCGGTTGCAGTGCCACCGATAAT‑CGTATTTTCAGCCGGAACCAGC‑‑C‑C‑C‑CCCCTcgcg                                                                                    <  2:3270162/88‑1 (MQ=255)                     cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCg                                                                            <  2:3308453/100‑1 (MQ=255)                       gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGcc                                                                          >  2:675593/1‑100 (MQ=255)                          tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTcc                                                                       >  1:3190196/1‑100 (MQ=255)                           cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCCGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTccc                                                                      >  2:2965832/1‑100 (MQ=255)                           cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTccc                                                                      >  1:2432894/1‑100 (MQ=255)                            ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTcccc                                                                     <  1:610384/100‑1 (MQ=255)                            ttGATCCGGCACATTCAGGGTTCCGGTACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑C‑CAC‑CCCCTCGCGAATAgctgcccgccc                                                                     >  2:1589393/1‑90 (MQ=255)                             tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCCATAAT‑CGTATTTTCAGCCGTAACCA‑ACCC‑C‑C‑CCCCTCGCGAATAGTCGCCTcccc                                                                     <  2:3601499/100‑1 (MQ=255)                              gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAg                                                                   <  1:2197643/100‑1 (MQ=255)                               aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGt                                                                  <  1:2262847/100‑1 (MQ=255)                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                >  2:1024117/1‑100 (MQ=255)                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                >  1:3410618/1‑100 (MQ=255)                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                              >  2:191880/1‑100 (MQ=255)                                      acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  2:2642472/99‑1 (MQ=255)                                       ccttCAGGGTTCCCATTCCGGTGGCCATGCCCCCCCTTATCCGTATTTCCGCCCTTACCCA‑CC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  2:2416476/97‑1 (MQ=255)                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  1:1993985/100‑1 (MQ=255)                                         ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                          >  1:1544518/1‑98 (MQ=255)                                         ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                        >  2:466197/1‑100 (MQ=255)                                         ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                        >  2:3218088/1‑100 (MQ=255)                                           cAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                      <  1:92946/100‑1 (MQ=255)                                               gTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                  <  2:1623851/100‑1 (MQ=255)                                                 tCCGATACCGGTGGCATCGCCACCGATAAT‑CGTATTTTCAGCCGTAACCACC‑‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                <  1:282294/100‑1 (MQ=255)                                                 acgattatcGGTGGCAATGCCACCGGTAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑CAC‑C‑CCCCTCGCGAATAGTCGGCTCCCCAGTTCTTTGATTATAAATTTg                                                >  1:3417477/9‑100 (MQ=25)                                                   gattatcGGTGGCAATGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑C‑C‑CACCCCTCGCGAATAATCGCCTCCACAGTTCCTTGATTACCCATTTGaa                                              >  1:2222485/7‑100 (MQ=255)                                                    attatcGGTGGCAATGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGATCCTTGCTTACCCATTTGaaa                                             >  2:3645500/6‑100 (MQ=255)                                                     aTACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                            <  1:2423035/100‑1 (MQ=255)                                                     aTACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                            <  1:1657350/100‑1 (MQ=255)                                                      tatcGGTGGCAATGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                           >  1:2531812/4‑100 (MQ=255)                                                       atcGGTGGCAATGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          >  2:2566663/3‑100 (MQ=255)                                                        ccGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  2:1631801/1‑100 (MQ=255)                                                         cGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCCATTTGAAATTCaa                                        >  1:1120439/1‑100 (MQ=255)                                                           gTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                      <  1:1616162/100‑1 (MQ=255)                                                           gTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                      >  1:1767955/1‑100 (MQ=255)                                                            tGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACTA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg                                     <  2:3346595/100‑1 (MQ=255)                                                             ggCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                       <  2:2465697/97‑1 (MQ=255)                                                             ggCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                    >  1:2793571/1‑100 (MQ=255)                                                              gCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGan                                   <  1:1441483/100‑2 (MQ=255)                                                               cATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:1706592/100‑1 (MQ=255)                                                               cATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:1134159/100‑1 (MQ=255)                                                                   gCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AGCC‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                                >  2:2485437/1‑99 (MQ=255)                                                                   gCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGTCTCCCCAATTCCTTGATTACCAATTTCAAATTCAATTGAGGa                                 >  2:3368337/1‑97 (MQ=255)                                                                   gCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                              >  1:896704/1‑100 (MQ=255)                                                                   gCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                              >  2:74655/1‑100 (MQ=255)                                                                    acACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTACTTGATTACCCATTTGAACTTCACTTGAggcctca                             >  1:779750/2‑95 (MQ=255)                                                                     cACCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCCATTTGAAATTCAATTGAGGAATCat                            >  2:718216/1‑100 (MQ=255)                                                                      aCCGATAAT‑CGTATTTTCAGCCGTAACCA‑ACCC‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                            <  2:972339/100‑1 (MQ=255)                                                                      aCCGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCTCCTCCCCAGTTCCTTGATTCCCAATTTGAAATTCAATTGAGGAATCatt                           >  2:1037808/1‑100 (MQ=255)                                                                        cGATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATGattat                         >  1:940737/1‑100 (MQ=255)                                                                         gATAAT‑CGTATTTTCAGCCGTAACCA‑ACCC‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                         <  1:1254488/100‑1 (MQ=255)                                                                         gATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattatt                        >  1:22790/1‑100 (MQ=255)                                                                         gATAAT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattatt                        >  1:1926536/1‑100 (MQ=255)                                                                            aaT‑CGTTTTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     <  1:973978/100‑1 (MQ=255)                                                                            aaT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     <  1:986574/100‑1 (MQ=255)                                                                            aaT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATCCAATTGAGGAATCATTAttttt                     >  2:331073/1‑100 (MQ=255)                                                                             aT‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  1:2025028/1‑100 (MQ=255)                                                                              t‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   >  1:465048/1‑100 (MQ=255)                                                                              t‑CGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  1:1788984/100‑1 (MQ=255)                                                                                cGTATTTTCAGCCGTAACCA‑ACCC‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  1:762230/1‑99 (MQ=255)                                                                                cGTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  >  1:1663748/1‑100 (MQ=255)                                                                                 gTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  1:1615375/1‑100 (MQ=255)                                                                                 gTATTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  1:1053774/1‑100 (MQ=255)                                                                                   aTTTTCAGCCGTAACCATC‑‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               <  1:1235732/100‑1 (MQ=255)                                                                                   aTTTTCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGACGAATCATTATTTTTTATTAg               <  2:1439744/100‑1 (MQ=255)                                                                                     tttCAGCCGTAACCA‑ACCC‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:2109442/100‑1 (MQ=255)                                                                                     tttCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             <  2:777163/100‑1 (MQ=255)                                                                                     tttCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             >  1:3595250/1‑100 (MQ=255)                                                                                     tttCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             >  2:2490425/1‑100 (MQ=255)                                                                                      ttCAGCCGTAACCA‑ACCC‑C‑C‑CCCCTCGCGAATAGTCGCCCCCCCCGTTCCTTGATTACCAATTTGGAATTCAATTTAGGAATCATTATTTTTTTTTAGcc             >  1:3714320/1‑100 (MQ=255)                                                                                      ttCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCa            >  1:135788/1‑100 (MQ=255)                                                                                      ttCAGCCGTAACCA‑AC‑C‑C‑C‑CCCCCCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGAAAATCAATTTAGGAAACATTATTTTTTATTTGcc             >  1:191728/1‑99 (MQ=255)                                                                                          gCCGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  2:1177078/1‑99 (MQ=255)                                                                                           ccGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc       <  1:3490784/100‑1 (MQ=255)                                                                                            cGTAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc      >  2:128790/1‑100 (MQ=255)                                                                                              tAACCA‑AC‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    <  2:2417643/100‑1 (MQ=255)                                                                                               aacca‑ac‑c‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCGAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    <  2:2566692/99‑1 (MQ=255)                                                                                               aacca‑ac‑c‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    >  2:35583/1‑99 (MQ=255)                                                                                                acca‑ac‑c‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGCATCATTATTTTTTATTAGCCACTCGCCAccc  >  1:647466/1‑100 (MQ=255)                                                                                                      |                                                                                                   GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAT‑CGTATTTTCAGCCGTAACCA‑A‑‑C‑C‑C‑CCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597208‑1597400

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑