Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,313	(C)7→8	pseudogene (675/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,306	1	.	C	93.2%	249.7 / 12.1	75	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (31/38);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

CCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCG  >  NC_000913/1597161‑1597451                                                                                                                                                    |                                                                                                                                                  ccATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑ccccc                                                                                                                                              <  1:1029003/151‑1 (MQ=255) ccATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑ccccc                                                                                                                                              >  2:1033794/1‑151 (MQ=255)   attattGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑ccccccc                                                                                                                                            >  2:870460/1‑151 (MQ=255)      attGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑cccccccct                                                                                                                                          >  2:1122298/1‑148 (MQ=255)         ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTcgcg                                                                                                                                       <  3:238118/151‑1 (MQ=255)           aaaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGa                                                                                                                                      >  2:1069417/1‑150 (MQ=255)             atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAAt                                                                                                                                    >  2:104732/1‑150 (MQ=255)             atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAAt                                                                                                                                    >  2:90547/1‑150 (MQ=255)             atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATa                                                                                                                                   <  2:1017633/151‑1 (MQ=255)               atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGt                                                                                                                                 >  2:487447/1‑151 (MQ=255)                 tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTc                                                                                                                                <  2:228470/150‑1 (MQ=255)                 agACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCg                                                                                                                               >  2:1014658/2‑151 (MQ=255)                  gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGc                                                                                                                              <  1:516086/151‑1 (MQ=255)                     tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTc                                                                                                                           >  4:133586/1‑151 (MQ=255)                      gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTcc                                                                                                                          >  2:109652/1‑151 (MQ=255)                           aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                                                                     >  2:878635/1‑151 (MQ=255)                             tcccGAATTACCCATTATGGAGTCTCCGTACCGTGTTGACAGAGCCTTGATCTTGCACATCAAGGGTTCCAATCCCGGTGGCATTGCCACGAATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGtt                                                                                                                    <  1:711887/148‑1 (MQ=255)                             tACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGtt                                                                                                                    <  2:83550/150‑1 (MQ=255)                             tACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGtt                                                                                                                    <  1:976810/150‑1 (MQ=255)                              aCCGAATTATCCATTATAGAGTCTGCGTACCGTGATGACAGCGTCTTGATCCTGCACATTCAGGGTTCCGATACTGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGtt                                                                                                                    >  2:934396/1‑149 (MQ=255)                                cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                                                                                <  2:372049/151‑1 (MQ=255)                                  aaaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTGCCGATACCGGTGGCATTGCCACCGATAATCGTGTTTTCAGCCGTAACCA‑ATCCCCGCCTCGCGAATAGGCGCCTCCCCATTTCCTTg                                                                                                               <  2:952537/150‑1 (MQ=255)                                    atatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                                                                            >  2:789642/1‑151 (MQ=255)                                       tCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                                                                         >  4:403828/1‑151 (MQ=255)                                         cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATCCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCGCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                                                                           <  3:322810/147‑1 (MQ=255)                                            tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                                                                     <  2:240907/150‑1 (MQ=255)                                             ataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                                                    <  2:743586/150‑1 (MQ=255)                                             ataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                                                                   >  1:1217540/1‑151 (MQ=255)                                              taGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                                                                   >  2:587745/1‑150 (MQ=255)                                                 agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                                                               >  3:292077/1‑151 (MQ=255)                                                 agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                                                               <  1:878594/151‑1 (MQ=255)                                                  gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                                                                                 <  4:307291/148‑1 (MQ=255)                                                  gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                                                                              <  1:697146/151‑1 (MQ=255)                                                     tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                                                                            >  2:965976/1‑150 (MQ=255)                                                      gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                                                                           <  4:178030/150‑1 (MQ=255)                                                          aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg                                                                                        >  2:821136/1‑149 (MQ=255)                                                               gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAt                                                                                  <  4:296066/150‑1 (MQ=255)                                                               gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                                                                                 >  2:415812/1‑151 (MQ=255)                                                                  gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                                                                              >  2:1255154/1‑151 (MQ=255)                                                                   aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                                                                             <  1:960151/151‑1 (MQ=255)                                                                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                                                                            <  1:693448/151‑1 (MQ=255)                                                                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                                                                            <  2:580016/151‑1 (MQ=255)                                                                      gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttt                                                                          >  1:43308/1‑151 (MQ=255)                                                                       agTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                                         <  4:364546/151‑1 (MQ=255)                                                                       agTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                                                                        <  1:1061208/151‑1 (MQ=255)                                                                            tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAGC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                                                                      <  4:105311/149‑1 (MQ=255)                                                                            tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                                       <  2:529884/148‑1 (MQ=255)                                                                              aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                                                                  <  2:188793/151‑1 (MQ=255)                                                                                ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                                                >  1:995792/1‑151 (MQ=255)                                                                                ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                                                <  1:279999/150‑1 (MQ=255)                                                                                   gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCa                                                               <  1:848874/149‑1 (MQ=255)                                                                                   gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt                                                             <  2:823192/151‑1 (MQ=255)                                                                                     acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc                                                            >  2:674113/1‑150 (MQ=255)                                                                                        ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCCCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc                                                         >  1:389635/1‑150 (MQ=255)                                                                                          cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAcc                                                      <  1:595219/151‑1 (MQ=255)                                                                                          cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAcc                                                      <  2:312945/151‑1 (MQ=255)                                                                                                 ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTc                                                >  2:68730/1‑150 (MQ=255)                                                                                                  cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCa                                               >  2:1255232/1‑150 (MQ=255)                                                                                                      aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGa                                           <  4:10447/150‑1 (MQ=255)                                                                                                        cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACaa                                        <  4:201816/151‑1 (MQ=255)                                                                                                             gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                   <  2:472465/151‑1 (MQ=255)                                                                                                               aTTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGcc                                 >  3:42776/1‑151 (MQ=255)                                                                                                                    cACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTAc                            >  2:336279/1‑151 (MQ=255)                                                                                                                         aTAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGccc                       <  1:743090/151‑1 (MQ=255)                                                                                                                         aTAATCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGccc                       <  1:91320/151‑1 (MQ=255)                                                                                                                           aaTCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGc                     <  1:628273/151‑1 (MQ=255)                                                                                                                           aaTCGTATTTTCAGCCGTAACCA‑A‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGcc                    >  2:1171162/1‑151 (MQ=255)                                                                                                                            aTCGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGcc                    >  1:1264990/1‑151 (MQ=255)                                                                                                                              cGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCt                   <  4:54753/150‑1 (MQ=255)                                                                                                                              cGTATTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTg                  >  4:134238/1‑151 (MQ=255)                                                                                                                                 aTTTTCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATa               <  2:566226/151‑1 (MQ=255)                                                                                                                                   tttCAGCCGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGc             >  1:531189/1‑151 (MQ=255)                                                                                                                                    ttCAGCGGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGcc            >  2:227453/1‑151 (MQ=255)                                                                                                                                      cAGCCGTAACCA‑A‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTaa         <  1:1171374/151‑1 (MQ=255)                                                                                                                                         ccGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTTATAGCCTaaaa       <  2:344896/151‑1 (MQ=255)                                                                                                                                          cGTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAAt      <  1:405813/151‑1 (MQ=255)                                                                                                                                           gTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATg     >  2:7274/1‑151 (MQ=255)                                                                                                                                           gTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATg     <  1:870634/151‑1 (MQ=255)                                                                                                                                           gTAACCA‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATg     <  1:1033989/151‑1 (MQ=255)                                                                                                                                             aacca‑accCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGc   >  2:324967/1‑151 (MQ=255)                                                                                                                                              acca‑accCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcg  >  2:153059/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCG  >  NC_000913/1597161‑1597451

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