BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,906,635	C→T	V145M (GTG→ATG)	pstA ←	high‑affinity phosphate‑specific transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,906,635	0	C	T	100.0%	924.6 / NA	247	V145M (GTG→ATG)	pstA	high‑affinity phosphate‑specific transport system
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (109/138); total (109/138)

Only 100 of 257 total aligned reads displayed.

TGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGG  >  NC_000913/3906545‑3906724

 gTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCcatcac                                                                                        <  2:2514021/93‑1 (MQ=255)

  ttCTCGGTAGTGCGGATAAAAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATcaca                                                                                       >  1:6564931/1‑93 (MQ=255)

    ctcGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAAt                                                                                     <  2:4328937/93‑1 (MQ=255)

     tcGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATg                                                                                    <  1:5268154/93‑1 (MQ=255)

     tcGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATg                                                                                    <  1:2742412/93‑1 (MQ=255)

      cGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATgg                                                                                   <  2:5816298/93‑1 (MQ=255)

      cGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATgg                                                                                   <  2:5214420/93‑1 (MQ=255)

         tggtgCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGgtgt                                                                                >  2:4673161/1‑93 (MQ=255)

         tggtgCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGgtgt                                                                                >  2:5346645/1‑93 (MQ=255)

          ggtgCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTa                                                                               >  1:6403546/1‑93 (MQ=255)

          ggtgCGAATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTa                                                                               <  1:4729389/93‑1 (MQ=255)

             gCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTaaac                                                                            <  2:1734457/93‑1 (MQ=255)

             gCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTaaac                                                                            <  1:704167/93‑1 (MQ=255)

                gATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTaaacaaa                                                                         <  2:5394131/93‑1 (MQ=255)

                 aTAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTaaacaaac                                                                        >  1:5823037/1‑93 (MQ=255)

                    aCAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGa                                                                     >  1:3865885/1‑93 (MQ=255)

                    aCAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGa                                                                     >  1:1858349/1‑93 (MQ=255)

                     cAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGac                                                                    <  2:383255/93‑1 (MQ=255)

                     cAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGac                                                                    <  2:1612110/93‑1 (MQ=255)

                     cAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGac                                                                    <  2:1059997/93‑1 (MQ=255)

                      aaTCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGacc                                                                   >  2:3167629/1‑93 (MQ=255)

                        tCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGaccaa                                                                 <  2:4072978/93‑1 (MQ=255)

                        acGGCGCCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGaccag                                                                 <  1:4948309/92‑2 (MQ=255)

                        acGGCCCCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGaccaa                                                                 <  1:3496161/92‑1 (MQ=255)

                          ggCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGaccaacc                                                               <  1:1579258/93‑1 (MQ=255)

                            cACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACcac                                                             <  1:2328826/93‑1 (MQ=255)

                             aCCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACcaca                                                            >  1:58950/1‑93 (MQ=255)

                             aCCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACcaca                                                            >  2:4307488/1‑93 (MQ=255)

                             aCCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACcaca                                                            <  2:95535/93‑1 (MQ=255)

                              ccTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACaa                                                           <  2:5065734/93‑1 (MQ=255)

                              ccTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACaa                                                           >  2:3738004/1‑93 (MQ=255)

                               cTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAAt                                                          <  2:867824/93‑1 (MQ=255)

                               cTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAAt                                                          <  2:656690/93‑1 (MQ=255)

                               cTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAAt                                                          <  1:2468667/93‑1 (MQ=255)

                                 gCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCg                                                        >  2:272029/1‑93 (MQ=255)

                                 gCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCg                                                        >  1:5986398/1‑93 (MQ=255)

                                    acaacGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACg                                                     >  2:5980853/1‑93 (MQ=255)

                                      aacGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGgc                                                   >  1:4545713/1‑93 (MQ=255)

                                            aGGCCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagag                                             <  1:3594047/93‑1 (MQ=255)

                                             gggCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagaga                                            <  1:4610794/93‑1 (MQ=255)

                                             gggCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagaga                                            >  1:2685936/1‑93 (MQ=255)

                                             gggCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagaga                                            >  2:4564185/1‑93 (MQ=255)

                                             gggCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagaga                                            >  1:1623671/1‑93 (MQ=255)

                                             gggCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagaga                                            <  1:4464517/93‑1 (MQ=255)

                                                cAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCa                                         <  1:5173846/93‑1 (MQ=255)

                                                cAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCa                                         <  2:2694572/93‑1 (MQ=255)

                                                cAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCa                                         >  2:2033250/1‑93 (MQ=255)

                                                cAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCa                                         <  1:108363/93‑1 (MQ=255)

                                                  aTCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGa                                       >  2:6675975/1‑93 (MQ=255)

                                                   tCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGaa                                      >  2:1901301/1‑93 (MQ=255)

                                                    cACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAAt                                     <  2:4490516/93‑1 (MQ=255)

                                                     aCGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATg                                    >  2:1321542/1‑93 (MQ=255)

                                                      cgcccgccCAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGt                                   >  1:6311281/1‑93 (MQ=255)

                                                       gcccgcccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTc                                  <  1:339837/93‑1 (MQ=255)

                                                         ccgcccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGt                                >  2:6620087/1‑93 (MQ=255)

                                                         ccgcccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGt                                <  2:5840834/93‑1 (MQ=255)

                                                          cgcccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGtt                               <  2:3959138/93‑1 (MQ=255)

                                                            cccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTaa                             <  2:1142028/93‑1 (MQ=255)

                                                            cccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTaa                             <  2:3621411/93‑1 (MQ=255)

                                                            cccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTaa                             <  2:5326341/93‑1 (MQ=255)

                                                            cccAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTaa                             <  1:1275038/93‑1 (MQ=255)

                                                              cAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGCATGTCGTTAATg                           >  1:153145/1‑93 (MQ=255)

                                                              cAGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATg                           <  1:4175217/93‑1 (MQ=255)

                                                               aGCCGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGa                          >  2:1202457/1‑93 (MQ=255)

                                                                  cGGAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGc                       >  2:3235517/1‑93 (MQ=255)

                                                                   ggAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCg                      >  2:3277961/1‑93 (MQ=255)

                                                                   ggAGAAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCg                      >  1:6497887/1‑93 (MQ=255)

                                                                    gagaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGa                     <  1:1375407/93‑1 (MQ=255)

                                                                     agaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGaa                    <  1:303060/93‑1 (MQ=255)

                                                                     agaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGaa                    <  1:427179/93‑1 (MQ=255)

                                                                      gaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                   <  2:2473673/93‑1 (MQ=255)

                                                                      gaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                   >  2:3376847/1‑93 (MQ=255)

                                                                      gaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                   <  1:2849604/93‑1 (MQ=255)

                                                                      gaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                   >  2:831592/1‑93 (MQ=255)

                                                                      gaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                   <  2:1246859/93‑1 (MQ=255)

                                                                      gaAGTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                   >  1:4393725/1‑93 (MQ=255)

                                                                         gTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCAc                <  2:6125026/93‑1 (MQ=255)

                                                                         gTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCAc                >  1:1177296/1‑93 (MQ=255)

                                                                         gTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCAc                <  1:1001601/93‑1 (MQ=255)

                                                                         gTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCAc                <  2:4082285/93‑1 (MQ=255)

                                                                         gTGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCAc                >  1:2512899/1‑93 (MQ=255)

                                                                          tGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACt               <  1:5244160/93‑1 (MQ=255)

                                                                          tGCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACt               <  1:4048842/93‑1 (MQ=255)

                                                                           gCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACtt              <  2:2148850/93‑1 (MQ=255)

                                                                           gCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACtt              >  1:6160060/1‑93 (MQ=255)

                                                                           gCTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACtt              <  2:6485165/93‑1 (MQ=255)

                                                                            cTCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTc             <  1:2140118/93‑1 (MQ=255)

                                                                             tCCATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCt            <  1:4749599/93‑1 (MQ=255)

                                                                              ccATCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTg           >  1:3945467/1‑93 (MQ=255)

                                                                                aTCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgcc         >  1:621285/1‑93 (MQ=255)

                                                                                aTCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgcc         <  2:2905369/93‑1 (MQ=255)

                                                                                aTCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgcc         >  2:2963370/1‑93 (MQ=255)

                                                                                 tCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCCGAATGTCGTTAATGAAGCGAATCACTTCTgcca        >  1:5613545/1‑93 (MQ=255)

                                                                                 tCTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgcca        <  2:6092171/93‑1 (MQ=255)

                                                                                  cTGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgccag       <  1:644019/93‑1 (MQ=255)

                                                                                   tGCGCCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGGGCAGAGAGCAGAATGTCGTGAATGAAGCGAATCACTTCTgccagc      <  2:5953249/93‑1 (MQ=255)

                                                                                    gcgcCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgccagcc     >  2:6032009/1‑93 (MQ=255)

                                                                                    gcgcCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgccagcc     <  1:3361536/93‑1 (MQ=255)

                                                                                     cgcCATCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgccagcca    <  1:5424419/93‑1 (MQ=255)

                                                                                       ccatcaCAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAgg  >  1:2659036/1‑93 (MQ=255)

TGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGG  >  NC_000913/3906545‑3906724

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`