Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,483,925	A→G	S289G (AGC→GGC)	hrpA →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,483,925	0	A	G	100.0%	354.2 / NA	93	S289G (AGC→GGC)	hrpA	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/54);  total (39/54)

CCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTC  >  NC_000913/1483777‑1484068                                                                                                                                                     |                                                                                                                                                ccACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGc                                                                                                                                               >  1:755582/1‑151 (MQ=255)  cACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGcc                                                                                                                                              >  1:970592/1‑151 (MQ=255)  cACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGcc                                                                                                                                              <  2:47958/151‑1 (MQ=255)       taataaTGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGc                                                                                                                                         >  2:842908/1‑151 (MQ=255)            atGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACAt                                                                                                                                    >  2:348497/1‑151 (MQ=255)            atGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACAt                                                                                                                                    <  2:834229/151‑1 (MQ=255)             tGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACAtt                                                                                                                                   >  4:326971/1‑151 (MQ=255)               cgcCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCt                                                                                                                                 <  2:1036755/151‑1 (MQ=255)                gcCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTg                                                                                                                                >  2:1033784/1‑151 (MQ=255)                 ccGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTg                                                                                                                                >  1:784444/1‑150 (MQ=255)                 ccGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTg                                                                                                                                >  3:156233/1‑150 (MQ=255)                 ccGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGa                                                                                                                               <  2:952756/151‑1 (MQ=255)                 ccGATTATTGAAGACTCCGGTCGGACCTACCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTg                                                                                                                                <  2:1233753/150‑1 (MQ=255)                   gATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATc                                                                                                                             <  1:1116547/151‑1 (MQ=255)                    attattGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCt                                                                                                                            >  1:614487/1‑151 (MQ=255)                     ttattGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCtt                                                                                                                           >  1:298553/1‑151 (MQ=255)                      tattGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCttt                                                                                                                          >  4:210265/1‑151 (MQ=255)                        ttGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCttt                                                                                                                          <  1:160274/149‑1 (MQ=255)                         tGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTAt                                                                                                                        <  1:1089891/150‑1 (MQ=255)                          gAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCAATGGTTGAATAAGCCGATGTCACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATg                                                                                                                       <  2:386938/150‑1 (MQ=255)                           aaGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATg                                                                                                                       >  2:1274377/1‑149 (MQ=255)                           aaGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGa                                                                                                                      <  2:1262021/150‑1 (MQ=255)                                 ccGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAgcgg                                                                                                                  <  1:801560/148‑1 (MQ=255)                                 ccGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGa                                                                                                               <  3:44070/151‑1 (MQ=255)                                    gTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCg                                                                                                            <  2:209046/151‑1 (MQ=255)                                      cGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCtt                                                                                                                           <  2:853874/134‑1 (MQ=255)                                      cGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCtt                                                                                                                           >  1:854046/1‑134 (MQ=255)                                         aCCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGaa                                                                                                        >  4:425726/1‑150 (MQ=255)                                           cTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGTGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATc                                                                                                     >  2:190536/1‑151 (MQ=255)                                             aTCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCcg                                                                                                   >  2:566117/1‑151 (MQ=255)                                              tCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCcgc                                                                                                  >  1:932516/1‑151 (MQ=255)                                               ccGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCcgc                                                                                                  <  3:383742/150‑1 (MQ=255)                                                cGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGa                                                                                                <  4:216976/151‑1 (MQ=255)                                                 ggTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGAt                                                                                               <  4:422447/151‑1 (MQ=255)                                                    ggAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATAcc                                                                                            <  1:307516/151‑1 (MQ=255)                                                      aaGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATCccgc                                                                                          <  4:49560/151‑1 (MQ=255)                                                      aaGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATAccgc                                                                                          >  2:1217951/1‑151 (MQ=255)                                                      aaGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATAccgc                                                                                          <  2:521401/151‑1 (MQ=255)                                                      aaGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATAccgc                                                                                          <  4:379502/151‑1 (MQ=255)                                                       aGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATAccgcc                                                                                         >  1:197452/1‑151 (MQ=255)                                                         tGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATAccgcc                                                                                         <  2:1031780/149‑1 (MQ=255)                                                            gcTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATgc                                                                                    >  1:118507/1‑151 (MQ=255)                                                            gcTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATgc                                                                                    >  1:733155/1‑151 (MQ=255)                                                                       aTTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGc                                                                         >  2:297978/1‑151 (MQ=255)                                                                        ttgttgAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCt                                                                        >  2:1091898/1‑151 (MQ=255)                                                                                 aagCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACg                                                               >  1:543394/1‑151 (MQ=255)                                                                                 aagCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACGGAGTCAGGAAGGCCATGGCGACTTTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTAc                                                                <  2:575696/150‑1 (MQ=255)                                                                                  agCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGc                                                              >  2:406912/1‑151 (MQ=255)                                                                                   gCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGc                                                              <  1:400883/150‑1 (MQ=255)                                                                                    ccGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGcc                                                             >  1:315002/1‑150 (MQ=255)                                                                                    ccGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGcc                                                             >  3:203991/1‑150 (MQ=255)                                                                                    ccGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCa                                                            <  3:281808/151‑1 (MQ=255)                                                                                      gATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCAt                                                           <  4:443247/150‑1 (MQ=255)                                                                                        tGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATAcc                                                        >  2:705301/1‑151 (MQ=255)                                                                                         gACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCg                                                       >  3:14420/1‑151 (MQ=255)                                                                                           cacCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGa                                                      <  3:281115/150‑1 (MQ=255)                                                                                            acCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGaa                                                     >  3:440292/1‑150 (MQ=255)                                                                                             ccGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGaaa                                                    <  2:400160/150‑1 (MQ=255)                                                                                             ccGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAAt                                                   <  4:373861/151‑1 (MQ=255)                                                                                              cGAGCGCGATGAGTTGCAGGCGATTTTGGTCACCGTGGACGACCGGAGTCGGGAAGGCAAGGGCGACATTCGATTGTTTATGGGGGGGAAGCGGGAAATCGGCGTTACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATc                                                  <  2:899266/151‑1 (MQ=255)                                                                                              cGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATc                                                  >  1:456229/1‑151 (MQ=255)                                                                                               gAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCt                                                 >  1:1242764/1‑151 (MQ=255)                                                                                                 gcgcGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCtt                                                <  2:41068/150‑1 (MQ=255)                                                                                                      aTCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCt                                          <  1:232676/151‑1 (MQ=255)                                                                                                      aTCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCt                                          <  1:331195/151‑1 (MQ=255)                                                                                                       tCAGTTGCAGGCGATTTTTGACGCCGAAGACAAACTGGGTCAGGAAGGCCATGGCGACATTCTGATCTTTATAAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCt                                          <  4:45069/150‑1 (MQ=255)                                                                                                        cAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCttt                                        >  2:681958/1‑151 (MQ=255)                                                                                                          gTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTAt                                      <  3:53704/151‑1 (MQ=255)                                                                                                           ttGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATg                                     >  2:423979/1‑151 (MQ=255)                                                                                                            tGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATg                                     >  4:27835/1‑150 (MQ=255)                                                                                                               aGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGcgcg                                 >  1:475884/1‑151 (MQ=255)                                                                                                                ggCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCgg                                <  1:609137/151‑1 (MQ=255)                                                                                                                  cGATTTTTGACGCCGTGGACGAACGGAGCCAGGAAGGGCAGGGCGTCATTCTGATCTTTATGAGCGGCGGGGGGAAAACCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCt                              <  2:369332/151‑1 (MQ=255)                                                                                                                    aTTTTTGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                            <  4:242852/151‑1 (MQ=255)                                                                                                                     tttttGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                            <  4:221610/150‑1 (MQ=255)                                                                                                                     tttttGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCtt                             <  1:799980/149‑1 (MQ=255)                                                                                                                      ttttGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTc                           <  1:223097/150‑1 (MQ=255)                                                                                                                      ttttGACGCCGTAGACGAACGGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTc                           <  2:246173/150‑1 (MQ=255)                                                                                                                        ttGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGGGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGa                         <  2:784285/150‑1 (MQ=255)                                                                                                                         tGACGCCGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAAc                       >  1:493256/1‑151 (MQ=255)                                                                                                                              ccGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCg                   <  1:297106/150‑1 (MQ=255)                                                                                                                              ccGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCg                   <  2:362149/150‑1 (MQ=255)                                                                                                                              ccGTAGACGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCg                   <  1:459273/150‑1 (MQ=255)                                                                                                                              ccGTACACGAACTGAGTCGGGAAGGCAATGGCGACATTCTGACCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGa                  <  2:349943/151‑1 (MQ=255)                                                                                                                               cGTAGACGACCTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGaa                 <  1:681631/151‑1 (MQ=255)                                                                                                                                    aCGAACTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACaaa             >  2:660908/1‑150 (MQ=255)                                                                                                                                        aCTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAgg        >  4:277999/1‑151 (MQ=255)                                                                                                                                         cTGAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAgg        <  3:110681/150‑1 (MQ=255)                                                                                                                                           gAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAggg       <  4:167965/149‑1 (MQ=255)                                                                                                                                           gAGTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTa     <  2:1265746/151‑1 (MQ=255)                                                                                                                                             gTCAGGAGGGCCATGGCGACATTCTGCTCTTTATGAGCGGCGAGCGGAAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTAtt   <  2:755429/151‑1 (MQ=255)                                                                                                                                             gTCAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTAtt   <  3:14330/151‑1 (MQ=255)                                                                                                                                               cAGGAAGGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTc  <  4:28628/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCACTTTAATAATGCGCCGATTATTGAAGTCTCCGGTCGGACCTATCCGGTGGAAGTGCGCTATCGCCCGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTC  >  NC_000913/1483777‑1484068

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