Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,556,244	T→C	5.7%	D26G (GAC→GGC)	yjiC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,556,244	0	T	C	5.7%	231.6 / 3.9	71	D26G (GAC→GGC)	yjiC	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

AAGGCTCTTTACTAATATCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAA  >  NC_000913/4556101‑4556392                                                                                                                                                |                                                                                                                                                     aaGGCTCTTTACTAATATCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACt                                                                                                                                                 >  1:126736/1‑149 (MQ=255)      tctTTACTAATATCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTgg                                                                                                                                            >  1:445814/1‑149 (MQ=255)       ctTTACTAATATCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTggg                                                                                                                                           <  4:346819/149‑1 (MQ=255)             tAATATCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGtgatga                                                                                                                                     >  2:418020/1‑149 (MQ=255)               atatCAAGTTTTCTAATTTCTACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATgaag                                                                                                                                   <  2:618684/149‑1 (MQ=255)               atatCAAGTTTTCTAATTGCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATgaag                                                                                                                                   <  4:82718/149‑1 (MQ=255)                 atCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATgaagga                                                                                                                                 >  4:335306/1‑149 (MQ=255)                  tCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATgaaggaa                                                                                                                                <  1:510097/149‑1 (MQ=255)                    aaGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTgg                                                                                                                                            >  3:63850/1‑135 (MQ=255)                    aaGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTgg                                                                                                                                            <  4:63850/135‑1 (MQ=255)                    aaGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGt                                                                                                                              >  3:69356/1‑149 (MQ=255)                    aaGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGt                                                                                                                              >  2:20536/1‑149 (MQ=255)                      gTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGa                                                                                                                            <  1:356391/149‑1 (MQ=255)                            tAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCt                                                                                                                      >  1:353609/1‑149 (MQ=255)                            tAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCt                                                                                                                      >  1:403812/1‑149 (MQ=255)                             aaTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTg                                                                                                                     <  3:438718/149‑1 (MQ=255)                                 tCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGtaat                                                                                                                 <  2:326762/149‑1 (MQ=255)                                    aCAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTGACGTGTAAACTCACTGTCGACTATTGTGTGATGAAGGAAGTGAGCATCTGTaataat                                                                                                              >  2:409110/1‑149 (MQ=255)                                       aCGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTg                                                                                                           >  1:44605/1‑149 (MQ=255)                                          cTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGAct                                                                                                        >  3:432856/1‑149 (MQ=255)                                           tttttCCATCGTAATTATGGTTATTGTCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGActc                                                                                                       <  2:28378/149‑1 (MQ=255)                                                  aTCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAggg                                                                                                >  3:572922/1‑149 (MQ=255)                                                   tCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGa                                                                                                          >  1:453292/1‑138 (MQ=255)                                                   tCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGa                                                                                                          <  2:453292/138‑1 (MQ=255)                                                   tCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGc                                                                                               >  4:293855/1‑149 (MQ=255)                                                     gTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                             >  2:516447/1‑149 (MQ=255)                                                     gTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                             <  1:20536/149‑1 (MQ=255)                                                      tAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCAtt                                                                                            >  1:282853/1‑149 (MQ=255)                                                             ggTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                     <  4:172730/149‑1 (MQ=255)                                                             ggTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGGCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                     <  3:335306/149‑1 (MQ=255)                                                              gTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATg                                                                                    <  2:603461/149‑1 (MQ=255)                                                                   tGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATa                                                                               <  1:56004/149‑1 (MQ=255)                                                                   tGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATa                                                                               <  4:69356/149‑1 (MQ=255)                                                                     ccACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGa                                                                             >  3:606249/1‑149 (MQ=255)                                                                      cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGCCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                            >  3:230564/1‑117 (MQ=255)                                                                      cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGCCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                            <  4:230564/117‑1 (MQ=255)                                                                      cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGCCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                            <  4:354809/117‑1 (MQ=255)                                                                      cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGAATATTTTAATTTACGTGTAAACTCACTGCCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                            >  3:354809/1‑117 (MQ=255)                                                                       aCGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACa                                                                           >  2:545731/1‑149 (MQ=255)                                                                          aGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataa                                                                        <  1:516447/149‑1 (MQ=255)                                                                            tAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataata                                                                      >  2:598673/1‑149 (MQ=255)                                                                            tAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataata                                                                      <  2:17046/149‑1 (MQ=255)                                                                                     gAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCttt                                                             >  3:278899/1‑149 (MQ=255)                                                                                       agatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCCACTGTTAGCTGGTGAAGGAAGTGAGTATCTGTAATAATTTGACTCTGAAGGGCATTGCTTTATGGCATAGACATAATAGTTACCTTTTa                                                           <  2:445814/149‑1 (MQ=255)                                                                                          tagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGa                                                        <  2:353609/149‑1 (MQ=255)                                                                                           agatGTTTTCCGAATCCATGCGTTCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTa                                                                  <  3:507572/138‑1 (MQ=255)                                                                                           agatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTa                                                                  >  4:507572/1‑138 (MQ=255)                                                                                            gatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTa                                                           <  2:96889/144‑1 (MQ=255)                                                                                            gatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTa                                                           >  1:96889/1‑144 (MQ=255)                                                                                                  ttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCa                                                >  4:601831/1‑149 (MQ=255)                                                                                                  ttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCa                                                <  2:126736/149‑1 (MQ=255)                                                                                                      gAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGa                                            >  2:366418/1‑149 (MQ=255)                                                                                                      gAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGa                                            >  2:31394/1‑149 (MQ=255)                                                                                                       aaTCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGAt                                           <  2:44605/149‑1 (MQ=255)                                                                                                         tCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATg                                          <  1:474639/148‑1 (MQ=255)                                                                                                         tCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATg                                          >  2:474639/1‑148 (MQ=255)                                                                                                             tGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAg                                     >  2:2880/1‑149 (MQ=255)                                                                                                                   cATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAAt                               <  3:601831/149‑1 (MQ=255)                                                                                                                    atatTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATa                              <  3:142882/149‑1 (MQ=255)                                                                                                                     tatTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAttt                                                  >  4:94508/1‑128 (MQ=255)                                                                                                                     tatTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAttt                                                  <  3:94508/128‑1 (MQ=255)                                                                                                                      atTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGa                            <  1:418020/149‑1 (MQ=255)                                                                                                                         ttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaa                         <  2:282853/149‑1 (MQ=255)                                                                                                                           aaTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAt                       <  1:2880/149‑1 (MQ=255)                                                                                                                             tttACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAttt                     >  2:487373/1‑149 (MQ=255)                                                                                                                                       aaCTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAg                                                       >  4:302887/1‑105 (MQ=255)                                                                                                                                       aaCTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAg                                                       <  3:302887/105‑1 (MQ=255)                                                                                                                                            aCTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTg                <  4:435072/139‑1 (MQ=255)                                                                                                                                            aCTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTg                >  3:435072/1‑139 (MQ=255)                                                                                                                                              tGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgta              >  3:546774/1‑139 (MQ=255)                                                                                                                                              tGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgta              <  4:546774/139‑1 (MQ=255)                                                                                                                                                tCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATAAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGaa  <  1:409110/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     AAGGCTCTTTACTAATATCAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAA  >  NC_000913/4556101‑4556392

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