Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,197,340	T→C	100%	T65A (ACT→GCT)	ymfD ←	e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,197,340	0	T	C	95.6%	276.1 / 1.3	94	T65A (ACT→GCT)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (41/45);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GT‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATA  >  NC_000913/1197244‑1197434                                                                                                  |                                                                                                cATTTTTATAATATCGTTGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AAc                                                                                             <  1:130340/100‑1 (MQ=255) cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AAc                                                                                             >  1:487254/1‑100 (MQ=255)  aTTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTATAATAATTTGCTATATTTCAAGTTGACTTTTCGAGTCTAGAAAA‑GC‑AACt                                                                                            >  2:1528900/1‑100 (MQ=255)      ttATAATATCGTGTGACATAGTCAATAATTTTAGTTTTAATTTCTCTAATAATTTGCTCTCTTTCAATTTGCCTTTTCGAGTCTAGAAAA‑GT‑AACTtc                                                                                          <  2:795593/98‑1 (MQ=255)      ttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTtcat                                                                                        <  2:2407942/100‑1 (MQ=255)       tataATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTtcatc                                                                                       >  2:1629583/1‑100 (MQ=255)         taataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCaa                                                                                     >  2:81672/1‑100 (MQ=255)          aataTCGTGGGACATAGTCAATAAATTTAATTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCaaa                                                                                    >  2:465635/1‑100 (MQ=255)           atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAAt                                                                                   >  1:1266721/1‑100 (MQ=255)            tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAAtt                                                                                  <  2:2196042/100‑1 (MQ=255)             atCGTGGTACATAATCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAAttt                                                                                 <  2:105674/100‑1 (MQ=255)              tCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTa                                                                                >  1:1330647/1‑100 (MQ=255)               cGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTAc                                                                               <  1:1658894/100‑1 (MQ=255)               cGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTAc                                                                               >  1:2012143/1‑100 (MQ=255)                  gggACATAGTCAATTATTTTAGTTTTAATTCCTCTAATAATTTTCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGa                                                                            <  2:1353575/100‑1 (MQ=255)                   ggACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAt                                                                           >  2:2308789/1‑100 (MQ=255)                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGTTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAtt                                                                          <  1:1976524/100‑1 (MQ=255)                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAtt                                                                          <  2:2186870/100‑1 (MQ=255)                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAtt                                                                          >  2:1811205/1‑100 (MQ=255)                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAtt                                                                          <  1:1320235/100‑1 (MQ=255)                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAtt                                                                          <  1:110917/100‑1 (MQ=255)                    gACATAGTCAATAATTTTAGTTTTAATACCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGAtt                                                                          <  1:2240984/100‑1 (MQ=255)                      cATAGTCGATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTaa                                                                        >  1:1589901/1‑100 (MQ=255)                       aTAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GT‑AACTTCATCAAATTTACTGATTAAt                                                                       <  2:139209/100‑1 (MQ=255)                       aTAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAAt                                                                       <  2:872393/100‑1 (MQ=255)                        tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTATAGAAAA‑GC‑AACTTCATCAAATTTAATGATTAAtt                                                                      >  2:2661331/1‑100 (MQ=255)                         aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtc                                                                     >  2:1420621/1‑100 (MQ=255)                          gTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtca                                                                    >  2:1698954/1‑100 (MQ=255)                          gTCAATAATTTTAGTTTTAATTACTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtca                                                                    >  2:172788/1‑100 (MQ=255)                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtcatc                                                                  <  1:1565489/100‑1 (MQ=255)                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtcatc                                                                  <  1:2809600/100‑1 (MQ=255)                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGAATTTTCAAGTATAAAAAA‑AA‑AAAATAATAAAATTTAATGATtaattaat                                                                   >  1:1680489/1‑99 (MQ=255)                             aataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCCAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtcataa                                                                 >  1:1682103/1‑98 (MQ=255)                             aataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTGTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATtcatc                                                                  <  1:2873945/99‑1 (MQ=255)                              ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAg                                                                <  1:2001959/100‑1 (MQ=255)                               taatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAg                                                                >  2:1366835/1‑99 (MQ=255)                                  ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAata                                                            <  1:911675/100‑1 (MQ=255)                                      aGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTa                                                        >  2:2808712/1‑100 (MQ=255)                                       gTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTaa                                                       >  2:571849/1‑100 (MQ=255)                                       gTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCATTTCGAGTTTAAAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATAAAAATATCTaa                                                       >  2:1595250/1‑100 (MQ=255)                                        ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAg                                                      >  1:2782776/1‑100 (MQ=255)                                         tttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTTCTGATTAATTCATCAGAATATCTaa                                                       >  2:1180574/1‑98 (MQ=255)                                              ttCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                <  1:671484/100‑1 (MQ=255)                                               tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                               <  2:175567/100‑1 (MQ=255)                                               tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                               <  1:25470/100‑1 (MQ=255)                                                ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                              >  2:2961429/1‑100 (MQ=255)                                                 ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                             <  2:64307/100‑1 (MQ=255)                                                 ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                             >  2:2284425/1‑100 (MQ=255)                                                     aataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCa                                         <  1:278239/100‑1 (MQ=255)                                                        aatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCaaaa                                      <  1:433134/100‑1 (MQ=255)                                                        aatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATAAGAATATCTAAGTTTTCCGCAGCCaaaa                                      >  1:177188/1‑100 (MQ=255)                                                             gCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAg                                 <  1:396836/100‑1 (MQ=255)                                                              ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAg                                 >  2:1833062/1‑99 (MQ=255)                                                              ctctctTTCAAGATACCTTTTCGAATATAAAAAA‑GA‑AACTTCATCAAATTTACTGATTAATTCATAAGAATATCTAAGTTTTCCGCAGACAAAATaaa                                  >  2:2586470/1‑98 (MQ=255)                                                               tctctTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAggg                               >  2:364807/1‑100 (MQ=255)                                                                  tgatgaaGTTG‑CTTTTCTAGACTCGAAAA‑GGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                            >  1:1191085/6‑100 (MQ=25)                                                                  tgatgaaGTTG‑CTTTTCTAGACTCGAAAA‑GGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                            >  1:2040396/6‑100 (MQ=21)                                                                   gatgaaGTTG‑CTTTTCTAGACTCGAAAAGGC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                            >  2:338196/5‑99 (MQ=25)                                                                   gatgaaGTTG‑CTTTTCTAGACTCGAAAAGGC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTg                           >  2:1133384/5‑100 (MQ=25)                                                                    ttCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                          <  2:2550448/100‑1 (MQ=255)                                                                     tCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGc                         <  1:571852/100‑1 (MQ=255)                                                                     tCAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGc                         <  1:32242/100‑1 (MQ=255)                                                                      cAAGTTGCCTTTTTTTTTTTTTAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgtc                         >  2:1558150/1‑97 (MQ=255)                                                                      cAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGcc                        <  1:2271255/100‑1 (MQ=255)                                                                      cAAGTTGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGcc                        <  1:265940/100‑1 (MQ=255)                                                                         attGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                     >  2:1776016/2‑100 (MQ=255)                                                                           tGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtttt                   >  1:532744/1‑100 (MQ=255)                                                                           tGCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAACCAAAATCAAGGGCGTGGCCAtttt                   <  2:1559055/100‑1 (MQ=255)                                                                            gCCTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAttttt                  >  1:2546901/1‑100 (MQ=255)                                                                             ccTTTTCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                 <  2:1847450/100‑1 (MQ=255)                                                                               ttttCTAGTCTATAAAA‑GC‑ATTTTTATCTTATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAATTTAAGGGCGTGGTCATTTTtc                 <  2:2440540/98‑1 (MQ=255)                                                                               ttttCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc               <  1:2200907/100‑1 (MQ=255)                                                                               ttttCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc               <  1:671068/100‑1 (MQ=255)                                                                                tttCGAGTCTAGAAAA‑GT‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa              <  1:2488482/100‑1 (MQ=255)                                                                                tttCGAGTCTAGAAAA‑GC‑AACTTCATCAATTGTACTGATTAATTCATCAGAATATCTAAGTTTTCCGTAGCCAAAATTAAGGGCGTGGCCATTTGtctc               <  2:158780/99‑1 (MQ=255)                                                                                tttCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa              <  1:953254/100‑1 (MQ=255)                                                                                tttCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa              >  2:1001975/1‑100 (MQ=255)                                                                                tttCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa              <  1:2364674/100‑1 (MQ=255)                                                                                  tCGAGTCTAGAAAA‑GT‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAAt            >  2:1565741/1‑100 (MQ=255)                                                                                  tCGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAAt            <  1:35973/100‑1 (MQ=255)                                                                                   cGAGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg           <  2:2975517/100‑1 (MQ=255)                                                                                     aGTCTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAACCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt         >  1:133435/1‑100 (MQ=255)                                                                                        cTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt      <  1:1742211/100‑1 (MQ=255)                                                                                        cTAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTACGCAGCCAAAACCAAGGACGTGGCCATTTTTCTCAATgat         >  2:2555556/1‑95 (MQ=255)                                                                                         tAGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCTTCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTa     >  1:2377465/1‑100 (MQ=255)                                                                                          aGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTaa    >  2:1453719/1‑100 (MQ=255)                                                                                          aGAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTaa    >  1:2028645/1‑100 (MQ=255)                                                                                           gAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAATGGCGTGGCCATTTTTCTCAATGCTTATAat   <  2:1081174/100‑1 (MQ=255)                                                                                           gAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat   >  2:922666/1‑100 (MQ=255)                                                                                           gAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat   >  1:1885491/1‑100 (MQ=255)                                                                                           gAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat   <  1:196501/100‑1 (MQ=255)                                                                                           gAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat   <  1:2372899/100‑1 (MQ=255)                                                                                           gAAAA‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat   <  1:34071/100‑1 (MQ=255)                                                                                            aaaa‑GC‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata  <  1:1832308/100‑1 (MQ=255)                                                                                                  |                                                                                                CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑GT‑AACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATA  >  NC_000913/1197244‑1197434

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