Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,706	G→A	V453I (GTT→ATT)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,186,706	0	G	A	100.0%	295.1 / NA	77	V453I (GTT→ATT)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/40);  total (37/40)

AAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTT  >  NC_000913/4186560‑4186852                                                                                                                                                   |                                                                                                                                                   aagaagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGAttt                                                                                                                                                 <  2:264147/150‑1 (MQ=255)          tCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCat                                                                                                                                        >  1:29363/1‑150 (MQ=255)           cGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCatat                                                                                                                                      >  1:373402/1‑151 (MQ=255)            gTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCatat                                                                                                                                      >  2:1129017/1‑150 (MQ=255)            gTTTGGGATATACTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCAtata                                                                                                                                     <  2:186692/151‑1 (MQ=255)               tGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAAc                                                                                                                                   <  1:53652/150‑1 (MQ=255)                 ggATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGcc                                                                                                                                <  1:17457/151‑1 (MQ=255)                     atCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACt                                                                                                                            <  2:933501/151‑1 (MQ=255)                          ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCg                                                                                                                         <  2:620531/149‑1 (MQ=255)                          ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCg                                                                                                                         <  1:229223/149‑1 (MQ=255)                          ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGAt                                                                                                                       <  2:1151611/151‑1 (MQ=255)                          ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGAt                                                                                                                       <  1:611189/151‑1 (MQ=255)                          ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGAt                                                                                                                       <  1:813617/151‑1 (MQ=255)                            aCGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATgg                                                                                                                     >  2:186216/1‑151 (MQ=255)                               aaGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGa                                                                                                                  <  2:845477/151‑1 (MQ=255)                                aGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAc                                                                                                                 >  2:747079/1‑151 (MQ=255)                                 gTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                                                                                                <  1:178707/151‑1 (MQ=255)                                  ttATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                                                                                                >  2:1061223/1‑150 (MQ=255)                                   tATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCa                                                                                                               <  2:630547/150‑1 (MQ=255)                                   tATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCa                                                                                                               <  2:1202849/150‑1 (MQ=255)                                    aTCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                                                                                                >  2:425315/1‑148 (MQ=255)                                    aTCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGAATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                                                                                                >  2:72858/1‑148 (MQ=255)                                      ccGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATg                                                                                                           <  1:40812/151‑1 (MQ=255)                                                     aCTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTAcc                                                                                             <  1:654050/150‑1 (MQ=255)                                                     aCTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCg                                                                                            >  2:677370/1‑151 (MQ=255)                                                        gctgAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCt                                                                                          <  2:325947/150‑1 (MQ=255)                                                         ctgAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTg                                                                                         >  1:193893/1‑150 (MQ=255)                                                           gAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGAc                                                                                       <  2:1040555/150‑1 (MQ=255)                                                            aaCCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGc                                                                                     >  1:806657/1‑151 (MQ=255)                                                              ccGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTg                                                                                   <  2:652710/151‑1 (MQ=255)                                                                      cGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCa                                                                           <  2:821917/151‑1 (MQ=255)                                                                      cGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCa                                                                           <  2:778884/151‑1 (MQ=255)                                                                      cGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCa                                                                           <  2:560819/151‑1 (MQ=255)                                                                        aCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGc                                                                         >  1:931352/1‑151 (MQ=255)                                                                             gCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTGCCGCTGACGCTGGAAGCCCAGCTgg                                                                      <  2:882390/149‑1 (MQ=255)                                                                             gCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTg                                                                       <  1:125455/148‑1 (MQ=255)                                                                                ccGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAgc                                                                  <  2:424993/150‑1 (MQ=255)                                                                                       ggTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTgcgc                                                          <  1:111807/151‑1 (MQ=255)                                                                                        gTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTgcgc                                                          >  1:759593/1‑150 (MQ=255)                                                                                          aTCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCtg                                                        >  2:88868/1‑150 (MQ=255)                                                                                            ccAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCtgatg                                                     <  2:347286/151‑1 (MQ=255)                                                                                            ccAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCtgat                                                      <  2:864007/150‑1 (MQ=255)                                                                                             cAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCtgatg                                                     >  1:686763/1‑150 (MQ=255)                                                                                              aGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTgatgat                                                   >  1:411914/1‑151 (MQ=255)                                                                                                gCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGt                                                 >  1:307702/1‑151 (MQ=255)                                                                                                  aTTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCt                                               >  1:639011/1‑151 (MQ=255)                                                                                                  aTTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCt                                               <  1:1100460/151‑1 (MQ=255)                                                                                                   tttGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCt                                               >  1:373097/1‑150 (MQ=255)                                                                                                     tGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTAcc                                            >  2:532620/1‑151 (MQ=255)                                                                                                      gAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACCCCGACTTCGATGGTGACCAGTTGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTAcc                                            >  1:109814/1‑150 (MQ=255)                                                                                                       aaCCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACca                                           >  1:1018413/1‑150 (MQ=255)                                                                                                        aCCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACcaac                                         <  2:27262/151‑1 (MQ=255)                                                                                                          cGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACcaaca                                        >  1:960378/1‑150 (MQ=255)                                                                                                            gTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCaacaac                                      >  1:976821/1‑150 (MQ=255)                                                                                                            gTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCaacaac                                      >  2:764272/1‑150 (MQ=255)                                                                                                            gTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAacaaca                                     <  2:623976/151‑1 (MQ=255)                                                                                                              aCTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATc                                   >  2:484207/1‑151 (MQ=255)                                                                                                                 gATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTg                                >  1:127783/1‑151 (MQ=255)                                                                                                                 gATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTg                                <  2:1001935/151‑1 (MQ=255)                                                                                                                  aTCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGt                               >  1:966895/1‑151 (MQ=255)                                                                                                                     gAAGGTAAAGCTATCCAGCTGCAGCCGCTGATTTGGGCGGCATATAACGCCGACTTCGATGCGGACCAGATGGCTGTTCACGTAACGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTccc                            >  1:226857/1‑151 (MQ=255)                                                                                                                     gAAGGTAAAGCTATCCAACTGCACCCGCTGATTTGTGCGGCACATAACGCCGCCTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTccc                            >  1:56457/1‑151 (MQ=255)                                                                                                                       aGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcg                          <  2:770941/151‑1 (MQ=255)                                                                                                                         gTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggc                        <  2:45509/151‑1 (MQ=255)                                                                                                                          tAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggc                        <  2:1072953/150‑1 (MQ=255)                                                                                                                           aaaGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcga                      <  1:490807/151‑1 (MQ=255)                                                                                                                           aaaGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcg                       >  2:966414/1‑150 (MQ=255)                                                                                                                            aaGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAACCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggc                        >  1:573697/1‑148 (MQ=255)                                                                                                                               cTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacgg                  >  2:432663/1‑151 (MQ=255)                                                                                                                                 aTCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacggc                 <  1:1095823/150‑1 (MQ=255)                                                                                                                                   ccAGCTGCTCCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacggcg                <  1:586986/149‑1 (MQ=255)                                                                                                                                   ccAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacggcga               <  2:744534/150‑1 (MQ=255)                                                                                                                                           aCCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAatcat      >  1:514726/1‑151 (MQ=255)                                                                                                                                            cccGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAatcatc     >  2:1147756/1‑151 (MQ=255)                                                                                                                                             ccGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCg    <  2:509126/151‑1 (MQ=255)                                                                                                                                             ccGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCg    >  1:668990/1‑151 (MQ=255)                                                                                                                                               gCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGATGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGtt  >  1:645970/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTT  >  NC_000913/4186560‑4186852

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