Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	810,116	A→G	Y258C (TAC→TGC)	bioB →	biotin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	810,116	0	A	G	98.4%	232.3 / NA	62	Y258C (TAC→TGC)	bioB	biotin synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (33/28);  minor base C (0/1);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGC  >  NC_000913/809968‑810261                                                                                                                                                     |                                                                                                                                                  cTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCg                                                                                                                                                 >  2:320188/1‑151 (MQ=255)    gCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTgc                                                                                                                                              <  2:174624/151‑1 (MQ=255)              cGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGc                                                                                                                                    >  1:338860/1‑151 (MQ=255)                 cgccgcCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGTGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCgg                                                                                                                                 >  1:392294/1‑151 (MQ=255)                  gccgccGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCgg                                                                                                                                 >  2:359360/1‑150 (MQ=255)                   ccgccgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAc                                                                                                                               <  1:50559/151‑1 (MQ=255)                   ccgccgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAc                                                                                                                               >  1:155173/1‑151 (MQ=255)                          aaaGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGc                                                                                                                         >  2:5545/1‑150 (MQ=255)                                   cAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacga                                                                                                               <  2:192109/151‑1 (MQ=255)                                   cAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacga                                                                                                               <  2:447282/151‑1 (MQ=255)                                   cAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacg                                                                                                                >  2:374199/1‑150 (MQ=255)                                       cAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacgaac                                                                                                             >  2:472187/1‑149 (MQ=255)                                        aaCATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacgaac                                                                                                             >  1:751122/1‑148 (MQ=255)                                         aCATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGAc                                                                                                         >  1:462574/1‑151 (MQ=255)                                            tGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTc                                                                                                       >  2:425321/1‑150 (MQ=255)                                             gCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACTATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAg                                                                                                     <  2:419565/151‑1 (MQ=255)                                             gCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAg                                                                                                     <  1:738773/151‑1 (MQ=255)                                                 gtgaaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGa                                                                                                 <  1:546356/151‑1 (MQ=255)                                                   gaaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAtg                                                                                               <  2:921021/151‑1 (MQ=255)                                                      ggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGc                                                                                            <  1:969968/151‑1 (MQ=255)                                                       gtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCt                                                                                           >  1:809141/1‑151 (MQ=255)                                                         gaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCttt                                                                                         >  2:931769/1‑151 (MQ=255)                                                         gaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCttt                                                                                         >  1:466941/1‑151 (MQ=255)                                                           aaGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTAt                                                                                       >  1:493874/1‑151 (MQ=255)                                                            aGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATg                                                                                      <  2:634446/151‑1 (MQ=255)                                                               cACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGc                                                                                    >  1:545902/1‑150 (MQ=255)                                                               cACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCa                                                                                   <  2:315852/151‑1 (MQ=255)                                                                aCGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAg                                                                                  <  1:396036/151‑1 (MQ=255)                                                                 cgccgcTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAgg                                                                                 >  1:659297/1‑151 (MQ=255)                                                                         gCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAAc                                                                          <  2:148294/150‑1 (MQ=255)                                                                           cGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCg                                                                       <  2:365918/151‑1 (MQ=255)                                                                                aCGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGAtttt                                                                  >  1:838527/1‑151 (MQ=255)                                                                                  gatgatGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCt                                                                <  1:553231/151‑1 (MQ=255)                                                                                   atgatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCt                                                                >  2:7140/1‑150 (MQ=255)                                                                                     gatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACg                                                             >  1:428486/1‑151 (MQ=255)                                                                                       tgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGt                                                           >  1:290374/1‑151 (MQ=255)                                                                                          cGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAGCAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTg                                                         <  2:698227/150‑1 (MQ=255)                                                                                            aTGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                                                       >  1:895162/1‑150 (MQ=255)                                                                                               ccTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAct                                                   >  1:117204/1‑151 (MQ=255)                                                                                               ccTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAct                                                   >  2:251810/1‑151 (MQ=255)                                                                                                 tttGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAActg                                                  >  1:807214/1‑150 (MQ=255)                                                                                                    gATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAActgctg                                               <  2:200967/150‑1 (MQ=255)                                                                                                      ttttATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGAcc                                            >  1:523211/1‑151 (MQ=255)                                                                                                        ttATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCAc                                          >  1:856795/1‑151 (MQ=255)                                                                                                         tATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCAc                                          <  1:610074/150‑1 (MQ=255)                                                                                                              gCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGaa                                    <  1:499053/151‑1 (MQ=255)                                                                                                                   aTTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg                               <  1:128377/151‑1 (MQ=255)                                                                                                                       cGGTCGCGCGGATCATGATGCCAACCTCTTCCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAACCCCGATTTTCTCCGGTTGCAAACTGCTGACCACGCCGAATCCGgaag                            <  2:378470/150‑1 (MQ=255)                                                                                                                        ggTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGgaagaa                          <  1:391341/151‑1 (MQ=255)                                                                                                                           cgcgcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGgaaga                           >  1:152039/1‑147 (MQ=255)                                                                                                                            gcgcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATa                      <  1:17302/151‑1 (MQ=255)                                                                                                                             cgcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATa                      >  2:128292/1‑150 (MQ=255)                                                                                                                              gcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATaaa                    <  2:155152/151‑1 (MQ=255)                                                                                                                                ggATCATGATGCCACCCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGa                  <  2:114369/151‑1 (MQ=255)                                                                                                                                    cATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTg              <  1:982387/151‑1 (MQ=255)                                                                                                                                     atgatgCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTg              <  1:320227/150‑1 (MQ=255)                                                                                                                                          gCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTg        >  1:507015/1‑151 (MQ=255)                                                                                                                                             aaCCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTc     >  2:274710/1‑151 (MQ=255)                                                                                                                                               ccTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTcc    <  1:931853/150‑1 (MQ=255)                                                                                                                                                ctctTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCTCGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCg   <  1:5548/150‑1 (MQ=255)                                                                                                                                                ctctTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGc  >  1:148747/1‑151 (MQ=255)                                                                                                                                                 tctTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGc  >  1:374805/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGC  >  NC_000913/809968‑810261

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