Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,541,448	A→C	N348H (AAC→CAC)	feoB →	ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,541,448	0	A	C	100.0%	234.9 / NA	64	N348H (AAC→CAC)	feoB	ferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (30/34);  total (30/34)

ATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGC  >  NC_000913/3541304‑3541595                                                                                                                                                 |                                                                                                                                                    aTCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCacac                                                                                                                                                 <  2:327660/149‑1 (MQ=255)         cggcggGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGcc                                                                                                                                        >  1:115198/1‑150 (MQ=255)          ggcggGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCAc                                                                                                                                      >  2:321534/1‑151 (MQ=255)              gggCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGt                                                                                                                                  <  1:701667/151‑1 (MQ=255)               ggCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTg                                                                                                                                 <  1:257439/151‑1 (MQ=255)               ggCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTg                                                                                                                                 >  2:334140/1‑151 (MQ=255)                 cGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGcc                                                                                                                               <  2:936432/151‑1 (MQ=255)                 cGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGcc                                                                                                                               <  2:326176/151‑1 (MQ=255)                    tACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGc                                                                                                                             <  2:361288/150‑1 (MQ=255)                     aCAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCa                                                                                                                            <  2:447899/150‑1 (MQ=255)                      cAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGa                                                                                                                          <  2:454925/151‑1 (MQ=255)                        gCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGAtt                                                                                                                        <  1:714749/151‑1 (MQ=255)                         ccGCTGTTTGACGGCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTg                                                                                                                       >  1:260096/1‑151 (MQ=255)                              gTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCat                                                                                                                   <  1:186269/150‑1 (MQ=255)                                 tGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCatgatg                                                                                                               >  1:888504/1‑151 (MQ=255)                                  gACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGt                                                                                                              >  1:879872/1‑151 (MQ=255)                                     gTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTAcc                                                                                                           >  1:876985/1‑151 (MQ=255)                                      tCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCt                                                                                                          >  2:251649/1‑151 (MQ=255)                                       cGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTg                                                                                                         <  1:364791/151‑1 (MQ=255)                                            ccGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTcc                                                                                                     <  1:791472/150‑1 (MQ=255)                                            ccGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTc                                                                                                      <  2:673646/149‑1 (MQ=255)                                                ggCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCtcc                                                                                                >  2:67807/1‑151 (MQ=255)                                                 gcgcTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCtcct                                                                                               >  2:807482/1‑151 (MQ=255)                                                        ttGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTg                                                                                         <  2:200139/150‑1 (MQ=255)                                                        ttGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTg                                                                                         <  2:870311/150‑1 (MQ=255)                                                         tgtgCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGa                                                                                        <  1:353015/150‑1 (MQ=255)                                                                 gTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCgg                                                                               <  2:106794/151‑1 (MQ=255)                                                                   aTTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGt                                                                             <  2:846/151‑1 (MQ=255)                                                                    ttCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGt                                                                             >  2:210703/1‑150 (MQ=255)                                                                        aTGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATg                                                                        <  2:824777/151‑1 (MQ=255)                                                                        aTGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATg                                                                        >  1:449991/1‑151 (MQ=255)                                                                          ggATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGgc                                                                      >  1:238854/1‑151 (MQ=255)                                                                              tGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGt                                                                  >  2:403127/1‑151 (MQ=255)                                                                               ggCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTg                                                                 <  1:816678/151‑1 (MQ=255)                                                                               ggCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTg                                                                 <  2:585851/151‑1 (MQ=255)                                                                                gCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTgc                                                                <  2:182727/151‑1 (MQ=255)                                                                                   acacGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTgcgg                                                              <  2:11620/150‑1 (MQ=255)                                                                                       gCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGttt                                                         >  1:400387/1‑151 (MQ=255)                                                                                           cACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTtgtg                                                      >  1:700395/1‑150 (MQ=255)                                                                                              ttCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATgg                                                  >  2:285297/1‑151 (MQ=255)                                                                                               tCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATgg                                                  >  1:806192/1‑150 (MQ=255)                                                                                               tCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATgg                                                  <  2:347017/150‑1 (MQ=255)                                                                                                   ggACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGt                                             >  2:379634/1‑151 (MQ=255)                                                                                                         gCTGACTATCTTCCTCGCCCAGAGCATGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATg                                       <  1:438441/151‑1 (MQ=255)                                                                                                           tGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGc                                      >  2:428966/1‑150 (MQ=255)                                                                                                            gACTATCTTCCTCGCCCAGGGGCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGGCTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCa                                     <  2:322962/150‑1 (MQ=255)                                                                                                            gACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCa                                     <  1:925542/150‑1 (MQ=255)                                                                                                             aCTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAgg                                   <  1:665859/151‑1 (MQ=255)                                                                                                                aTCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGgcgc                                >  1:710271/1‑151 (MQ=255)                                                                                                                   ttCCTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTgg                             >  2:892779/1‑151 (MQ=255)                                                                                                                     ccTCGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTggg                            <  2:680425/150‑1 (MQ=255)                                                                                                                        cGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCtt                         <  1:136868/150‑1 (MQ=255)                                                                                                                        cGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCt                          <  2:751030/149‑1 (MQ=255)                                                                                                                        cGCCCAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTg                        >  1:226440/1‑151 (MQ=255)                                                                                                                            cAGGGCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCgg                    >  1:381695/1‑151 (MQ=255)                                                                                                                              gggCCTGGGTGGCGGCATTCACTCCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCgg                    <  2:832131/149‑1 (MQ=255)                                                                                                                              gggCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCgggg                  >  1:243726/1‑151 (MQ=255)                                                                                                                               ggCCTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGa                 <  2:548444/151‑1 (MQ=255)                                                                                                                                 ccTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGaaa               <  2:850550/151‑1 (MQ=255)                                                                                                                                 ccTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGaa                >  2:70566/1‑150 (MQ=255)                                                                                                                                 ccTGGGTGGCGGCATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGaa                >  1:913340/1‑150 (MQ=255)                                                                                                                                         gcggcATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTtgt       >  2:41969/1‑151 (MQ=255)                                                                                                                                           ggcATTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGc     >  2:518987/1‑151 (MQ=255)                                                                                                                                              aTTCACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGc  >  2:573109/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    ATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGC  >  NC_000913/3541304‑3541595

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