Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,193:1	(C)6→7	100%	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	94.4%	294.2 / 10.8	90	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (5/0);  new base C (45/40);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 6.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281                                                                                                     |                                                                                                 gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                   <  1:906728/100‑1 (MQ=255)     ctcTTCGACTTCACCGTCACTACACATAGCTTCAAATTCTTCCCGCATAGTGATCGTATCCTGCCGCTATTGCAAAGAAGAAGACTATGTTGAAAT‑‑cccc                                                                                             <  2:1590746/98‑1 (MQ=255)     cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cccc                                                                                             <  2:1392447/100‑1 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccc‑cc                                                                                         <  2:1360407/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccc‑cc                                                                                         <  1:1735795/100‑2 (MQ=255)          cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccc‑cctg                                                                                       >  1:294339/1‑97 (MQ=255)           gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccc‑cctgt                                                                                      <  1:2042552/100‑5 (MQ=255)            aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccc‑cctgtt                                                                                     >  2:1814021/1‑95 (MQ=255)            aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCC‑CTGTTc                                                                                     >  1:1229101/1‑100 (MQ=255)             cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTc                                                                                     >  2:2277438/1‑100 (MQ=255)             cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTc                                                                                     <  1:2054992/100‑1 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGaa                                                                                >  1:1596742/1‑100 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGaa                                                                                <  1:2356916/100‑1 (MQ=255)                       aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGa                                                                           <  1:1131916/100‑1 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGAc                                                                          >  2:885191/1‑100 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGAc                                                                          <  1:40346/100‑1 (MQ=255)                         ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACt                                                                         <  2:1039908/100‑1 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAGGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTa                                                                        >  2:1714230/1‑99 (MQ=255)                            acaTAGCTTCAAATTCTTCCCAGATAGTCTTTGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATg                                                                            <  2:2227634/94‑1 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                     >  2:654428/1‑100 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCC‑CTGTTCTGGAAAATGACTAGCCt                                                                    >  1:661350/1‑100 (MQ=255)                                 gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGtt                                                                 <  2:1352206/100‑1 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2181426/1‑100 (MQ=255)                                   ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:1688809/1‑100 (MQ=255)                                   ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:411769/100‑1 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:980972/1‑100 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:819574/1‑100 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:1958649/1‑100 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:500198/1‑100 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTACCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGCTCTGGCAAATGACTAACCTGTTTgat                                                             >  2:792801/1‑98 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:815092/1‑100 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  2:1698472/100‑1 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:1914937/100‑1 (MQ=255)                                         ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:2100033/100‑1 (MQ=255)                                           cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  1:801692/1‑100 (MQ=255)                                           cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:876600/1‑100 (MQ=255)                                            ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  1:2354876/1‑100 (MQ=255)                                            ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:2216904/1‑100 (MQ=255)                                            ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:1147482/1‑100 (MQ=255)                                             tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                     >  2:1829705/1‑100 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:26239/1‑100 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  2:550032/1‑100 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:696170/1‑100 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1762124/100‑1 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1759734/100‑1 (MQ=255)                                                   aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               <  1:202680/100‑1 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGTAGTATTGCtgttg                                              <  1:1787337/100‑1 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2122095/100‑1 (MQ=255)                                                     aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  1:1808561/100‑1 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:1201838/100‑1 (MQ=255)                                                        cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:1835432/1‑100 (MQ=255)                                                        cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:17318/1‑100 (MQ=255)                                                         ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:886227/1‑100 (MQ=255)                                                         ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:1802744/100‑1 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:2387553/100‑1 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  2:502474/1‑100 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:452058/1‑100 (MQ=255)                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  1:609166/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:263796/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:2359897/100‑1 (MQ=255)                                                                 ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 <  1:120758/100‑1 (MQ=255)                                                                 ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:194639/1‑100 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:427828/1‑100 (MQ=255)                                                                       ccATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCCGTACTGCTGTCGTTGATTCCGATAATACTGa                           >  2:216526/1‑100 (MQ=255)                                                                       ccATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:70042/1‑100 (MQ=255)                                                                       ccATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:2068968/1‑100 (MQ=255)                                                                          ttGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  1:1863711/100‑1 (MQ=255)                                                                          ttGCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  1:2029515/100‑1 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      <  1:1814036/100‑1 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:469941/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:1978691/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:1909893/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:794111/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:972358/1‑100 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:51812/100‑1 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:2186539/1‑100 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:70043/100‑1 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  1:125892/1‑100 (MQ=255)                                                                                   agaagaCTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  2:2057324/1‑99 (MQ=255)                                                                                       gaCTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  1:1201345/1‑100 (MQ=255)                                                                                       gaCTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  1:784668/1‑100 (MQ=255)                                                                                       gaCTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           <  2:2098334/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:1005178/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          >  1:736707/1‑100 (MQ=255)                                                                                         cTATGTTGAAATCCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:1888214/100‑1 (MQ=255)                                                                                         cTATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:1957476/100‑1 (MQ=255)                                                                                          tATGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        <  1:2277459/100‑1 (MQ=255)                                                                                            tGTTGAAATC‑CCCCC‑CTGTTCTGGAACATGACTAGACCGTTTGGTGCAGTGTTGCTGTAGTGGATTCCGATAATGCTGGTGCGGCAGGTGAGTGTGGAAc      >  2:2311899/1‑100 (MQ=255)                                                                                            tGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  1:2045350/100‑1 (MQ=255)                                                                                            tGTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  2:2356591/100‑1 (MQ=255)                                                                                             gTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:353200/1‑99 (MQ=255)                                                                                             gTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  1:1454013/100‑1 (MQ=255)                                                                                             gTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:689539/1‑100 (MQ=255)                                                                                             gTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:454876/1‑100 (MQ=255)                                                                                             gTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  2:45203/100‑1 (MQ=255)                                                                                             gTTGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  2:2354855/100‑1 (MQ=255)                                                                                               tGAAATC‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  2:2382581/1‑100 (MQ=255)                                                                                                 aaaT‑‑CCCCCTCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGCTTGTCGAACGTGc  >  2:2295181/1‑99 (MQ=255)                                                                                                     |                                                                                                 GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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