Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)6→7	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.7%	271.8 / ‑2.5	76	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (37/38);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281                                                                                                     |                                                                                                gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                  <  2:2066300/100‑1 (MQ=255) gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                  <  2:1157360/100‑1 (MQ=255)  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                  >  1:1086614/1‑99 (MQ=255)      gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑c                                                                                           >  1:552799/1‑100 (MQ=255)       cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑cc                                                                                          >  1:768123/1‑100 (MQ=255)       cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑cc                                                                                          <  2:644010/100‑1 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑ccc                                                                                         <  1:1386110/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑ccc                                                                                         <  2:1668346/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTACAAAGGAGAAGACTATGTTGAAAT‑cccc‑ccc                                                                                         <  2:1496052/100‑2 (MQ=255)          cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑ccctg                                                                                       >  1:936458/1‑97 (MQ=255)           gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑ccctgt                                                                                      >  2:1472714/1‑96 (MQ=255)           gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc‑ccctgt                                                                                      >  1:83999/1‑96 (MQ=255)             cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCCTTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTc                                                                                     <  1:1894162/100‑1 (MQ=255)               tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTg                                                                                   <  1:157287/100‑1 (MQ=255)                cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGCAGACTATGTTGAAATCCCCC‑CCTGTTCTgg                                                                                  >  1:1924447/1‑100 (MQ=255)                 aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGa                                                                                 <  1:373799/100‑1 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCAGATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGaa                                                                                <  1:1002032/100‑1 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGaa                                                                                >  2:885414/1‑100 (MQ=255)                     tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAAt                                                                             <  2:1815314/100‑1 (MQ=255)                     tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAAt                                                                             >  1:1366556/1‑100 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATg                                                                            <  2:2150328/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATg                                                                            <  1:1531187/100‑1 (MQ=255)                       aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGa                                                                           >  2:338876/1‑100 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGAc                                                                          <  1:1519661/100‑1 (MQ=255)                         ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACt                                                                         <  2:641964/100‑1 (MQ=255)                          tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACt                                                                         >  2:1453691/1‑99 (MQ=255)                          tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTa                                                                        >  1:281618/1‑100 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAg                                                                       <  1:1078375/100‑1 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGcc                                                                     <  2:369057/100‑1 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGcc                                                                     >  1:1626132/1‑100 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGcc                                                                     <  1:1755868/100‑1 (MQ=255)                               tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTg                                                                   >  1:1705372/1‑100 (MQ=255)                                aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGt                                                                  >  1:1934693/1‑100 (MQ=255)                                 gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:512861/1‑100 (MQ=255)                                   ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTa                                                                        >  1:372045/1‑91 (MQ=255)                                   ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:1979730/1‑100 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:384460/1‑100 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:347038/100‑1 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGAGTATGTTGAAAT‑CCCCGCCTGTTCTAGAAAATGACTAGCCTGTTTGGt                                                             >  1:1499747/1‑100 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:1204684/100‑1 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           >  1:2162911/1‑100 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:894021/100‑1 (MQ=255)                                          tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  1:862493/100‑1 (MQ=255)                                             tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                     <  2:1168794/100‑1 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:1982889/100‑1 (MQ=255)                                              cccACATAGTCATCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:1170109/1‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:1679582/1‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:944848/1‑100 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  1:1595010/100‑1 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:700857/100‑1 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:501109/100‑1 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  1:86457/1‑100 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:562067/100‑1 (MQ=255)                                                              aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  1:850919/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:935456/100‑1 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:790063/100‑1 (MQ=255)                                                                ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCCGTTGTTGATTCCGATa                                  >  1:1725812/1‑100 (MQ=255)                                                                 ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  1:1890650/1‑98 (MQ=255)                                                                  tgctgaCATTGCAAAGGAGAAGACTCTGTTGAACTCCCCC‑CCTGTTCTCGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                <  2:1215701/100‑1 (MQ=255)                                                                     tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:478133/100‑1 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:138743/1‑100 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:800546/1‑100 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:658812/1‑100 (MQ=255)                                                                         aTTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         <  2:119692/100‑1 (MQ=255)                                                                         aTTGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  2:1582840/1‑100 (MQ=255)                                                                           tGCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:21622/100‑1 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:1272983/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:1744682/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:944849/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:761331/1‑100 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:1230097/100‑1 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:915802/1‑100 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  2:1487802/1‑100 (MQ=255)                                                                                  gagaAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  1:1650940/100‑1 (MQ=255)                                                                                  gagaAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  2:552800/100‑1 (MQ=255)                                                                                      agaCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  2:1934703/100‑1 (MQ=255)                                                                                      agaCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:786309/1‑100 (MQ=255)                                                                                       gaCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  2:2203654/1‑100 (MQ=255)                                                                                        aCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:1170111/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:768124/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          >  1:1157150/1‑100 (MQ=255)                                                                                         cTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:1728577/1‑100 (MQ=255)                                                                                         cTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  2:2061463/100‑1 (MQ=255)                                                                                         cTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  2:1366558/100‑1 (MQ=255)                                                                                            tGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:2089005/1‑100 (MQ=255)                                                                                            tGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:1940220/1‑100 (MQ=255)                                                                                                gAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTTATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:1220296/100‑1 (MQ=255)                                                                                                gAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:695412/1‑100 (MQ=255)                                                                                                     |                                                                                                GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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