Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	852,223	T→C	I125M (ATA→ATG)	opgE ←	OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	852,223	0	T	C	100.0%	233.5 / NA	62	I125M (ATA→ATG)	opgE	OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/31);  total (31/31)

GATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGG  >  NC_000913/852075‑852369                                                                                                                                                     |                                                                                                                                                   gATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCAt                                                                                                                                                  >  1:943708/1‑151 (MQ=255)    aGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTAc                                                                                                                                               >  1:474011/1‑151 (MQ=255)      ccAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACTCTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGc                                                                                                                                             <  2:651145/151‑1 (MQ=255)       cAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTC‑TTTTTTCGCCTGATTGCCATTACGCCg                                                                                                                                           <  1:1024262/151‑1 (MQ=255)            tCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAgcg                                                                                                                                       >  2:676949/1‑151 (MQ=255)               ttCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAgc                                                                                                                                        <  1:460119/147‑1 (MQ=255)                tCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAgcgc                                                                                                                                      <  2:856396/148‑1 (MQ=255)                 cAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCaaa                                                                                                                                   <  2:262021/150‑1 (MQ=255)                      tAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAaaataa                                                                                                                              <  1:994329/150‑1 (MQ=255)                        aaTACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAaaataaat                                                                                                                            >  1:980663/1‑150 (MQ=255)                             gCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                      <  1:1016396/151‑1 (MQ=255)                              ccACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGa                                                                                                                     <  1:51960/151‑1 (MQ=255)                                aCCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGAcc                                                                                                                   >  1:190233/1‑151 (MQ=255)                                  cccAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGAcc                                                                                                                   >  1:76011/1‑149 (MQ=255)                                   ccaggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAgg                                                                                                                >  1:261930/1‑151 (MQ=255)                                    caggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGcc                                                                                                                       <  1:536666/143‑1 (MQ=255)                                    caggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGcc                                                                                                                       >  2:536635/1‑143 (MQ=255)                                                   aCGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGaa                                                                                                <  2:399503/151‑1 (MQ=255)                                                   aCGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGa                                                                                                 >  1:575088/1‑150 (MQ=255)                                                     gggAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGaaa                                                                                               >  2:603297/1‑150 (MQ=255)                                                                   aTGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTc                                                                                >  1:755677/1‑151 (MQ=255)                                                                     gTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCg                                                                               <  1:141488/150‑1 (MQ=255)                                                                       cGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTc                                                                            <  2:825255/151‑1 (MQ=255)                                                                           acacTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                           <  1:836068/148‑1 (MQ=255)                                                                           acacTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                           >  2:836022/1‑148 (MQ=255)                                                                            cacTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGActc                                                                       >  1:5841/1‑151 (MQ=255)                                                                                gAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCtt                                                                     <  1:92031/149‑1 (MQ=255)                                                                                    aCCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCatat                                                               <  2:2646/151‑1 (MQ=255)                                                                                        gACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtggt                                                           <  2:706116/151‑1 (MQ=255)                                                                                           acaGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtggtgg                                                         >  1:917629/1‑150 (MQ=255)                                                                                            caGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtggtgg                                                         >  1:362902/1‑149 (MQ=255)                                                                                            caGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAg                                                       <  2:941726/151‑1 (MQ=255)                                                                                             aGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAg                                                       >  2:293692/1‑150 (MQ=255)                                                                                              ggCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTg                                                     <  1:62321/151‑1 (MQ=255)                                                                                                 tGTCGAAGGTATTCCAGCGCGCTTTTTCAGTCTCTTTTTTTCGCCTGATGGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCa                                                   >  1:513079/1‑150 (MQ=255)                                                                                                     taaGGTATTCCTGCGCGCTTTTTCATTATCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCTACCAGGTTGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCGTCCATATTGGTGGAGTGCAATGc                                               >  2:169503/2‑150 (MQ=255)                                                                                                     gAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGcc                                              <  1:921118/151‑1 (MQ=255)                                                                                                     gAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAAGCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGct                                              <  2:161301/151‑2 (MQ=255)                                                                                                       aGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGcc                                              >  2:851704/1‑149 (MQ=255)                                                                                                       aGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGcc                                              >  2:807982/1‑149 (MQ=255)                                                                                                         gTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGc                                           >  1:879811/1‑150 (MQ=255)                                                                                                           aTTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAAt                                        <  2:413246/151‑1 (MQ=255)                                                                                                              ccAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGt                                      >  1:385779/1‑150 (MQ=255)                                                                                                               cAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAt                                    >  2:434847/1‑151 (MQ=255)                                                                                                               cAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAt                                    >  2:609808/1‑151 (MQ=255)                                                                                                                aGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAtt                                   <  2:266413/151‑1 (MQ=255)                                                                                                                     gcTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGTAATGCCAGCAATGTATTAAAAt                              <  2:463319/151‑1 (MQ=255)                                                                                                                     gcTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAt                              <  1:584520/151‑1 (MQ=255)                                                                                                                        ttttCATTCTCTTTTTTTCGCCTGATTGCCATTGCGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAtt                             <  1:413750/149‑1 (MQ=255)                                                                                                                        ttttCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTa                            <  2:424781/150‑1 (MQ=255)                                                                                                                           tCATTCTCTTTTTTTCGCCTGCTTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGTTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGCAATTAGgcg                        >  2:1026229/1‑151 (MQ=255)                                                                                                                              ttCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATATGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcgcg                      >  2:256370/1‑150 (MQ=255)                                                                                                                                  ctTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGcca                   <  2:737991/149‑2 (MQ=255)                                                                                                                                  ctTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCg                   <  2:362876/149‑1 (MQ=255)                                                                                                                                    ttttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCCGCAATGTATTGAAATTAGgcgcgaaa                   <  2:543758/147‑4 (MQ=255)                                                                                                                                      ttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAgga              <  2:563389/150‑2 (MQ=255)                                                                                                                                      ttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGcc             >  2:898352/1‑151 (MQ=255)                                                                                                                                       tttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGcc             <  1:24416/150‑1 (MQ=255)                                                                                                                                        ttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTg           >  1:441028/1‑151 (MQ=255)                                                                                                                                            ccTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGc        >  2:85628/1‑150 (MQ=255)                                                                                                                                                gttGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCAGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCCCCGTAGGCCTGAGCaa      >  2:175503/2‑148 (MQ=255)                                                                                                                                                 ttGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAgg  >  2:620202/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGG  >  NC_000913/852075‑852369

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