Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,987,858	G→C	11.9%	intergenic (‑100/+15)	yecJ ← / ← azuC	DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,858	0	G	C	11.9%	173.2 / 12.8	65	intergenic (‑100/+15)	yecJ/azuC	DUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base C (3/5);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCT  >  NC_000913/1987713‑1988006                                                                                                                                                  |                                                                                                                                                     ccACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGc                                                                                                                                                  <  1:1273597/150‑1 (MQ=255)  cacaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACACCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCa                                                                                                                                                 <  1:424789/150‑1 (MQ=255)   acaTCAGAAACCAATTCCTGATCTGCATTCAGCGGCTTGGACATCGTAAACTCCTGAAAGCATTTTGAAAGCCGTTATACTAGCGACTGCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAg                                                                                                                                                <  1:24664/150‑1 (MQ=255)       cAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGttt                                                                                                                                            <  1:925846/150‑1 (MQ=255)          aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCTCCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGtttttt                                                                                                                                         <  1:860084/150‑1 (MQ=255)          aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGtttttt                                                                                                                                         <  1:1027502/150‑1 (MQ=255)           acccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTAATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTg                                                                                                                                        <  1:668383/148‑1 (MQ=255)            aCCAATTCCTGATCGGCACTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCATGACACAAAAAAACCTGCTCCGGCAGGtttttttg                                                                                                                                       <  1:659288/150‑3 (MQ=255)               aaTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCAACCGTCATCAGGACACAAAAAAACCTGCCGGAGCTGCTTTTTTGTTAt                                                                                                                                    >  2:1086951/1‑150 (MQ=255)                 ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCg                                                                                                                                  <  2:745614/150‑1 (MQ=255)                 ttCCTGATCGGCATGCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGtttttttgtgtcc                                                                                                                                  <  1:668349/150‑8 (MQ=255)                  tCCTTATCTGCATTCATCGGCTGGTATATCGTAACCTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGGCACATCCACCGTCATCAGGACACAAAAAAACCTTCCGGAGCAGGTTTTTTGTTATCgg                                                                                                                                 <  2:61759/150‑1 (MQ=255)                   ccTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTAACATCCACCTACATCAGGACACAAAAAAACCTGCCGGATCAGGTTTTTAGTGATCGGa                                                                                                                                >  2:556691/1‑150 (MQ=255)                     tGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTGGAGCAGGTTTTTTGTTATCGGAAc                                                                                                                              >  1:218636/1‑150 (MQ=255)                     tGATCGGCAATCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAAc                                                                                                                              <  1:1049477/150‑1 (MQ=255)                      gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACa                                                                                                                             >  1:1262830/1‑150 (MQ=255)                      gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACa                                                                                                                             <  2:28057/150‑1 (MQ=255)                       aTCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACat                                                                                                                            >  1:546033/1‑150 (MQ=255)                               tCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTg                                                                                                                    <  1:396198/150‑1 (MQ=255)                                cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCCCCGTCATCAGGACACAAAAAAACCTGCCGGCGCAGGTGTTTTGTTATCGGAACATATTGCCTgg                                                                                                                   >  2:588891/1‑150 (MQ=255)                                cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTgg                                                                                                                   >  1:590671/1‑150 (MQ=255)                                  gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCg                                                                                                                 <  2:289053/150‑1 (MQ=255)                                      cTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGCCTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                             >  1:1108907/1‑150 (MQ=255)                                         ggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTc                                                                                                          <  2:894478/150‑1 (MQ=255)                                           aCATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                                        >  1:601324/1‑150 (MQ=255)                                            cATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                                       <  1:245946/150‑1 (MQ=255)                                             aTCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTg                                                                                                      >  2:1054444/1‑150 (MQ=255)                                                  aaaCTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGt                                                                                                 >  2:971529/1‑150 (MQ=255)                                                      tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGTTACGTCTTTGAACGTCTTg                                                                                             >  1:418684/1‑150 (MQ=255)                                                       ccTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                            >  1:1107726/1‑150 (MQ=255)                                                       ccTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                            <  2:998513/150‑1 (MQ=255)                                                        cTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCa                                                                                           <  1:10736/150‑1 (MQ=255)                                                             aGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGt                                                                                      <  1:1044658/150‑1 (MQ=255)                                                                  tttGATAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGTAACATATTGCATGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTg                                                                                 <  1:853378/150‑1 (MQ=255)                                                                    tGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGaa                                                                               <  1:746682/150‑1 (MQ=255)                                                                      aaaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATTGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACa                                                                             <  1:835551/150‑1 (MQ=255)                                                                       aatccGTTATAGTAGCTACTTCATATCTTCAGCGATAGTCACATGCTCCGTCATAAGGACACAAATAAACCTTCCTGAGCAGGTTTTTTGTTATCGGAACATAATGCGTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACAt                                                                            <  1:179342/147‑1 (MQ=255)                                                                        aGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATa                                                                           >  2:291705/1‑150 (MQ=255)                                                                           cGTTATAGTAGCGAATTCACATATTCAGAGATAGTCACATCCACCGTAATCAGTACACAAAAAAACCTGCCTGAGCAGGTTTTTTGTTATCGTAACATATTGCCTGGCGGTACGTCTTTGATCGTCTTGCAATAGTTATTGAACATACtt                                                                        <  1:400402/150‑1 (MQ=255)                                                                               ataGTAGCGACTTCACATCTTCATCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                    <  1:724060/150‑1 (MQ=255)                                                                               ataGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTTCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                    >  2:855422/1‑150 (MQ=255)                                                                                  gtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                                 >  2:821981/1‑150 (MQ=255)                                                                                       gACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTg                                                            >  2:103018/1‑150 (MQ=255)                                                                                          ttCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTgcgc                                                         <  1:971529/150‑1 (MQ=255)                                                                                          ttCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTgcgc                                                         <  1:399062/150‑1 (MQ=255)                                                                                                tCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTc                                                   <  2:509997/150‑1 (MQ=255)                                                                                                              cacaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTAACGGAACATATTCCCTGGCGGTACGTCTTTGATCGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTccc                                     >  2:93645/1‑149 (MQ=255)                                                                                                                caTCCACCGTCATCAGGACCCAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGGACGTCGTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                   <  1:824484/150‑1 (MQ=255)                                                                                                                caTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                   <  1:884723/150‑1 (MQ=255)                                                                                                                caTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGATATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTAATCGCGGCACTCCGAc                                   >  2:496907/1‑150 (MQ=255)                                                                                                                   ccACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAt                                >  1:412835/1‑150 (MQ=255)                                                                                                                    cACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAtt                               >  2:1032222/1‑150 (MQ=255)                                                                                                                            tcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTCCGGCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATAc                       <  1:1108323/150‑1 (MQ=255)                                                                                                                               ggACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGTATTTTGCGCAGTTGCATCGCTGCACTCCGACCATTTGTTATACAgg                    <  1:777025/150‑1 (MQ=255)                                                                                                                                  cacaAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtgt                 <  2:418684/150‑1 (MQ=255)                                                                                                                                     caaaaaaCCTGCTCCGGCAGGTTTTTTGTTATGGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCTGCACTCCGACCATTTGTTATACTGGTGTTAt              >  2:202032/2‑150 (MQ=255)                                                                                                                                     caaaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCATTAGTTATTGAACATACTTTTCAGGATTTTGCGCGGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt              >  1:31514/2‑150 (MQ=255)                                                                                                                                      aaaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt             >  2:910430/1‑150 (MQ=255)                                                                                                                                      aaaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCCCTCCGACCATTTGTTATACAGGTGTTAtt             >  1:1189832/1‑150 (MQ=255)                                                                                                                                      aaaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt             >  2:884006/1‑150 (MQ=255)                                                                                                                                       aaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg            >  1:577992/1‑150 (MQ=255)                                                                                                                                         aaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTc          >  1:182117/1‑150 (MQ=255)                                                                                                                                              tGCCGGAGCAGGTTTTTTGTTATTGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGCCCATTTGTTATACAGGTGTTTTTGTCTTTgc     >  2:342222/1‑150 (MQ=255)                                                                                                                                               gCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTgcg    >  2:430627/1‑150 (MQ=255)                                                                                                                                                 cGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCt  >  1:920593/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCT  >  NC_000913/1987713‑1988006

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