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	95.5%	295.6 / 7.9	88	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (45/39);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638088‑4638280                                                                                                     |                                                                                              gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                >  2:401609/1‑100 (MQ=255) gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                >  1:2823085/1‑100 (MQ=255) gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                >  1:2206742/1‑100 (MQ=255) gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                >  1:1996204/1‑100 (MQ=255)    tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc                                                                                            <  2:1629638/100‑1 (MQ=255)     cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc                                                                                           <  2:720925/100‑1 (MQ=255)     cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc                                                                                           <  2:750469/100‑1 (MQ=255)      gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc                                                                                          >  1:1059341/1‑100 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccc                                                                                        <  2:1359460/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccc                                                                                        <  2:1522718/100‑2 (MQ=255)          cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAACGGAGAAGACTATGTTGAAAT‑CCCCCCTGt                                                                                      >  1:2301676/1‑100 (MQ=255)            aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccctgtt                                                                                    >  2:1447330/1‑95 (MQ=255)            aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccctgtt                                                                                    >  2:2665734/1‑95 (MQ=255)              ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCt                                                                                   <  2:2717171/100‑1 (MQ=255)               tcgccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTg                                                                                  <  2:124915/97‑1 (MQ=255)               tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTTCAAAGTAGAAGACTATGTGGAAATCCCCCCCTGTTCTg                                                                                  <  1:1666321/100‑1 (MQ=255)                cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTgg                                                                                 >  2:1601183/1‑100 (MQ=255)                    gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGaaaa                                                                             <  2:2617980/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATg                                                                           <  1:2709001/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATg                                                                           <  1:737994/100‑1 (MQ=255)                      cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATg                                                                           <  2:2250898/100‑1 (MQ=255)                      cACTTCACATACATTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATg                                                                           <  1:1761807/100‑1 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGAc                                                                         >  2:1027714/1‑100 (MQ=255)                         ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACt                                                                        <  1:2843358/100‑1 (MQ=255)                         ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACt                                                                        >  1:829798/1‑100 (MQ=255)                          tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAg                                                                      >  1:2608319/1‑100 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAg                                                                      >  2:2155916/1‑100 (MQ=255)                            acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGc                                                                     >  2:1340089/1‑100 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    >  2:2890465/1‑100 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  2:1482813/100‑1 (MQ=255)                               tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTg                                                                  <  1:2074125/100‑1 (MQ=255)                               tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTg                                                                  <  2:482594/100‑1 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                               >  1:1116432/1‑100 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                               <  1:364828/100‑1 (MQ=255)                                  cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                               >  1:402094/1‑100 (MQ=255)                                    tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             >  1:510020/1‑100 (MQ=255)                                       aaTTCTTCCCACATAGTCTTGGCATCCTGCTGCCATTTCAAAGGAGAGGACGATGTGGAAATCCCCCCCTGTTCTGTAAAATGACTAGCCTGTGTGGTGc                                                          >  1:1693051/1‑100 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                          <  1:1867472/100‑1 (MQ=255)                                       aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                          <  1:675746/100‑1 (MQ=255)                                        aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                          <  1:1719201/99‑1 (MQ=255)                                           cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                      >  2:824269/1‑100 (MQ=255)                                             tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                    >  2:693792/1‑100 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                   >  2:2174911/1‑100 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                  >  1:398476/1‑100 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                  >  1:664854/1‑100 (MQ=255)                                                nacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  1:158201/2‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  2:71124/1‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  1:421156/1‑100 (MQ=255)                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGACTGTTTGGTGCAGGATTg                                                   >  1:358122/1‑98 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                >  2:1397594/1‑100 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                <  1:666394/100‑1 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                >  2:2507603/1‑100 (MQ=255)                                                   aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                             >  1:1874228/1‑100 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           >  2:1383626/1‑100 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           >  1:2610380/1‑100 (MQ=255)                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  1:978062/100‑1 (MQ=255)                                                       tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                          <  1:1601158/100‑1 (MQ=255)                                                        cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                         >  2:2153689/1‑100 (MQ=255)                                                         ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATATTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        <  2:2042664/100‑1 (MQ=255)                                                          tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                       >  1:1704269/1‑100 (MQ=255)                                                           cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGACTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                      >  1:2169394/1‑100 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     >  1:1477023/1‑100 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     >  2:2703596/1‑100 (MQ=255)                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                    <  1:2526772/100‑1 (MQ=255)                                                                     tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  1:2201457/100‑1 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                           <  1:825476/100‑1 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTACTGCTGTTGTTGATTCCGATAATGCTg                           >  2:1477622/1‑100 (MQ=255)                                                                        cATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                         <  2:1683138/100‑1 (MQ=255)                                                                           tGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                      >  2:2604626/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    >  1:1570330/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    <  2:2796085/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    <  2:2792939/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    >  1:15558/1‑100 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    <  1:2884433/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    >  1:2882693/1‑100 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                  >  1:1404629/1‑100 (MQ=255)                                                                               aaGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                  >  1:2802818/1‑100 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                <  2:2084765/100‑1 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                >  2:1274657/1‑100 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                >  2:1148172/1‑100 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                >  1:2091800/1‑100 (MQ=255)                                                                                  gagaAGACTATGTTGAGATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg               <  1:208122/100‑1 (MQ=255)                                                                                  gagaAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg               <  1:1993069/100‑1 (MQ=255)                                                                                  gagaAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg               >  1:2192433/1‑100 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt             >  2:800223/1‑100 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt             >  1:1500979/1‑100 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt             >  1:1422701/1‑100 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt             <  1:1340071/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         <  2:2806516/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         <  2:509619/100‑1 (MQ=255)                                                                                        aCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         >  2:1446982/1‑100 (MQ=255)                                                                                        aCTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg        <  2:2608352/100‑1 (MQ=255)                                                                                         cTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg        <  2:2706264/100‑1 (MQ=255)                                                                                          tATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa       <  1:401604/100‑1 (MQ=255)                                                                                            tGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc     <  2:664860/100‑1 (MQ=255)                                                                                             gTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    >  1:1060359/1‑100 (MQ=255)                                                                                             gTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    <  2:2206769/100‑1 (MQ=255)                                                                                             gTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    <  1:292185/100‑1 (MQ=255)                                                                                               tGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg  >  2:1859558/1‑100 (MQ=255)                                                                                               tGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg  <  2:829804/100‑1 (MQ=255)                                                                                                     |                                                                                              GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638088‑4638280

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: ‑