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	98.6%	297.4 / ‑4.6	87	S4P (TCC→CCC)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (49/34);  minor base A (1/1);  total (50/37)
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

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281                                                                                                      |                                                                                                gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑t                                                                                                  <  2:1616515/100‑1 (MQ=255)  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c                                                                                               >  2:545867/1‑100 (MQ=255)  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACCATGTTGAAA‑T‑‑c                                                                                               <  1:1524342/100‑1 (MQ=255)     cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑cccn                                                                                            <  2:1560971/100‑2 (MQ=255)      ncTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑ccccc                                                                                           >  1:1318795/2‑100 (MQ=255)       cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑cccccc                                                                                          >  2:2204245/1‑100 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑ccccccc                                                                                         <  2:198818/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑ccccccc                                                                                         <  2:870865/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑ccccccc                                                                                         >  2:1142657/1‑99 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑ccccccc                                                                                         <  2:395416/100‑2 (MQ=255)         tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑ccccccct                                                                                        >  1:506118/1‑98 (MQ=255)                cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTGCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTc                                                                                     <  1:2327692/97‑1 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGCGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGa                                                                                 >  1:1385293/1‑99 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGaa                                                                                >  1:1019453/1‑100 (MQ=255)                   cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑CCCCCCTGTTCTGGaaaa                                                                              <  1:2374299/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATg                                                                            <  2:40605/100‑1 (MQ=255)                       aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGa                                                                           >  1:90999/1‑100 (MQ=255)                       aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGa                                                                           >  2:886980/1‑100 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCCTGTTCTGGAAAATGa                                                                           >  2:1051861/1‑100 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGa                                                                           >  1:2648026/1‑99 (MQ=255)                          tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTa                                                                        >  1:768508/1‑100 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGcc                                                                     <  1:1641051/100‑1 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGcc                                                                     >  2:2384583/1‑100 (MQ=255)                                 gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 <  1:1787915/100‑1 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2299871/1‑100 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:1600610/1‑100 (MQ=255)                                   ttCAAATTCTTCCCACATAGTCCTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGCCTCGCCTGTTTg                                                               >  2:2087308/1‑100 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTGTCATGCTGCCATTGTAAAAGAGAAGACTGTGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAAAGAATAg                                                                       >  1:2556082/1‑91 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGATAAGACTATGTTGAAA‑TA‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  1:2863638/99‑1 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:1410286/100‑1 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             <  2:772668/100‑1 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             <  2:1877699/100‑1 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:2483597/1‑100 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  2:354260/1‑100 (MQ=255)                                        aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          <  1:190906/100‑1 (MQ=255)                                         ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:232253/100‑1 (MQ=255)                                         ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:2568356/100‑1 (MQ=255)                                         ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:2244202/100‑1 (MQ=255)                                           cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:167044/1‑100 (MQ=255)                                            atctcacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      <  1:2312694/96‑1 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:1061759/100‑1 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TA‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  2:2605611/1‑100 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:1576801/1‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtagat                                                  >  2:2312709/1‑96 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:2106712/100‑1 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  2:2723514/1‑100 (MQ=255)                                                  caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGAGTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  2:1772568/100‑1 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:1849841/100‑1 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2207535/100‑1 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:127825/1‑100 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:1788482/1‑100 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:77482/100‑1 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:2582257/100‑1 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:962837/100‑1 (MQ=255)                                                         ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  2:269575/1‑100 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAGC‑‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCGGTATTGCTGTTGTTGATTcc                                      <  1:1764202/100‑1 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:171959/1‑100 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:2082164/1‑100 (MQ=255)                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTTAAC‑TC‑CCCCCCTTTTCTGTAAAATGACTAGCCTGTTTTGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:876509/99‑1 (MQ=255)                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:1814486/1‑100 (MQ=255)                                                              aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  1:1302107/100‑1 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:1762591/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:2861791/100‑1 (MQ=255)                                                                 ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  2:1320937/1‑100 (MQ=255)                                                                 ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:375682/1‑100 (MQ=255)                                                                   gctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  1:1399339/1‑99 (MQ=255)                                                                   gctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  2:1108433/1‑100 (MQ=255)                                                                    ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:1477991/100‑1 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1411865/1‑100 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:469232/1‑100 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  1:2790159/100‑1 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:1316620/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:2367001/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATCGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:1020104/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGGGA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGTTGCAGTATTGCTGTTGTTGATTCCGATCATGCTGAtt                         <  1:547940/96‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTCTTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:53097/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:2082177/100‑1 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:911258/100‑1 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:403702/100‑1 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  1:1409770/1‑100 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  1:2321377/100‑1 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCCGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  >  1:556880/1‑99 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:2144107/100‑1 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGGTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:2041803/1‑100 (MQ=255)                                                                                    gaagaCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:2299886/100‑1 (MQ=255)                                                                                      agaCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:2401019/1‑100 (MQ=255)                                                                                      agaCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:1849462/1‑100 (MQ=255)                                                                                      agaCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:440297/1‑100 (MQ=255)                                                                                      agaCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:1253617/1‑100 (MQ=255)                                                                                       gaCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  1:1394540/1‑100 (MQ=255)                                                                                        aCTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:2199807/100‑1 (MQ=255)                                                                                         cTATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  2:710636/1‑100 (MQ=255)                                                                                          tATGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        >  1:142542/1‑100 (MQ=255)                                                                                           aTGTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  1:2201612/100‑1 (MQ=255)                                                                                             gTTGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:67812/1‑100 (MQ=255)                                                                                              ttGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2720642/1‑100 (MQ=255)                                                                                               tGAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  2:1327736/1‑100 (MQ=255)                                                                                                gAAA‑TC‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGCTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  1:1196779/1‑100 (MQ=255)                                                                                                      |                                                                                                GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  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 < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑