Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,194,516	A→G	K551K (AAA→AAG)	yqiK →	PHB family membrane protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,194,516	0	A	G	100.0%	259.4 / NA	71	K551K (AAA→AAG)	yqiK	PHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/31);  total (40/31)

ACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAA  >  NC_000913/3194368‑3194662                                                                                                                                                     |                                                                                                                                                   aCCTGGCTGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCTTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTATCCTCAACAACTCCCGTCGAAGAAAAGGc                                                                                                                                                  <  2:968461/151‑1 (MQ=255) aCCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGc                                                                                                                                                  <  2:432669/151‑1 (MQ=255) aCCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGc                                                                                                                                                  <  1:927873/151‑1 (MQ=255) aCCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGc                                                                                                                                                  <  1:716175/151‑1 (MQ=255)            cAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAccc                                                                                                                                       <  2:945622/151‑1 (MQ=255)            cAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAccc                                                                                                                                       <  1:92349/151‑1 (MQ=255)            cAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAccc                                                                                                                                       <  1:84317/151‑1 (MQ=255)            cAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAccc                                                                                                                                       <  1:720720/151‑1 (MQ=255)            cAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAccc                                                                                                                                       <  1:596788/151‑1 (MQ=255)              aGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                                                                      >  2:1177040/1‑150 (MQ=255)              aGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                                                                      >  1:793390/1‑150 (MQ=255)              aGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATCGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAAATCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                                                                      >  1:193553/1‑150 (MQ=255)                cATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCctc                                                                                                                                     <  2:215084/149‑1 (MQ=255)                     tCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                                                                              <  2:270250/151‑1 (MQ=255)                        gCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaat                                                                                                                           >  2:203030/1‑151 (MQ=255)                         ccGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTCGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGa                                                                                                                             <  2:419196/148‑1 (MQ=255)                              ctctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGCAAAGGCAGAGTAAACCCTCTTCCAGGCATAATcc                                                                                                                      >  2:799556/1‑150 (MQ=255)                              ctctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                                                                      >  2:854077/1‑150 (MQ=255)                              ctctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                                                                      >  2:895271/1‑150 (MQ=255)                              ctctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                                                                      >  1:631241/1‑150 (MQ=255)                              ctctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                                                                      >  1:557490/1‑150 (MQ=255)                               tctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                                                                    <  1:780653/151‑1 (MQ=255)                                ctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCCGGAATAATcc                                                                                                                      >  2:967255/1‑148 (MQ=255)                                ctctTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                                                                      >  1:862780/1‑148 (MQ=255)                                         cacacaGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaaa                                                                                                          >  1:647241/1‑151 (MQ=255)                                         cacacaGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaaa                                                                                                          <  1:329393/151‑1 (MQ=255)                                           cacaGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCCCCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGGAAACCCTCCTCCAGGAATAATCCCCGGCCCCAAAAg                                                                                                         >  1:557849/1‑150 (MQ=255)                                                  aCCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                                                                                 <  2:1036720/151‑1 (MQ=255)                                                   ccGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtt                                                                                                >  1:898106/1‑151 (MQ=255)                                                   ccGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtt                                                                                                >  2:1194864/1‑151 (MQ=255)                                                   ccGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtt                                                                                                >  2:628473/1‑151 (MQ=255)                                                        gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGCAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                                                                           >  1:245642/1‑151 (MQ=255)                                                        gATTGACTCCTTGCTCAATGAAATTGGCGTTTAAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                                                                           >  1:768813/1‑151 (MQ=255)                                                         aTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                                                                           <  1:887609/150‑1 (MQ=255)                                                             aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCTGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATc                                                                                       <  1:256568/150‑1 (MQ=255)                                                             aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                                                                        >  2:628357/1‑149 (MQ=255)                                                             aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATc                                                                                       <  1:780070/150‑1 (MQ=255)                                                             aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATc                                                                                       <  1:622174/150‑1 (MQ=255)                                                               tCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                                                                    <  1:122203/151‑1 (MQ=255)                                                                            aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                                                                       >  1:879220/1‑151 (MQ=255)                                                                                tGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTGTTATATCTGGGTTCCGCTAAAAAGTGc                                                                   <  2:472888/151‑1 (MQ=255)                                                                                  gCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaa                                                                 >  1:262569/1‑151 (MQ=255)                                                                                     tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt                                                               <  1:1195351/150‑1 (MQ=255)                                                                                     tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg                                                              <  2:429389/151‑1 (MQ=255)                                                                                           gcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATc                                                           <  1:1236839/148‑1 (MQ=255)                                                                                              gcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt                                                      <  2:191527/150‑1 (MQ=255)                                                                                                  aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgc                                                 <  2:631138/151‑1 (MQ=255)                                                                                                    tGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcag                                               <  1:389784/151‑1 (MQ=255)                                                                                                     ggcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTACAAAGTGCAAATGCTCAGGTTGTTgcag                                               >  2:591672/1‑150 (MQ=255)                                                                                                       cggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAt                                            <  2:190297/151‑1 (MQ=255)                                                                                                         gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAttt                                          >  1:332020/1‑151 (MQ=255)                                                                                                         gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAttt                                          >  2:423572/1‑151 (MQ=255)                                                                                                              gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGt                                     >  2:1079540/1‑151 (MQ=255)                                                                                                              gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGATGTTGCAGCATTTTGCGt                                     >  2:642025/1‑151 (MQ=255)                                                                                                                     cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                                                                        <  1:711081/109‑1 (MQ=255)                                                                                                                     cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                                                                        >  2:710768/1‑109 (MQ=255)                                                                                                                     cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACgg                                >  2:454964/1‑149 (MQ=255)                                                                                                                     cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACgg                                >  1:942427/1‑149 (MQ=255)                                                                                                                     cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                               >  1:1023462/1‑150 (MQ=255)                                                                                                                     cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGGGCAATTGATCAGGTTGTTGAACCCTTTTGCGTGCCGGc                               >  1:147881/1‑150 (MQ=37)                                                                                                                       cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAACAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCg                            >  2:771709/1‑151 (MQ=255)                                                                                                                             cTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                      >  1:654246/1‑151 (MQ=255)                                                                                                                               caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGaatg                     >  1:243037/1‑148 (MQ=255)                                                                                                                               caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAgg                     >  2:496237/1‑150 (MQ=255)                                                                                                                               caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGc                    <  1:758015/151‑1 (MQ=255)                                                                                                                                     cTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaaa                >  1:754293/1‑149 (MQ=255)                                                                                                                                          gTCGCAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGCTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAATTGCAt          >  2:82855/1‑150 (MQ=255)                                                                                                                                          gTCGAAGAAAAGGCAGAGAAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCAt          >  2:960703/1‑150 (MQ=255)                                                                                                                                              aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag     >  1:311900/1‑151 (MQ=255)                                                                                                                                              aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag     <  2:738241/151‑1 (MQ=255)                                                                                                                                              aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag     >  2:268736/1‑151 (MQ=255)                                                                                                                                                 aaaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaaa  >  1:38770/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAA  >  NC_000913/3194368‑3194662

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