Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,118,637	T→G	6.5%	Q137P (CAG→CCG)	yceJ ←	putative cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,118,637	0	T	G	6.5%	218.5 / 5.4	73	Q137P (CAG→CCG)	yceJ	putative cytochrome b561
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base G (3/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCG  >  NC_000913/1118488‑1118780                                                                                                                                                      |                                                                                                                                                cATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGATGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCt                                                                                                                                                 <  2:764535/150‑1 (MQ=255)  aTAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCtg                                                                                                                                                <  1:182609/150‑1 (MQ=255)      tCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTgcg                                                                                                                                            <  1:1476602/150‑1 (MQ=255)       cAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgccc                                                                                                                                           <  1:1076207/150‑3 (MQ=255)               tacTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                                                                                                   <  2:1349574/149‑1 (MQ=255)               gacTTTCCCAGCATGCGCTGCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                                                                                                   <  1:980353/150‑1 (MQ=255)                 cTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGc                                                                                                                                 >  1:1595073/1‑150 (MQ=255)                    tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGcc                                                                                                                                          <  2:558674/138‑1 (MQ=255)                    tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGcc                                                                                                                                          >  1:558674/1‑138 (MQ=255)                    tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                              >  2:550034/1‑150 (MQ=255)                    tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                              <  1:365938/150‑1 (MQ=255)                           aTGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAGCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTcgcg                                                                                                                       <  1:926509/150‑1 (MQ=255)                            tGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCgg                                                                                                                      <  1:824396/150‑1 (MQ=255)                               gcTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGa                                                                                                                   >  1:670934/1‑150 (MQ=255)                                cTTCAGAGTGTCGTCTTTATCGATGAAATGATGGTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGAc                                                                                                                  <  2:149252/150‑1 (MQ=255)                                cTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGAc                                                                                                                  <  1:274552/150‑1 (MQ=255)                                cTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGAc                                                                                                                  >  1:1419088/1‑150 (MQ=255)                                  tCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGt                                                                                                                <  2:391841/150‑1 (MQ=255)                                   cAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTc                                                                                                               <  2:484877/150‑1 (MQ=255)                                      agTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCaaa                                                                                                            <  2:334580/150‑1 (MQ=255)                                          tcgtcTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAg                                                                                                        <  2:320644/150‑1 (MQ=255)                                           cgtcTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGc                                                                                                       <  2:1406107/150‑1 (MQ=255)                                            gtcTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGcc                                                                                                      <  2:1439749/150‑1 (MQ=255)                                                ttATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCaaaa                                                                                                  <  2:1151284/150‑1 (MQ=255)                                                    cGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGc                                                                                              >  1:456246/1‑150 (MQ=255)                                                     gATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCt                                                                                             >  1:584241/1‑150 (MQ=255)                                                     gATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCt                                                                                             >  1:1300752/1‑150 (MQ=255)                                                     gATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCt                                                                                             >  2:745249/1‑150 (MQ=255)                                                     gATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCt                                                                                             >  1:984621/1‑150 (MQ=255)                                                      aTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCTAAAATGCGGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                            <  2:1612615/150‑1 (MQ=255)                                                         aaaTGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATc                                                                                         >  1:39050/1‑150 (MQ=255)                                                         aaaTGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATc                                                                                         <  1:1238227/150‑1 (MQ=255)                                                          aaTGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCg                                                                                        <  2:1236738/150‑1 (MQ=255)                                                            tgatgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGt                                                                                      <  2:581795/150‑1 (MQ=255)                                                             gatgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGtt                                                                                     >  1:585456/1‑150 (MQ=255)                                                                 cTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTAcc                                                                                 >  2:25813/1‑150 (MQ=255)                                                                 cTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTAcc                                                                                 >  2:379346/1‑150 (MQ=255)                                                                   tCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCAt                                                                               >  2:275265/1‑150 (MQ=255)                                                                     aGGGCCATAAATCCGTGCATAACGGACAGTACTACGCCGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCg                                                                             <  2:821693/150‑1 (MQ=255)                                                                            tAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTg                                                                      <  2:260661/150‑1 (MQ=255)                                                                              aaTCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAGACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTgag                                                                    <  2:476410/150‑1 (MQ=255)                                                                              aaTCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTgag                                                                    >  2:384141/1‑150 (MQ=255)                                                                                 ccGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATc                                                                 <  1:1020570/150‑1 (MQ=255)                                                                                   gTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAg                                                               <  2:123221/150‑1 (MQ=255)                                                                                       aTAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAg                                                           >  1:1247682/1‑150 (MQ=255)                                                                                       aTAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAg                                                           <  2:203851/150‑1 (MQ=255)                                                                                          aCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCTGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAcccc                                                        >  2:5817/1‑149 (MQ=255)                                                                                             gACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTg                                                     >  1:161340/1‑150 (MQ=255)                                                                                               cAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGAt                                                   >  1:1136470/1‑150 (MQ=255)                                                                                               cAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGAt                                                   >  2:1002908/1‑150 (MQ=255)                                                                                                     tacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGAt                                             <  1:199640/150‑1 (MQ=255)                                                                                                      acgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTCACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATg                                            <  2:670934/150‑1 (MQ=255)                                                                                                       cgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATgg                                           >  1:95122/1‑150 (MQ=255)                                                                                                        gacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGc                                          <  2:538727/150‑1 (MQ=255)                                                                                                           gCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGaa                                       <  1:550034/150‑1 (MQ=255)                                                                                                             tCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGaaaa                                     >  1:28449/1‑150 (MQ=255)                                                                                                             tCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGaaaa                                     <  1:745249/150‑1 (MQ=255)                                                                                                              ccACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGCTGCCGATGGCGAAAAg                                    >  1:1371330/1‑150 (MQ=255)                                                                                                              ccACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAg                                    >  2:414635/1‑150 (MQ=255)                                                                                                               cACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGc                                   <  1:675116/150‑1 (MQ=255)                                                                                                                  gcgAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAAt                                <  2:547983/150‑1 (MQ=255)                                                                                                                    gAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAg                              <  2:1020084/150‑1 (MQ=255)                                                                                                                       ccAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTa                           >  1:58939/1‑150 (MQ=255)                                                                                                                       ccAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTa                           >  1:1268952/1‑150 (MQ=255)                                                                                                                       ccAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTa                           >  2:1209170/1‑150 (MQ=255)                                                                                                                             aTGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGc                     >  1:1202944/1‑150 (MQ=255)                                                                                                                              tGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCg                    >  1:669350/1‑150 (MQ=255)                                                                                                                               ctggcgCACAGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga                   >  1:906384/6‑150 (MQ=255)                                                                                                                                   ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCACTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATg               >  2:1256408/1‑150 (MQ=255)                                                                                                                                    gCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACCTCAAACCAGCCAAAAACGCTGATCGGTTTACAATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGAGAAAAGCAATAGGTACAATGCGAGATGt              >  2:559355/1‑150 (MQ=255)                                                                                                                                     cacaGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGGTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATCGGTACCATGCGCGATGtg             >  2:1059694/1‑149 (MQ=255)                                                                                                                                     cacaGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGAAAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTc             >  2:202690/1‑150 (MQ=255)                                                                                                                                         ggCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGc         >  1:784442/1‑150 (MQ=255)                                                                                                                                           cAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGctc       <  1:278073/150‑1 (MQ=255)                                                                                                                                                tCTGCCGGTGCGCCAGCGTGGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCg  <  2:550760/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                CATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCG  >  NC_000913/1118488‑1118780

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