Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,724,413	A→T	5.5%	S96C (AGC→TGC)	yfiM →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,724,413	0	A	T	5.5%	314.5 / 7.1	92	S96C (AGC→TGC)	yfiM	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (40/46);  new base T (2/3);  total (42/50)
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.80e-01

ACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAG  >  NC_000913/2724265‑2724561                                                                                                                                                     |                                                                                                                                                     acaGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCaa                                                                                                                                                      >  3:656012/1‑149 (MQ=255)   aGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGc                                                                                                                                                    >  3:271249/1‑149 (MQ=255)        cAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCgg                                                                                                                                               <  2:672118/149‑1 (MQ=255)           gggATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCta                                                                                                                                            >  3:315737/1‑149 (MQ=255)             gATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTAGCCGGTGCTTGCACCGGcgaca                                                                                                                                          <  3:53913/149‑5 (MQ=255)             gATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGACCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCACGCACCGGCtata                                                                                                                                          <  4:577195/149‑1 (MQ=255)              aTGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATAc                                                                                                                                         >  1:463811/1‑149 (MQ=255)                   ccGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                    <  1:82240/149‑1 (MQ=255)                         tCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAAc                                                                                                                              <  1:198008/149‑1 (MQ=255)                            cAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGa                                                                                                                           >  3:370864/1‑149 (MQ=255)                                   aTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGTCCCGTCa                                                                                                                    >  3:586947/1‑149 (MQ=255)                                   aTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCa                                                                                                                    >  1:509091/1‑149 (MQ=255)                                   aTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCa                                                                                                                    <  4:601734/149‑1 (MQ=255)                                        tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTaaa                                                                                                               <  2:420573/149‑1 (MQ=255)                                        tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTaaa                                                                                                               <  1:72431/149‑1 (MQ=255)                                         ggATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAg                                                                                                              >  2:456569/1‑149 (MQ=255)                                          gATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGa                                                                                                             >  2:390458/1‑149 (MQ=255)                                           aTTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                            >  3:639533/1‑149 (MQ=255)                                           aTTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                            <  2:244512/149‑1 (MQ=255)                                            ttGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACg                                                                                                           >  2:390886/1‑149 (MQ=255)                                            ttGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACg                                                                                                           >  2:443843/1‑149 (MQ=255)                                            ttGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACg                                                                                                           >  4:663843/1‑149 (MQ=255)                                              gATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCa                                                                                                         >  4:327165/1‑149 (MQ=255)                                                  ttCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATccc                                                                                                     >  3:377161/1‑149 (MQ=255)                                                  ttCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGCCGCCTGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATccc                                                                                                     <  2:189592/149‑1 (MQ=255)                                                       tGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttcc                                                                                                >  2:697099/1‑149 (MQ=255)                                                       tGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCACCcccttcc                                                                                                >  1:511025/1‑149 (MQ=255)                                                        gTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttccc                                                                                               >  3:672051/1‑149 (MQ=255)                                                        gTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttccc                                                                                               >  2:534550/1‑149 (MQ=255)                                                         tCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttccct                                                                                              >  1:586774/1‑149 (MQ=255)                                                             tttGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGc                                                                                          >  1:608708/1‑149 (MQ=255)                                                             tttGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGc                                                                                          <  2:414258/149‑1 (MQ=255)                                                                gggggCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGAt                                                                                       >  2:432330/1‑149 (MQ=255)                                                                    gCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAg                                                                                   <  4:271249/149‑1 (MQ=255)                                                                       tCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCAt                                                                                >  3:268970/1‑149 (MQ=255)                                                                       tCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCAt                                                                                >  1:249119/1‑149 (MQ=255)                                                                           aatagCTGTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGAGGTCGCCGGTGCAAGCACCGGCTATACCGCCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  4:420666/146‑1 (MQ=255)                                                                           aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  3:528737/149‑1 (MQ=255)                                                                           aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  1:224904/149‑1 (MQ=255)                                                                           aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  1:390886/149‑1 (MQ=255)                                                                           aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  1:237090/149‑1 (MQ=255)                                                                           aaGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  1:429830/149‑1 (MQ=255)                                                                             gagCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            <  1:472556/147‑1 (MQ=255)                                                                             gagCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCaaa                                                                            >  2:472556/1‑147 (MQ=255)                                                                                cTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGTCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAc                                                                       <  1:588904/149‑1 (MQ=255)                                                                                 tttGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCa                                                                                                                                >  1:82715/1‑91 (MQ=255)                                                                                 tttGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCa                                                                                                                                <  2:82715/91‑1 (MQ=255)                                                                                      gATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGaaa                                                                 <  4:654005/149‑1 (MQ=255)                                                                                      gATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGaaa                                                                 <  1:620080/149‑1 (MQ=255)                                                                                            cgccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccac                                                           <  2:242820/149‑1 (MQ=255)                                                                                            cgccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccac                                                           <  1:328921/149‑1 (MQ=255)                                                                                                 cGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACgg                                                      >  2:439400/1‑149 (MQ=255)                                                                                                    aGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTACACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCa                                                   <  4:469834/149‑1 (MQ=255)                                                                                                         gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGAt                                                                                       <  2:510502/108‑1 (MQ=255)                                                                                                         gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGAt                                                                                       >  1:510502/1‑108 (MQ=255)                                                                                                         gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTa                                                                                    <  1:238255/111‑1 (MQ=255)                                                                                                         gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTa                                                                                    >  2:238255/1‑111 (MQ=255)                                                                                                              tggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTcac                                         <  1:612317/149‑1 (MQ=255)                                                                                                               ggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTcaca                                        <  2:337026/149‑1 (MQ=255)                                                                                                               ggagctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTcaca                                        <  2:588229/149‑1 (MQ=255)                                                                                                                     ggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCag                                  >  3:269087/1‑149 (MQ=255)                                                                                                                      gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCaga                                 >  4:144230/1‑149 (MQ=255)                                                                                                                      gaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCaga                                 >  4:439664/1‑149 (MQ=255)                                                                                                                       aaGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCagaa                                <  4:315737/149‑1 (MQ=255)                                                                                                                          gATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAgaagaa                             <  1:390458/149‑1 (MQ=255)                                                                                                                          gATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAgaagaa                             <  2:188137/149‑1 (MQ=255)                                                                                                                           aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                            <  4:337036/149‑1 (MQ=255)                                                                                                                           aTTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                            <  1:434982/149‑1 (MQ=255)                                                                                                                            tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccac                                                           <  3:331747/117‑1 (MQ=255)                                                                                                                            tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccac                                                           >  4:331747/1‑117 (MQ=255)                                                                                                                            tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAgg                           >  2:601495/1‑149 (MQ=255)                                                                                                                            tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAgg                           >  4:602664/1‑149 (MQ=255)                                                                                                                            tttGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAgg                           >  2:41430/1‑149 (MQ=255)                                                                                                                               ggCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTg                        <  1:41430/149‑1 (MQ=255)                                                                                                                                gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGt                       <  4:15838/149‑1 (MQ=255)                                                                                                                                  cTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGttt                     >  1:301117/1‑149 (MQ=255)                                                                                                                                    gggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTcc                   <  3:22019/149‑1 (MQ=255)                                                                                                                                     ggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCa                  <  4:673149/149‑1 (MQ=255)                                                                                                                                      gATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCAt                 <  4:322955/149‑1 (MQ=255)                                                                                                                                       aTGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCAtt                >  3:371716/1‑149 (MQ=255)                                                                                                                                       aTGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCAtt                <  1:66417/149‑1 (MQ=255)                                                                                                                                           cGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTa            <  1:34757/149‑1 (MQ=255)                                                                                                                                           agCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTa            >  4:55777/2‑149 (MQ=255)                                                                                                                                           agCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTa            >  1:702812/2‑149 (MQ=255)                                                                                                                                            gCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt           <  2:702812/149‑1 (MQ=255)                                                                                                                                            gCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt           >  3:233862/1‑149 (MQ=255)                                                                                                                                            gCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt           <  2:592130/149‑1 (MQ=255)                                                                                                                                            gCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGA‑CCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAtt          <  1:530534/149‑1 (MQ=255)                                                                                                                                                gTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTCGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAt       <  3:464019/149‑1 (MQ=255)                                                                                                                                                  gCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAATCGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAttt     <  4:377161/149‑1 (MQ=255)                                                                                                                                                  gCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGa        >  1:401352/1‑146 (MQ=255)                                                                                                                                                  gCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGa        <  2:401352/146‑1 (MQ=255)                                                                                                                                                  gCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGAttt     >  2:132238/1‑149 (MQ=255)                                                                                                                                                    aaGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCa   >  1:533856/1‑149 (MQ=255)                                                                                                                                                     aGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAg  >  2:319297/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCGCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAG  >  NC_000913/2724265‑2724561

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