Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,777,363	T→C	100%	F83F (TTT→TTC)	yfjZ →	CP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,777,363	0	T	C	98.0%	315.4 / ‑4.0	96	F83F (TTT→TTC)	yfjZ	CP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (54/39);  total (56/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 = 9.93e-01

TATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGA  >  NC_000913/2777266‑2777456                                                                                                  |                                                                                              tATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCAc                                                                                             <  2:175381/100‑1 (MQ=255)  aTTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGACCAATGCGTCACCCTGTACCACAACGGTTTCAcc                                                                                            >  1:1681678/1‑100 (MQ=255)    ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTg                                                                                          <  1:1896746/100‑1 (MQ=255)      ccACATTTTATCAGCCAGATAGAGTCGATGCTGAGCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                        <  2:1352294/100‑1 (MQ=255)      ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                        >  1:2229405/1‑100 (MQ=255)       cacaTTTTATCAGCCAGATAGAGTCGATGCTTACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGa                                                                                       <  2:2548230/100‑1 (MQ=255)        acaTTTTGTCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGaa                                                                                      >  1:1348306/1‑100 (MQ=255)         caTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACGCTGTACCACAACGGTGTCACCTGCGa                                                                                       >  2:1100010/1‑98 (MQ=255)           ccgatctaGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGcc                                                                                   <  1:167773/93‑1 (MQ=255)            tttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGATCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCg                                                                                  <  2:1807732/100‑1 (MQ=255)            tttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCg                                                                                  >  1:916323/1‑100 (MQ=255)               aTCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATa                                                                               >  2:329431/1‑100 (MQ=255)                 cAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATAct                                                                             >  1:1374369/1‑100 (MQ=255)                  aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGccagctcgg                                                                            >  2:167771/1‑93 (MQ=255)                  aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActc                                                                            >  1:365466/1‑100 (MQ=255)                  aGCCAGATAGAGTAGATGCTGACCACTGGTGAACTGAACCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATAct                                                                             >  2:1581781/1‑99 (MQ=255)                   gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActct                                                                           >  2:1833766/1‑100 (MQ=255)                    ccAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCtt                                                                          <  1:1066911/100‑1 (MQ=255)                      aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTgg                                                                        >  1:2440529/1‑100 (MQ=255)                      aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTgg                                                                        <  2:2515489/100‑1 (MQ=255)                       gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                       <  2:1315262/100‑1 (MQ=255)                       gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                       <  1:356024/100‑1 (MQ=255)                       gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                       <  2:911547/100‑1 (MQ=255)                        aTAGAGTCGATGCTGACCACTTGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTa                                                                      <  2:1305864/100‑1 (MQ=255)                            agTGGATGCTGACGACTGGTGTACTGAATCTGCGCCATGCCCAATGCGTCATTTTGTTGCATAACGGTTTGACCTGCGACGCGGGTACTTTTGGTAGtt                                                                   >  2:1325377/1‑99 (MQ=255)                              tCGATGGTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                                <  1:1026726/100‑1 (MQ=255)                               cGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCgg                                                               >  1:2093404/1‑100 (MQ=255)                                gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                              <  2:976574/100‑1 (MQ=255)                                 aTGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGTAGCCGATACTCTTGGTAGTTGCGGCt                                                             <  1:1558057/100‑1 (MQ=255)                                   gCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTAc                                                           >  2:133234/1‑100 (MQ=255)                                    cTGACCACTGGTGCACTGAATCTACGACATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCAGCGAAGCCGATACTCCTGGTAGTCGCGGCTACg                                                          >  2:395654/1‑100 (MQ=21)                                    cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                          <  1:2215622/100‑1 (MQ=255)                                    cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                          >  2:1992822/1‑100 (MQ=255)                                    cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTCCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                          >  2:857508/1‑100 (MQ=255)                                     tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                         >  2:195535/1‑100 (MQ=255)                                     tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                         >  2:418173/1‑100 (MQ=255)                                        ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                      <  1:1957047/100‑1 (MQ=255)                                           cTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACAt                                                   >  1:1275840/1‑100 (MQ=255)                                           cTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACAt                                                   >  1:914014/1‑100 (MQ=255)                                            tGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATc                                                  >  1:2485754/1‑100 (MQ=255)                                             ggTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCg                                                 >  2:196252/1‑100 (MQ=255)                                              gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGc                                                <  1:1695148/100‑1 (MQ=255)                                               ttaaCTGAATCCCCGCCATGCCCACTGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCg                                                 <  2:2480386/96‑1 (MQ=255)                                               tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCt                                               >  1:500075/1‑100 (MQ=255)                                               tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTATCACCTGCGAAGCCGATCCTCTTGGTAGTTGCGTCTACGTATACATCGc                                                >  2:2201121/1‑99 (MQ=21)                                               tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTATACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCt                                               >  1:1801379/1‑100 (MQ=255)                                                gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                              >  1:2512007/1‑100 (MQ=255)                                                 aaCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                              <  1:624038/99‑1 (MQ=255)                                                 aaCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                             >  2:1814143/1‑100 (MQ=255)                                                  aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATCCTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                            >  1:51045/1‑100 (MQ=255)                                                  aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                            >  1:922968/1‑100 (MQ=255)                                                  aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                            >  1:1256005/1‑100 (MQ=255)                                                  aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                            >  2:14442/1‑100 (MQ=255)                                                   ccgAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                           <  1:2494794/98‑1 (MQ=255)                                                    tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTa                                          >  2:795177/1‑100 (MQ=255)                                                     gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                         >  2:914751/1‑100 (MQ=255)                                                       aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                       >  1:1224777/1‑100 (MQ=255)                                                        tCCCCGCCATGCCCAATGCGTCACCCTGTACCACCACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                      <  2:886987/100‑1 (MQ=255)                                                        tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                      >  1:1151933/1‑100 (MQ=255)                                                         ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                     >  2:872053/1‑100 (MQ=255)                                                         ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTGGGTAGATGCGGCTATGTATACATCGCTGTTTACCCCa                                     >  2:1729072/1‑100 (MQ=21)                                                          cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccccc                                    >  2:52561/1‑98 (MQ=255)                                                          cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCAc                                    >  1:1298733/1‑100 (MQ=255)                                                           ccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACt                                   <  1:1020781/100‑1 (MQ=255)                                                           ccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACt                                   <  1:1301140/100‑1 (MQ=255)                                                           ccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACt                                   >  1:2503235/1‑100 (MQ=255)                                                             gCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCa                                 <  1:2115421/100‑1 (MQ=255)                                                             gCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCa                                 >  2:1074639/1‑100 (MQ=255)                                                               cATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAAc                               <  1:2505398/100‑1 (MQ=255)                                                                 tGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                             <  1:133236/100‑1 (MQ=255)                                                                 tGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                             <  2:974789/100‑1 (MQ=255)                                                                 tGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                             <  1:364827/100‑1 (MQ=255)                                                                   cccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                           <  2:2184681/100‑1 (MQ=255)                                                                    ccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                          <  1:2023020/100‑1 (MQ=255)                                                                    ccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                          >  2:1867327/1‑100 (MQ=255)                                                                     cAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaac                         <  2:54946/100‑1 (MQ=255)                                                                      aaTGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaact                        >  1:259304/1‑100 (MQ=255)                                                                            tCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                  <  1:2074755/100‑1 (MQ=255)                                                                               cccTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                >  2:1149932/1‑99 (MQ=255)                                                                               cccTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACg               >  2:2291241/1‑100 (MQ=255)                                                                               cccTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACg               >  1:395061/1‑100 (MQ=255)                                                                               cccTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACg               <  2:1224769/100‑1 (MQ=255)                                                                                 cTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGAGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag             >  1:1724031/1‑100 (MQ=255)                                                                                     aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca         <  1:105239/100‑1 (MQ=255)                                                                                     aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca         >  2:2482969/1‑100 (MQ=255)                                                                                      ccACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa        >  1:1099318/1‑100 (MQ=255)                                                                                      ccACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa        >  1:1031283/1‑100 (MQ=255)                                                                                       cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag       >  1:307910/1‑100 (MQ=255)                                                                                       cacaACGGGTTCACCTGCGACGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag       <  2:597909/100‑1 (MQ=255)                                                                                        acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc      >  1:1347576/1‑100 (MQ=255)                                                                                        acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc      <  1:1239163/100‑1 (MQ=255)                                                                                        acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc      >  1:2158479/1‑100 (MQ=255)                                                                                         caACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca     <  2:1251987/100‑1 (MQ=255)                                                                                           aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTGCCCCACTCAACGCTAACTAATTGCACGAGagcaagca     <  2:74382/98‑1 (MQ=255)                                                                                           aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATg   <  1:2197074/100‑1 (MQ=255)                                                                                            cGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGa  >  1:1535104/1‑100 (MQ=255)                                                                                                  |                                                                                              TATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGA  >  NC_000913/2777266‑2777456

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