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	89.4%	254.2 / 27.0	85	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (5/4);  new base C (39/37);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCT  >  NC_000913/4638089‑4638282                                                                                                     |                                                                                                 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGTCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑c                                                                                                 <  2:2185260/100‑1 (MQ=255)   tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc                                                                                               >  2:897978/1‑100 (MQ=255)    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑cccc                                                                                              <  1:1660537/100‑1 (MQ=255)    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑cccc                                                                                              <  1:317035/100‑1 (MQ=255)     tcTTCTACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                             <  2:396169/99‑1 (MQ=255)     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCATAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                             >  2:1116635/1‑100 (MQ=255)     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGTGAAGACTATGTTGAAA‑T‑ccccc                                                                                             <  2:309999/100‑1 (MQ=255)     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                             >  2:981515/1‑100 (MQ=255)      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑c                                                                                           >  1:228659/1‑100 (MQ=255)       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑cc                                                                                          <  1:1536760/100‑2 (MQ=255)       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑cc                                                                                          <  2:1409668/100‑2 (MQ=255)         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑cctg                                                                                        <  2:705226/100‑4 (MQ=255)         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑cctg                                                                                        <  1:1832000/100‑4 (MQ=255)          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑cctgt                                                                                       >  1:1605082/1‑96 (MQ=255)            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCAATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc‑cctgtt                                                                                      >  2:212260/1‑94 (MQ=255)             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCC‑CTGTTCTg                                                                                    <  1:1822419/100‑1 (MQ=255)              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTg                                                                                    <  2:1880360/100‑1 (MQ=255)               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTgg                                                                                   >  2:321584/1‑100 (MQ=255)                aCCGTCACGTCACATAGCTTCAAATTCTTCCCACATAGTTTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAGC‑CCCCC‑CTGTTCTgt                                                                                   <  2:1568562/99‑2 (MQ=255)                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGaa                                                                                 <  2:2170490/100‑1 (MQ=255)                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGaaa                                                                                <  2:345153/100‑1 (MQ=255)                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGaaa                                                                                <  1:466747/100‑1 (MQ=255)                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGaaa                                                                                <  1:2146068/100‑1 (MQ=255)                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGaaa                                                                                <  2:842696/100‑1 (MQ=255)                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGaaat                                                                               >  2:1461439/1‑99 (MQ=255)                    tCACTTCACATAGCTTCAAATGCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAAt                                                                              <  1:742810/100‑1 (MQ=255)                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGAc                                                                           >  1:216811/1‑100 (MQ=255)                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTa                                                                         >  2:1055761/1‑100 (MQ=255)                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                      >  2:625153/1‑100 (MQ=255)                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                      <  2:251637/100‑1 (MQ=255)                             aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                      <  2:1218069/99‑1 (MQ=255)                              tAGCTTCAAATTCTTCCCACATAGACTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTg                                                                    >  2:314328/1‑100 (MQ=255)                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGtt                                                                  >  2:784493/1‑100 (MQ=255)                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCATGtt                                                                  >  2:1744862/1‑100 (MQ=255)                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                              <  2:1144484/100‑1 (MQ=255)                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGAATATGTTGAAA‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCATGTTTGGt                                                              >  1:1332233/1‑100 (MQ=255)                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                            >  2:510048/1‑100 (MQ=255)                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                           <  2:1891511/100‑1 (MQ=255)                                            tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGG‑AAATGACTAGCCTGTTTGGTGCAGTATTg                                                     <  1:693969/100‑1 (MQ=255)                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                     >  1:1263293/1‑100 (MQ=255)                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGCGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                    >  1:331816/1‑100 (MQ=255)                                              ccAAATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTTGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                      <  2:400856/98‑1 (MQ=255)                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                   <  2:1396506/100‑1 (MQ=255)                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                   >  1:218871/1‑100 (MQ=255)                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                  <  1:591100/99‑1 (MQ=255)                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                                >  1:1059013/1‑100 (MQ=255)                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                                <  1:759892/100‑1 (MQ=255)                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                              >  2:1851730/1‑100 (MQ=255)                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                              <  1:1630929/100‑1 (MQ=255)                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                              <  2:935396/100‑1 (MQ=255)                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgtn                                             <  2:477786/100‑2 (MQ=255)                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                             >  1:823841/1‑100 (MQ=255)                                                     gTCTTCGTATCCTGCTGCCATTACAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                             >  1:1222404/1‑100 (MQ=255)                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                            >  2:1581369/1‑100 (MQ=255)                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                            <  2:1635054/100‑1 (MQ=255)                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                           <  1:1196200/100‑1 (MQ=255)                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                          >  1:1077559/1‑100 (MQ=255)                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                      <  2:1483684/100‑1 (MQ=255)                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                  >  2:112488/1‑100 (MQ=255)                                                                 tgctgcTATTGCAAAGGAGAAGACTGTGTTGAAA‑TCCCCCC‑CTGTTCTTGAAAATGACTAGCCTGTTTGGTTCAGTATTGCTGTTGTTGATTCCGATAAt                                 >  2:1019046/1‑100 (MQ=255)                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                   <  2:1041624/98‑1 (MQ=255)                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAAGAACTAGCCTGTTTGGTGCAGTATTGCTTTTGTTGATTCCGATAATg                                >  2:361622/1‑100 (MQ=255)                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATgtt                              <  1:1520065/100‑3 (MQ=255)                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCA‑CTGTTCTGGAAAATCACTAGCCTGTTTGGTGCAGAATTGCTGTTGTTGATTCCGATAATGCTGa                            >  2:1914161/1‑100 (MQ=255)                                                                     gcCATTGCAAAGGAGAAGACTATATTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                               <  2:2225039/98‑1 (MQ=255)                                                                      ccATTTCAAAGGAGAAGACTATGTTTATA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                            <  2:468504/100‑1 (MQ=255)                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                            >  2:2090125/1‑100 (MQ=255)                                                                       cATTGCAAAGGATAAGTCTATGTTTAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                            <  2:1272965/99‑1 (MQ=255)                                                                       cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                            >  1:583864/1‑99 (MQ=255)                                                                       cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                           >  1:83310/1‑100 (MQ=255)                                                                          tGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                        >  2:696640/1‑100 (MQ=255)                                                                          tGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                        <  2:840464/100‑1 (MQ=255)                                                                          tGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                        >  2:1670620/1‑100 (MQ=255)                                                                           ncAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                       >  2:390301/2‑100 (MQ=255)                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                       >  2:1605225/1‑100 (MQ=255)                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                        >  1:28703/1‑99 (MQ=255)                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                      >  1:80506/1‑100 (MQ=255)                                                                             aaaGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                     <  1:390300/100‑1 (MQ=255)                                                                              aaGGAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                    <  1:891359/100‑1 (MQ=255)                                                                                ggAGAAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCTATAATGCTGATTCGGCAGGt                  <  2:610666/100‑1 (MQ=255)                                                                                 gagaAGACTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                 >  2:1114087/1‑100 (MQ=255)                                                                                   gaagaCTATGTTGAAA‑TCCCCCC‑CTTTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt               <  2:1806538/100‑1 (MQ=255)                                                                                   gaagaCTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt               <  1:225554/100‑1 (MQ=255)                                                                                   gaagaCTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt               <  2:1477356/100‑1 (MQ=255)                                                                                   gaagaCTATGTTGAAA‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt              >  2:1050215/1‑100 (MQ=255)                                                                                     agaCTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg             >  2:205022/1‑100 (MQ=255)                                                                                      gaCTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTTTTGATTCCGATAATGCTGATTCGGCAGGTGTTTGt            >  2:38934/1‑100 (MQ=255)                                                                                       aCTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc           <  1:481873/100‑1 (MQ=255)                                                                                        cTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc           >  1:1481571/1‑99 (MQ=255)                                                                                        cTATGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg          >  2:145892/1‑100 (MQ=255)                                                                                        cTATGTTGAAA‑T‑CCCCCACTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGATGTTGTTGATTCCGATAATGTTGATTCGGAAGGTGATTGt            >  2:294625/1‑98 (MQ=255)                                                                                           tGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa         >  1:2138452/1‑98 (MQ=255)                                                                                           tGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc       >  2:2107930/1‑100 (MQ=255)                                                                                           tGTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc       <  1:2211805/100‑1 (MQ=255)                                                                                            gTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc       >  2:2224986/1‑99 (MQ=255)                                                                                            gTTGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg      <  2:923594/100‑1 (MQ=255)                                                                                             ttGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATACTGCTGATTCGGCAGGTGATTGTCGAACg      >  2:758395/1‑99 (MQ=255)                                                                                             ttGAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt     >  2:1665797/1‑100 (MQ=255)                                                                                               gAAA‑TCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc   <  2:1059011/100‑1 (MQ=255)                                                                                                aaa‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCTATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCt  <  2:1332227/99‑1 (MQ=255)                                                                                                     |                                                                                                 CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCT  >  NC_000913/4638089‑4638282

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