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	97.7%	308.4 / 0.3	87	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (46/39);  total (47/40)
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

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281                                                                                                     |                                                                                                gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                  >  1:273278/1‑100 (MQ=255)  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                >  2:279890/1‑100 (MQ=255)    tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc                                                                                              >  2:445623/1‑100 (MQ=255)    tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc                                                                                              <  2:2717133/100‑1 (MQ=255)       cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑ccc                                                                                          <  1:2512408/100‑1 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccc                                                                                         <  2:978623/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccc                                                                                         <  1:1706802/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccc                                                                                         >  1:898241/1‑99 (MQ=255)               tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTg                                                                                   <  1:2746547/100‑1 (MQ=255)                cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTgg                                                                                  >  2:429557/1‑100 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaa                                                                                >  1:2219289/1‑100 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaa                                                                                >  2:520659/1‑100 (MQ=255)                  ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaa                                                                                <  2:2181260/100‑1 (MQ=255)                    gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaaaa                                                                              >  2:1885719/1‑100 (MQ=255)                     tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAAt                                                                             <  1:1203585/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATg                                                                            <  1:576767/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATg                                                                            <  2:58953/100‑1 (MQ=255)                       aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGa                                                                           <  2:2479741/100‑1 (MQ=255)                        cTTCACATAGCTTCAACTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGAc                                                                          <  1:1691312/100‑1 (MQ=255)                         tccgatctgCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTTTTCTGGAAAATGACt                                                                         <  1:2698890/92‑1 (MQ=255)                          tgagaTAGGTTCAACTTCTTGACACATAGTCCTCGTATCCTGCTGCCAATGCAAAGGAGAAGACTATGTTGAAAT‑CCCGCCCTGTTCTGGAAAATGACTa                                                                        <  1:2645450/96‑1 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAg                                                                       >  1:1769794/1‑100 (MQ=255)                            acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGc                                                                      <  2:1151629/100‑1 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGcc                                                                     >  1:55082/1‑100 (MQ=255)                                 gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTTTTCTGGAAAATGACTAgatcgga                                                                 >  2:2698914/1‑94 (MQ=255)                                 gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:2074650/1‑100 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                <  2:1769806/100‑1 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:969529/1‑100 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:737108/1‑100 (MQ=255)                                   ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:983320/1‑100 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  2:95935/100‑1 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:1011155/100‑1 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  2:5230/100‑1 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:455473/1‑100 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  2:2711539/1‑100 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  2:775674/1‑100 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:419217/1‑100 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:520657/100‑1 (MQ=255)                                          tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  2:1815903/100‑1 (MQ=255)                                           cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:2875199/1‑100 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:1502196/1‑100 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTCTTGc                                                   >  2:528320/1‑100 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:768337/100‑1 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  1:2468370/1‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:16294/1‑100 (MQ=255)                                                  caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  2:1043525/100‑1 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              >  2:797455/1‑100 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTCTGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCATGTTTGGTGCAGTATTGCTGttgttg                                           >  1:1266115/1‑100 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:2818710/100‑1 (MQ=255)                                                        cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTACTGTTGTTGa                                          >  2:1781753/1‑100 (MQ=255)                                                          tCGTATCCTGCTGCCAGTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  1:1256204/1‑100 (MQ=255)                                                           cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  2:983324/100‑1 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:1885705/100‑1 (MQ=255)                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  1:1612527/1‑100 (MQ=255)                                                              aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGCCTCGCCTGTCTGGTGCAGTATTGCTGTCGTTGAtctcga                                    >  1:2020355/1‑95 (MQ=255)                                                              aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  2:2293225/1‑100 (MQ=255)                                                              aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  2:2383942/100‑1 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:2396142/100‑1 (MQ=255)                                                                 ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:1483084/1‑100 (MQ=255)                                                                  tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  1:2236558/1‑100 (MQ=255)                                                                  tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                <  2:1275239/100‑1 (MQ=255)                                                                   gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  1:279889/100‑1 (MQ=255)                                                                    ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              >  2:1657230/1‑100 (MQ=255)                                                                    ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  1:606727/100‑1 (MQ=255)                                                                         aTTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        >  1:12256/1‑100 (MQ=255)                                                                          ttGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  1:509827/100‑1 (MQ=255)                                                                          ttGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        >  2:1108069/1‑100 (MQ=255)                                                                           tGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  1:413798/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2810165/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:2639813/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2571639/1‑100 (MQ=255)                                                                            gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2134205/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:2543444/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:939160/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:1502208/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTCGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:1840828/1‑100 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:2327131/100‑1 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:969533/100‑1 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCATTATTGCTGTTGTTGATTCCGATAACGCTGATTCGGCAg                   <  1:1781741/100‑1 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  2:898244/100‑1 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCCC‑CCCTGTTCTGGATAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  1:1907650/1‑100 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  2:2142524/1‑100 (MQ=255)                                                                                      agaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:968328/1‑100 (MQ=255)                                                                                      agaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:775671/100‑1 (MQ=255)                                                                                       gaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           <  2:1115223/100‑1 (MQ=255)                                                                                       gaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           <  2:2468392/100‑1 (MQ=255)                                                                                         cTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  2:2020571/1‑100 (MQ=255)                                                                                         cTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  2:2236576/100‑1 (MQ=255)                                                                                         cTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:2132040/1‑100 (MQ=255)                                                                                         cTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  2:852697/100‑1 (MQ=255)                                                                                          tATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        >  1:1703818/1‑100 (MQ=255)                                                                                           aTGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  1:364716/100‑1 (MQ=255)                                                                                            tGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  2:223655/100‑1 (MQ=255)                                                                                              ttGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:783438/1‑100 (MQ=255)                                                                                                gAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:1594470/1‑100 (MQ=255)                                                                                                     |                                                                                                GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  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: ‑